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Abstract 



Evidence for pair production of a new particle consistent with the Standard Model 
top quark has been reported recently by groups studying proton anti-proton col- 
lisions at 1.8 TeV center of mass energy at the Fermi National Accelerator Lab- 
oratory. In this paper we review the history of the search for the top quark in 
electron positron and proton anti-proton collisions. We report on a number of pre- 
cise electro- weak measurements and the value of the top quark mass which can be 
extracted from these measurements within the context of the Standard Model. We 
review the theoretical predictions for top quark production and the dominant back- 
grounds. We describe the collider and the detectors that were used to measure the 
pair production process and follow on to describe the data from which the existence 
of the top quark is evinced. Finally, we present possible measurements that could be 
made in the future with more data, measurements of quantities that would confirm 
the nature of this particle, the details of its production in hadron collisions, and its 
decay properties. 
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I. INTRODUCTION 

The intensive experimental efforts in the search for the heaviest fundamental fermion 
culminated in 1995 with the discovery of the top quark in proton-antiproton annhilations 
at the Fermilab Tevatron Collider. Observation of the top quark is the latest in a long 
series of triumphs for the Standard Model of particles and fields. The top quark is the last 
fundamental fermion and the next-to-last fundamental particle predicted by the Standard 
Model. Only the Higgs boson remains unobserved. 

The search for the top quark started in the late seventies, soon after discovery of the 
companion bottom-quark. It has been a long and arduous process because the top quark 
turned out to be much more massive than was originally expected. The mass of the top 
quark is remarkably large, approximately 200 times larger than the mass of the proton and 
40 times higher than the mass of the next lightest quark. Whether this property of the top 
quark is a mere accident or a manifestation of a deeper physical process is an unanswered 
question in particle physics. 

The discovery of the top quark has been made possible by the technological progresses 
in high energy physics in the past fifteen years. In particular, the development of proton- 
antiproton colliders, pioneered first at CERN and then at Fermilab has been a crucial in- 
gredient in the discovery of the top. 

In this article we review the discovery of the top quark as well as the developments that 
led to it. In Section [H| we discuss the top quark within the framework of the Standard 
Model. While the top quark mass is a free parameter in the Standard Model, its value 
enters in calculations of a number of electroweak observables. The top-mass dependence 
of the theoretical predictions is in general rather weak. However, the accuracy of many of 
these measurements is now such that meaningful constraints on the top quark mass can be 
obtained by comparing them with theoretical predictions. These constraints constitute a 
test of the predictive power of the Standard Model, and are reviewed in Section |ITT[ The top 
production mechanisms in proton-antiproton collisions, and the experimental signatures of 
top events that are crucial to the understanding of the experimental results are discussed in 
Section fV| and [V]. Early searches for the top quark are described in Section [Vl[ The Tevatron 



Collider, whose remarkable performance played a very important role in the discovery of 
the top quark, will be described briefly in Section [VII| . The data that finally led to the 
discovery of the top quark are reviewed in Section |V111| . The value of the top mass is 
of fundamental importance, and it is needed to complete precise tests of the Standard 
Model. Furthermore, from an experimental point of view, the techniques developed to 
measure the top mass are new and particularly interesting. The top mass measurement is 
described in Section fX[ Historically, discoveries of new leptons or quarks have opened up 
new fields of inquiry which have enhanced our understanding of elementary particles and 
their interactions. Consequently, this article concludes in Section |X| with a discussion of the 
experimental prospects for top physics. 
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II. THE ROLE OF THE TOP QUARK IN THE STANDARD MODEL 



Quarks and leptons constitute the basic building blocks of matter in the Standard Model 
(SM). There are three generations of quarks and leptons in the model, with identical quan- 
tum numbers, but different masses. Within each generation, quarks and leptons appear in 
pairs, (see Fig. [I]). The left-handed quarks form weak isospin doublets, with the Q= +2/3 
and Q= —1/3 quarks having weak isospin I 3 = +1/2 and —1/2 respectively. 
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FIG. 1. Leptons and quarks in SU(2)xU(l) (Standard Model). Also shown are the values for the 
SU(2) weak isospin (I3), U(l) weak hypercharge (Y), and electric charge (Q, in units of the electron 
charge). The subscripts L and R refer to the left and right-handed components respectively. 

The tau-lepton (r) was the first particle of the third generation to be discovered (Perl et 
al, 1975). A short time later, in 1977, the T was discovered at Fermilab (Herb et al, 1977) 
as a resonance in the fi + n~ invariant mass spectrum in the reaction p + nucleon — > /i + /i~+ 
X. This resonance was interpreted as a bb bound state (the T), which subsequently decays 
into muon pairs. As will become abundantly clear in the remainder of this paper, the top 
signature in hadron collisions is much more complicated. 

In the past fifteen years, a tremendous amount of experimental data on the properties 
of the 6-quark and of fe-flavored hadrons have become available, mostly from experiments at 
e + e~ colliders. Both the charge and the weak isospin of the bottom-quark are by now well 
established (Q& = —1/3, and I 3 = —1/2). 

The value of the charge was first inferred from measurements of the T leptonic width 
(Berger et al, 1978; Darden et al, 1978; Bienlein et al, 1978) at the DORIS e + e~ storage 
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FIG. 2. Feynman diagram for T — > e + e . 





Model 1 


Model 2 


Model 3 


Experimental value 


1.05 keV 


1.07 keV 


1.15 ± 0.20 keV 


1.22 ± 0.03 keV 





TABLE I. The experimental value of T(T — ► e + e~) compared with theoretical expectations from 
quark-antiquark potential models assuming Q& = —1/3. Model 1 : Krasemann and Ono, 1979. 
Model 2 : Buchmuller, Grunberg and Tye, 1980. Model 3 : Voloshin and Zakharov, 1980. The 
quoted experimental value is from the compilation of the Particle Data Group, corrected by 7 % for 
consistency in comparison with potential model calculations, which only include the lowest order 
Born term (Montanet et al, 1994). 



ring. This width is proportional to the square of the charge of the 5-quark, (see Fig. 2), 
and can be quantitatively estimated from heavy quark-antiquark potential models, see Ta- 
ble H (Quigg and Rosner, 1977; Eichten and Gottfried, 1977; Rosner, Quigg, and Thacker, 
1978; Krasemann and Ono, 1979; Buchmuller, Grunberg and Tye, 1980; Voloshin and Za- 
kharov, 1980). The charge assignment was subsequently confirmed by measurements of R 
= (j(e + e~ — > hadrons) /cr(e + e _ — > fi + fi~). At lowest order, and ignoring resonance effects, 
R = J2quarks ^Q 2 q i where the factor of three arises from the fact that quarks come in three 
colors. The sum is over all quarks that can be produced, i.e. all quarks with mass below 
one-half the center of mass energy of the e + e~ system. Above threshold for bb production, 
the value of R was found to increase by 0.36 ± 0.09 ± 0.03 (Rice et al, 1982), in agreement 
with the expectations of 3Qj=l/3. The pole mass of the 6-quark is estimated to be in the 
range 4.5 - 4.9 GeV/c 2 from knowledge of the T and 6-meson masses (Montanet et al, 1994). 

The weak isospin of the 6-quark was first extracted from the forward-backward asym- 
metry (Afb) in e + e~ — > bb. This asymmetry is defined in terms of the 6-quark production 
cross section a{b) as 

_ a{b,6> 90°) -a(b,6 < 90°) 
FB ~ a(b,6> 90°) + a(b,6 < 90°) 

where 9 is the polar angle of the 6-quark in the e + e~ center of mass as measured from 
the direction of flight of the e~. The asymmetry originates from the coupling of the Z to 
fermions, which in the Standard Model depends on the weak isospin through a term in the 
lagrangian of the form f^^gv — QAl^Z^f, where the vector and axial couplings gy and g a 
are given by 
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FIG. 3. An example of a fermion triangle diagram which could cause an anomaly; g{ is the fermion 
axial coupling to the Z, and Q/ is the fermion charge. 

h ~ 2Q sin 2 6 W 
® V 2 sin 9w cos 6w 
h 

2 sin 9w cos 9w 

and 9w is the Weinberg angle. The first measurement of A F b was performed in the mid- 
eighties (A FB = -22.8 ±6.0 ±2.5% at ^ = 34.6 GeV; Bartel et al, 1984), and was found to 
be consistent with the Standard Model prediction (A FB = —25%) assuming I 3 = —1/2 for 
the weak isospin of the 6-quark. Alternative isospin assignments (e.g. 1=0) for the bottom- 
quark were also found to be inconsistent with the observed suppression of flavor-changing 
neutral current decays of b- mesons. If the 6-quark formed a weak-isospin singlet, and if 
there were only five quarks (u, d, c, s, b), then it can be shown that BR(B— >X1 + 1~) > 0.12 
BR(B^Xlzz) ~ 0.026 (Kane and Peskin, 1982). This was soon found to be inconsistent with 
the first upper limits placed on flavor changing neutral currents in 6-decays, BR(B— »X1 + 1~) 
< 0.008 at 90% C.L. (Matteuzzi et al, 1983). 

The I3 = —1/2 isospin of the bottom quark implies the existence of an additional quark, 
the top quark, as the third-generation weak isospin partner of the bottom quark. Fur- 
thermore, the existence of such a third generation quark doublet, in conjunction with the 
presence of three lepton generations, ensures the necessary cancellations in diagrams con- 
tributing to triangle anomalies. For the electro-weak theory to be renormalizable, the sum 
over fermions for diagrams such as the one displayed in Fig. 3 should vanish (see for example 
Leader and Pedrazzi, 1982). The contributions to this diagram for each fermion in the the- 
ory is proportional to N c g{Qj, where the factor N c = 3 is the number of colors and applies 
to quarks only. Hence, the contribution from a lepton isodoublet exactly cancels that of a 
quark isodoublet. With three lepton generations, the existence of a third quark isodoublet, 
whose members are the top and bottom quarks, results in the desired cancellation of triangle 
anomalies. 

Measurements of the Z width at the LEP and SLC colliders rule out the existence of 
a 4th generation neutrino with mass M u < M^/2 (Montanet et al, 1994). Unless the 4th 
generation neutrino is very massive, no additional generations are allowed in the context of 
the Standard Model. The top quark is therefore the last fermion expected in the Standard 
Model. Only the Higgs boson is left to be discovered in order to complete the particle and 
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field content of the minimal Standard Model. 

While the Standard Model predicts the charge and weak isospin of the top quark (Q = 
2/3, and I3 = 1/2), its mass remains a free parameter. As we will discuss in Section [IX|, 
the recent top mass measurements yield M top = 175 ± 8 GeV/c 2 , a factor of 40 higher than 
the mass of the second heaviest fundamental fermion (the 6-quark). The reason for such 
a high mass for the top quark is a mystery of the Standard Model. It does however occur 
quite naturally in local supersymmetric theories where the electro- weak symmetry is broken 
through radiative corrections (Ibanez and Lopez, 1983; J. Ellis et ai, 1983; Alvarez- Gaume, 
Polchinski, and Wise, 1983). 

The value of the top mass enters in the calculation of radiative corrections to a large 
number of electro-weak observables. As we will discuss in Section |TTT|, the level of precision 
achieved in these measurements is good enough that a comparison between the measured 
top mass and the calculation of electro-weak radiative corrections provides a stringent test 
of electro-weak theory and is sensitive to physics beyond the Standard Model. 



III. TOP MASS AND PRECISION ELECTRO-WEAK MEASUREMENTS 

Over the past few years, a number of very precise measurements at e + e~ colliders have 
been performed using large samples of Z events. At CERN the four LEP experiments 
(Aleph, Delphi, L3, and Opal) have collected of order two million Z events each, and at 
the SLC collider at SLAC the SLD experiment has collected a data sample of order one 
hundred thousand events. Although the SLAC data sample is considerably smaller than 
the LEP data sample, it was possible to polarize the beams at the SLC and this lead to a 
competitive measurement of the Weinberg angle. The level of accuracy of the LEP and SLC 
measurements, as well as of measurements of the W^-boson mass at the Tevatron, is now 
such that the data are not well described by tree-level theoretical calculations and radiative 
corrections must be included. 

If we assume that the Standard Model can be used to correctly calculate higher order 
electro-weak processes, we can infer the top quark mass by comparing these calculations to 
precise measurements (Altarelli, Kleiss, and Verzegnassi, 1989). We can then check whether 
this value is consistent with the directly measured top quark mass. The free parameters in 
the model are the weak coupling constant Gp ermi , the electromagnetic coupling constant a, 
the Weinberg angle, 9 W , the mass of the Higgs boson, M Higgs , the strong coupling constant 
a s , the masses of the six quarks, the masses of the six leptons, and the four quark mixing 
parameters which determine the CKM matrix. 

At lowest order the masses of the weak intermediate vector bosons can be determined 
completely from the first three of these parameters, G Fermi, ®, and 6w- The best measured 
Standard Model parameters are a, G Fermi and M Z - Using these three parameters we can, at 
lowest order predict several measurable quantites. However, when higher order corrections 
are considered, the fermions and Higgs masses enter into the calculations. The top mass 
plays a particularly large role in these radiative corrections due to the large mass difference 
beween the top quark and its weak isospin partner, the bottom quark. The dependence of 
these radiative corrections on the top quark mass contains terms quadratic in M top of the 
form Mj op /M|, whereas the Higgs mass dependence is logarithmic (\n(M Higgs /M z ))- Fig. 4 
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shows two higher order diagrams involving top quark and Higgs radiative corrections which 
modify observables at the Z resonance. 



A. Neutral current experiment measurments 

At LEP and SLC, the line shape and asymmetries at the Z have been precisely mea- 
sured and can be compared with theoretical prediction. The uncertainties on most of these 
quantities are smaller from the LEP measurements due to the size of the LEP data samples. 
These quantities are: 

• The total width of the Z, T z . 

• The value of the hadronic cross section at the Z peak, a^ ad = ^ Tee R had . 

m z z 

• The ratio of the hadronic to leptonic widths, Ri = T had /Tu. 

• The forward-backward asymmetry in Z — > 11 decays, A / B = jA e Af where Af = 
2 9vf9Af anc | ^ e leptons (I) include e, a and r. 

• A T as defined above. This is obtained from measurements of the r polarization defined 
as P r = where a R and a L are 
right and left-handed r's respectively. 



as P r = a a B R+ a a L L where (Jr and are the r pair cross sections for the production of 



• A e , as defined above, i.e. A e = 2 l Ve9 i e . 

• The forward-backward asymmetry for decays, Z — > 66 and Z — > cc, at the Z pole 
mass, A°p b B and yl^ 



c 



• The value of sin 2 O^ff = |(1 — ^) from the hadronic charge asymmetry, < Qfb > 
which is the forward-backward asymmetry measured from Z — > qq decays. The charge 
of the outgoing quark is determined using a statistical weighting method. 

p&t) pec 

• The ratios Rb = rha £. ons and R c = rha f rons ■ 

1 z 1 z 

• The left-right asymmetry, A LR = (o" L — (J R )/{a L + cr^), where o" L and o"^ are the 
production cross sections for Z bosons at the Z pole energy with left-handed and 
right-handed electrons respectively. This measurement has been performed only at 
the SLC, since the beams are unpolarized at LEP. 

H 





FIG. 4. Examples of radiative corrections to the Z mass involving top quarks or Higgs boson loops. 
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In Fig. [| and Fig. |6| we show the measurements of the variables listed above at LEP (The 
LEP Collaborations, 1995), compared with Standard Model expectations as a function of 
the top quark mass. The vertical bands are the measurements, where the width of the bands 
include one standard deviation uncertainties. The cross-hatched bands show the theoretical 
predictions taking into account the uncertainty on the Higgs mass (the inner bands) and 
the uncertainty on the value of the strong coupling constant a s (the outer bands). In 
calculating the theoretical uncertainties, the Higgs mass is varied from 60 to 1000 GeV/c 2 
and a s is varied within the interval a s (Mj) = 0.123±0.006 (Bethke, 1995). The dependence 
of the Z width measurements on a s enters through radiative diagrams involving gluons, an 
example of which is shown in Fig. 7. The LEP measurements are consistent with theoretical 
predictions, except for the ratios R& and R c , which will be discussed in Section [III.A.2| . 



1. Asymmetries at the Z 

The asymmetries measured at LEP, A°p l B , A T , A e , A°p b B , A°p C B and < Qfb > are effectively 
measurements of sin 2 O^fl- Despite the lower statistics SLC Z sample, the A^r measurement 
from SLD provides a competitive measurement of sin 2 O^ff- sin 2 O^ff has also been measured 
in vN experiments at lower center of mass energies (Abramowicz et al., 1986; Blondel et 
al, 1990; Allaby et al, 1986 and 1987; Arroyo et al, 1994). All of these measurements of 
sin 2 6 l ^Jf are found to be consistent within their respective uncertainties and are combined. 
The combined result from LEP is sin 6 l e Jf = 0.23186 ± 0.00034, while the corresponding 
result from the SLD A LR measurement is sin 2 6%$ = 0.23049 ± 0.00050 (Woods, 1996). 
Figure ||] summarizes all the measurements of sin 2 9 l £ff. 
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Preliminary 




FIG. 5. Comparison of LEP measurements with Standard Model predictions as a function of Mt op - 
The experimental errors on the parameters are indicated as vertical bands. The cross-hatch pattern 
parallel to the axes indicates the variation of the Standard Model prediction with MHi ggs spanning 
the interval 60 < Mmggs < 1000 GeV/c 2 , and the diagonal cross-hatch pattern corresponds to a 
variation of a s (M§) within the interval a s (M§) = 0.123 db 0.006. The total width of the band 
corresponds to the linear sum of both uncertainties (The LEP Collaborations, 1995). 
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Preliminary 




FIG. 6. Comparison of LEP measurements with Standard Model predictions as a function of M( 0p 
as in the previous figure. For the ratios of the partial widths the variations with M.Hi ggs and 
a s (M|) nearly cancel. For the comparison of R& with the Standard Model the value R c has been 
fixed to its Standard Model prediction. To illustrate the impact of special vertex corrections to 
the Standard Model prediction for is also shown (The LEP Collaborations, 1995). 
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FIG. 7. One Feyman diagram illustrating a QCD correction to the Z width. 



2. R b 

The decay of the Z into bb is of particular interest, because it is sensitive through weak 
vertex corrections to the top quark mass, (see Fig. 9). The dependence of all the Z partial 
widths on the Higgs mass is due mostly to corrections to the Z propagator. Therefore, in 
Rb, the ratio of T b ^ to Y^ adrons , most of the Higgs and a s dependence cancels. This ratio 
then gives the only indirect measurement of the top mass independent of the Higgs mass. 
A comparison between the experimental measurement and the theoretical prediction is a 
particularly good test of the Standard Model and possibly the first place to look for new 
(non-Standard Model) physics. 

The value of Rb from the combinations of many measurements at LEP and SLC is higher 
than expected, whereas the value of R c is lower than expected. The measurement of Rb 
depends on what is assumed for R c because charm quarks are a background to the bottom 
quark signal in the data. Because the two are correlated, Rb is quoted either assuming a 
Standard Model value for R c , or the measured value. In Fig. [K] we show the LEP mea- 
surements of R c vs. Rb, together with the Standard Model theoretical prediction based on 
the direct measurement of the top quark mass performed at the Tevatron (see Section [TX| ). 
The disagreement between data and theory may be an indication of new physics beyond the 
Standard Model (see for example Altarelli et al. 1996; Wells and Kane, 1995). 



B. W mass 

The W vector boson mass (M w ) also depends on the top quark and Higgs masses through 
loop diagrams like those shown in Fig. 11, in which W —>■ tb — >W, or W — > WH — > W. A 
precise measurement of Mw constrains the top mass for a fixed Higgs mass. When combined 
with a precise measurement of the top mass, such a measurement can provide information on 
the Higgs mass, or, in the case of disagreement with theory, can signal the presence of new 
physics. Even with precise direct measurements of M top and Mw, however, the constraints on 
M Higgs are weak, because the Higgs mass dependence is only logarithmic. The present status 
of the W and top mass measurements, and their comparison with theory, is summarized in 
Fig. p. 
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FIG. 8. A summary of the measurements of sin 2 from LEP and SLC. The leftmost point is 



the value from A LR at SLD; sin 2 9 l £ JJ = 0.2305 ± 0.0005. The LEP average is sin 2 8%fi 
± 0.0004. (Woods, 1996). 
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FIG. 9. Diagrams showing t-quark corrections to Z — ► 66. 
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Preliminary 
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FIG. 10. Contours in the R& - R c plane derived from LEP data, corresponding to 68%, 95% and 
99.7% confidence levels assuming Gaussian systematic errors. The Standard Model prediction for 
M top = 180 ± 12 GeV/c 2 is also shown clS cL dot with an arrow through it. (This is the average of 
the 1995 CDF and DO top mass measurements; the current average is 175 ±8 GeV/c 2 ). The arrow 
points in the direction of increasing values of Mt op (The LEP Collaborations, 1995). 
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FIG. 11. Lowest order radiative corrections to the W mass involving top and bottom quarks and 
the Higgs. 
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CDF Measurements 

M„ = 80.41 ± 0.18 GeV/c 2 
M top = 176 ± 9 GeV/c z 

DO Measurements 

M w = 80.33 ± 0.27 GeV/c 2 
M lo = 170 ± 18 GeV/c 2 



CDF 



DO 



M T0P (GeV/c 2 ) 

FIG. 12. VF mass and top quark mass measurements from the Fermilab collider experiments (CDF 
and DO). The top mass values are from the full Tevatron data sets, with an integrated luminosity 
of ~ 100 pb _1 . The W mass values are derived from analyses of the first 15-20 pb~ x only. The 
lines are Standard Model predictions for four different Higgs masses (Flattum,1996). 



C. Global fits for M top and Mniggs 

Combining the indirect information from the neutral current experiments and the W 
mass measurement, a global fit for M top has been made by the LEP Electro- Weak Working 
Group (The LEP Collaborations, 1995). The fits are made with a s and M top as free parame- 
ters, since a s at the Z mass has a large uncertainty. The best predicted value for M top using 
data from LEP, SLC, the Fermilab collider W mass measurements and v N scattering data, 
is M top = 178 ± 8tH GeV/c 2 , with a s = 0.123 ± 0.004 ± 0.002 and with x 2 /d.o.f. = 28/14 
(where we have chosen M Higgs = 300 GeV/c 2 to quote the goodness of fit). The second un- 
certainty in this fit to the top mass comes from varying M#j 99S from 60 GeV/c 2 to 1 TeV/c 2 . 
The fit results are in good agreement with the directly measured values of a s and M top , 
a s (M z ) = 0.123 ± 0.006 (Bethke, 1995) and M top = 175 ± 8 GeV/c 2 , see Section 0. The 
variation of the fit-x 2 as a function of M top for three different choices of M#j ggs is displayed 



in Fig. 13 



In conclusion, all the neutral current data, as well as the W and top mass measurements 
are in agreement with each other, with the exception of the measurement of Rfe. The situation 
is nicely summarized in Fig. [Tl. In this figure, the correlation between R/ and R;, is due 
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to the fact that Ki depends on the total hadronic width, and hence on T(Z — > bb). Given 
the measured value of R/ and a s , and assuming Standard Model dependence of the partial 
widths on sin 2 Q l *fl for all but the 6-quarks, R/ constrains R& and sin 2 9 1 ^. These three 
measurements are compared with the Standard Model prediction given the measured top 
mass. The measured values of R& and R c are somewhat inconsistent with the combination 
of sin 2 9 l eff, R; and M top within the context of the Standard Model. 

Despite the large number of very precise measurements, there is still little information 
on the mass of the Higgs boson, although the data seems to prefer a low value for M Higgs , 
(see Fig. O). An estimate for the Higgs mass can be made using all the neutral current and 



hadron collider data as shown in Fig. |15| The best estimate of the Higgs mass is shown with 
and without the R& and R c measurements included. In all cases a Higgs mass less than 300 
GeV/c 2 is favored with large uncertainty. 
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FIG. 13. The x 2 curves for the Standard Model fit to the electro-weak precision measurements 
from LEP, SLD, CDF and DO (W mass only) and neutrino scattering experiments as a function 
of M^p for three different Higgs mass values spanning the interval 60 GeV/c 2 < M Higgs < 1000 
GeV/c 2 . The number of degrees of freedom is 14 (The LEP Collaborations, 1995). 
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FIG. 14. The combined LEP/SLD measurements of sin 2 O^j and R& assuming the Standard Model 
value of R c = 0.172 and the Standard Model prediction. Also shown is the constraint resulting 
from the measurement of R; on these variables, assuming a s {M%) = 0.123 ± 0.006, as well as the 
Standard Model dependence of light-quark partial widths on sin 2 O^ff (The LEP Collaborations, 
1995). 
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FIG. 15. Ax 2 = X 2 — Xmin vs ^-Higgs curves. The continuous line uses all the data from neutral 
currents and pp. The dashed line excludes the LEP and SLD and R c measurements. The dotted 
line excludes the SLD data (The LEP Collaborations, 1995). 
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IV. TOP QUARK PRODUCTION 



Because of its large mass, the top quark can only be observed directly in collider experi- 
ments, where sufficiently high center-of-mass energies (\/s) have been achieved. In electron- 
positron collisions, top quarks are produced in pairs via a photon or a Z. Since at lowest 
order this is a purely electro-weak process, the cross section and production kinematics can 
be precisely predicted. Today's highest energies e + e~ accelerators, LEP at CERN and SLC 
at SLAC, operate at y/s ~ M z m 91 GeV, and therefore the mass region M top > 46 GeV/c 2 
cannot be explored. Searches for top in e + e~ collisions will be briefly reviewed in Section [VI|. 
The HERA electron-proton collider at DESY has achieved a center-of-mass energy of ~ 310 
GeV. However the top production cross section at HERA is too small for observation and 
study of the top quark. 

Significantly higher center-of-mass energies have been achieved at hadron colliders. The 
pp collider at CERN (the SppS), which operated between 1981 and 1989, reached y/s = 
630 GeV; the pp collider at Fermilab (the Tevatron) came on line in 1987 with y/s = 1800 
GeV = 1.8 TeV; a new pp collider (LHC, y/s = 14 TeV) is under development at CERN, 
and is expected to begin operation in 2003. Until a very high energy e + e~ machine is built, 
top quark physics can be directly pursued only at hadron colliders. In this Section we will 
concentrate on the production of top quarks in pp collisions. Top quarks are produced by 
colliding partons (quarks, gluons, and antiquarks) from the proton and anti-proton. There- 
fore many aspects of our discussion also apply to other (e.g. pp) hadron-hadron collisions. 
Because the partons carry only a fraction of the momentum of the hadron, the center-of- 
mass energies of parton-parton collisions span a wide range of energies, see the discussion of 
parton luminosities in Section [IV. A . 

A. Production mechanisms and cross sections 

There are three mechanisms for top production in pp collisions : 

• Pair production of top quarks, pp — > tt + X. The leading order Feynman diagrams 
for this process are shown in Fig. 16. At higher order, gluon-quark scattering also 
contributes, tt pairs can also be produced through a Z or a photon, however the cross 
section is much smaller and we will not consider this possibility further. 

• Drell-Yan production of a W boson, with subsequent decay into tb, i.e. pp — >W + X; 
W — > tb, (see Fig. 17). Except for small contributions from off mass-shell W boson 
production, this mechanism only contributes for top masses smaller than M^— 

• Single top-quark production via H^-gluon fusion, see Fig 18. Photon-gluon and Z-gluon 
fusion are also allowed, with a much lower cross section. 

As will be illustrated in this Section, at Fermilab's Tevatron the strong-interaction pair 
production process (pp — > tt) is dominant over a wide range of top masses. For M top ~ 60 
GeV/c 2 , the top production rate from W — > tb is comparable to that of tt. For very high 
top mass, above M top rs 220 GeV/c 2 , the expected cross section for single top production 
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FIG. 16. Lowest order Feynman diagrams for production of tt pairs in pp collisions. 
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FIG. 17. Lowest order Feynman diagram for Drell-Yan production of tb pairs, pp — > W — > t6. 
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FIG. 18. Lowest order Feynman diagrams for production of a single top quark via VF-gluon fusion 
in pp collisions. When combining diagrams (a) and (c), care must be exercised to avoid double 
counting and to define the 6-quark distribution inside the proton in a consistent fashion (Heinson, 
Belyaev, and Boos, 1995; Carlson and Yuan, 1995). 



through VF-gluon fusion becomes larger than the pair production cross section due to the 
very high parton center-of-mass energy required to produce a ti pair (see discussion below). 

The pair production cross section for heavy quarks such as the top can be calculated in 
perturbative QCD. It factorizes as a product of the parton distribution functions inside the 
protons and the parton-parton point cross section, and is written as a sum over contributions 
from partons inside the proton and anti-proton (Collins, Soper, and Sterman, 1986) : 



a{pp -> ti) = J dxiFi(xi, /j 2 ) J dx j F j (x j , p 2 )aij(s, p 



2 M top ) 



The functions Fi and Fj are the number densities of light partons (quarks, antiquarks, 
and gluons) evaluated at a scale p in the proton and anti-proton; Xi and Xj are the momentum 
fractions of the incoming partons, i.e. parton i (j) has momentum X{P (—XjP), where P is 
the magnitude of the proton momentum in the center of mass frame (which in colliding beam 
experiments coincides with the lab frame); &ij is the point cross section for i + j —> ti and 
s = Ax the square of the center of mass energy of the parton-parton collision. 

The factorization and renormalization scale p is an arbitrary parameter with dimensions of 
energy, which is introduced in the renormalization procedure. The exact result for the cross 
section should be independent of the value of p. However, since calculations are performed 
to finite order in perturbative QCD, cross-section predictions are in general dependent on 
the choice of scale, which is usually taken to be of the order of M top . The sensitivity of 
perturbative calculations to reasonable variations in p is used to estimate the accuracy of 
the prediction. Parametrizations of the parton number densities (Fi and Fj) are extracted 
from fits to a large number of experimental results, mostly from deep inelastic scattering, 



see for example Fig. O. 



21 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

X X 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

X X 



FIG. 19. MT-B2 parametrization (Morfin and Tung, 1991) of xFi for valence u-quarks, xll v (x, /i 2 ), 
valence d-quarks, xD v (x, /i 2 ), gluons, xG(x,fi 2 ), and sea u- and d-quarks, xU{x,[i 2 ) + xD(x, /x 2 ). 
These are calculated at fi 2 = (20 GeV) 2 (solid line) and [? = (400 GeV) 2 (dashed line). The tt 
cross section is usually calculated with /i 2 = M 2 . 



The cross section for pp — > tt can also be written as (Eichten et al, 1984) : 

= ^-rr<7ij(s,ti ,M top ) 

where r = s/s and dLij/dr are the differential parton luminosities defined as : 

i 



_L_ J<^L [F J (x, / i 2 )F :; (r/x, / i 2 ) + F J (x, / i 2 )F J (r/x, A x 2 )]. 



dr ^ , wy 

The parton luminosities for quark-antiquark and gluon-gluon processes at Tevatron en- 
ergies are displayed in Fig. The sharp fall-off of these luminosities with increasing s, as 
well as the asymptotic 1/s dependence of a, result in predictions for the tt cross-sections 
which fall steeply as a function of M top . At <Js = 1.8 TeV the gg luminosity is larger than 
the qq luminosity up to s ~ (220 GeV) 2 . As a result, top pair production is dominated by 
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the gg — > tt process up to M top pa 90 GeV/c 2 . For higher top quark mass, qq initial states 
are the most important source of tt pairs. 



T 



C 



pp at Vs = 1 .8 TeV 



o 10 r 





=o 10 r 



-2 




t- -A 
-'10 r 



-5 : 
10 r 



-6 - 
10 r 

-7 : 



10 r 



-8 



10 r 



-9 : 
10 r 



-10 
10 



100 200 300 400 500 600 700 800 



V§~ (GeV) 



FIG. 20. Gluon-gluon (gg) and quark-antiquark (qq = uu + dd + ss) parton luminosities in pp 
collisions at y/s = 1.8 TeV. These are calculated using the MT-B2 parametrization of the parton 
distribution functions (Morfin and Tung, 1991) evaluated at a scale g 2 = s. 

The leading order (LO) cross section for producing a pair of heavy quarks in parton- 
parton collisions was calculated in the late 1970's (Gluck, Owens, and Reya, 1978; Corn- 
bridge, 1979; Babcock, Sivers, and Wolfram, 1978; Hagiwara and Yoshino, 1979; Jones and 
Wyld, 1978; Georgi et al, 1978). The full next-to-leading order (NLO) calculation was 
performed by Nason, Dawson, and Ellis in 1988, and shortly afterwards by Beenaker et al, 
1991. On the basis of their result, cross-section predictions were then made (Altarelli et 
al, 1988; Ellis, 1991) by convoluting the partonic cross-section with parametrizations of the 
parton distribution functions. 

There are two sources of uncertainty in such a calculation of the pp — > tt cross-section 
as a function of M top . As mentioned above, the first uncertainty is due to the nature 
of the perturbative QCD calculation for the partonic cross-section (a). The size of the 
uncertainty is customarily quantified by varying the arbitrary value of the scale /i by a 
factor of two around the top mass. Note that this is not a rigorous procedure, and it merely 
results in a reasonable estimate of the systematic uncertainty due to the missing higher 
order terms in the calculation. An additional uncertainty arises from the limited knowledge 
of the input parton distribution functions, and the assumed value of the QCD parameter 
h-QCD- The Aqot- dependence arises from the fact that the assumed value of Aqcd affects 
the /x 2 evolution of both a s and the quark and gluon distributions. In particular, the 
extraction of the gluon distribution from deep inelastic data also depends on Aq C d- The 
uncertainty on the cross-section calculation due to the parton distribution uncertainties are 
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Calculation 




Order 


a 


(tt) 




(1) Ellis, 1991 




NLO 


4.20^ 


-0.28 
-0.54 


pb 


(2) Laenen, Smith and van Neerven, 1994 


NLO + 


gluon resummation 


4.94;! 


-0.71 
-0.45 


pb 


(3) Berger and Contopanagos, 1995 


NLO + 


gluon resummation 


5.52^ 


-0.07 
-0.45 


pb 


(4) Catani et al, 1996 


NLO + 


gluon resummation 


4.75^ 


-0.63 
-0.68 


pb 



TABLE II. Calculations of the pp — ► ti cross-sections at Tevatron energies for M top = 175 GeV/c 2 . 
Note that these cross-section calculations use different sets of parton distribution functions (PDFs). 
The systematic uncertainties in (2) and (3) do not include the effects of varying the input PDFs. 



very hard to quantify. This uncertainty is usually estimated by studying the variations of the 
calculated cross-section using different parametrizations for the parton distribution functions 
and different values of Aqcd- As a result of these studies, the total theoretical uncertainty 
on the ti production cross-section at y^i =1.8 TeV is estimated to be of order ± 20%. The 
uncertainties due to the choice of scale and to the parton distribution assumptions are found 
to contribute approximately the same amount to the total uncertainty. 

The next-to-leading-order, O(af), predictions for the pair production cross-section have 
been subsequently refined by Laenen, Smith, and van Neerven (1992 and 1994). In their 
calculation the corrections due to initial state gluon brehmstrahlung, which are large near ti 
threshold, have been resummed to all orders in perturbative QCD, and have been included 
in the computation. This procedure introduces a new scale /j,q » Aqcd, where the resum- 
mation is terminated since the calculation diverges as /i — > 0, where nonperturbative effects 
are expected to dominate. Given that the corrections due to soft gluons have been shown 
to be positive at all orders in perturbative QCD for /i = M top , they estimate the lower limit 
on the tt cross-section as the sum of the full O(a^) prediction and the O(a^) soft gluon 
correction, using the conservative value of Aqcd — 105 MeV. Their best estimate of the 
cross section includes the full gluon resummation contributions, and the uncertainty arises 
mostly from the choice of /io, which is allowed to become as small as 0.05M top and 0.2M iop 
for the qq — > ti and gg — > ti channels respectively. 

A separate calculation of the ti cross-section including the perturbative resummation of 
gluon radiative corrections has become available in the past year (Berger and Contopanagos, 
1995). This calculation is based on Principal Value Resummation (PVR) techniques (Con- 
topanagos and Sterman, 1993 and 1994), and is independent of the arbitrary infrared cutoff 
fiQ. Theoretical uncertainties are estimated by varying the renormalization and factorization 
scale by a factor of two around the top mass. A more recent evaluation of the effects of 
gluon resummation suggests that its contribution is much smaller than previously thought 
(Catani et al, 1996). In Table [TI| we summarize the results of the various calculations of the 
pp — > ti cross-section at Tevatron energies for a top mass of 175 GeV/c 2 which, as we will 
discuss in Section |X|, corresponds to the directly measured value of the top mass. 

The expected top production cross-sections at Tevatron and SppS energies, for the three 
production mechanisms, are displayed in Fig. EH! and Fig. 1221. The Drell-Yan cross section, 



a(pp — > W — * tb), is calculated from the diagram shown in Fig. 17. The value of this cross 
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FIG. 21. Top production cross-sections in pp collisions at ^Js = 1.8 TeV. (a) pp — > ti from Laenen, 
Smith and van Neerveen, 1994. The band represents the estimated theoretical uncertainty; (b) 
sum of tb and tb from W decay (Drell-Yan); (c) sum of tb and ib from VF-gluon fusion. See text 
for details. 



section is normalized to the rate of W production through measurements of pp — > W — ■> ev 
(F. Abe et al, 1991b; Alitti, et ai, 1990a; Albajar et ai, 1991a), including corrections for 
the phase-space suppression of a tb pair and the finite W width. The tree-level W^-gluon 
fusion cross-section has been calculated by several authors (Willenbrock and Dicus, 1986; 
Yuan, 1990; Anselmo, van Eijk, and Bordes, 1992; Ellis and Parke, 1992; Bordes and van 
Eijk, 1993); the cross sections shown in Figs. |2l] and |22] are obtained using the PYTHIA 
Monte Carlo event generator (Sjostrand and Bengtsson, 1987). Since the I^-gluon fusion 
matrix element is calculated at tree level, the systematic uncertainties on the absolute rate 
prediction can be large, see for example Ellis and Parke, 1992. Here we have used the default 
PYTHIA scale /i 2 = 0.5 (M 2 X + M 2 2 ), where Ma and M t 2 are the transverse masses of the 
outgoing partons. Recently, a calculation of the next-to-leading order QCD corrections 
for the PF-gluon fusion process has been performed (Bordes and van Eijk, 1995), and the 
enhancement of the cross section over the Born-level result has been found to be of order 
30%. It is worth mentioning that, despite their apparent similarities, the Drell Yan and 
H^-gluon mechanisms are quite distinct. The higher order Drell-Yan diagram, qg — > qW*, 
W* — > tb (see Fig. |23) and the H^-gluon fusion diagrams (see Fig. 18) have the same initial 
and final state partons. However, in one case the W is space-like, and in the other case it is 
time-like. Furthermore, the tb pairs in the two processes are in different color states. 

As anticipated, the pair production cross-section at the Tevatron (see Fig. [21]) is dominant 
up to very high mass, except for the top mass region around 60 GeV/c 2 , where the cross 
sections for the ti and W —>■ tb processes are approximately equal. On the other hand, at 
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FIG. 22. Top production cross-sections in pp collisions at y/s = 0.63 TeV. (a) pp — > ti from Laenen, 
Smith and van Neerveen, 1994. The band represents the estimated theoretical uncertainty; (b) sum 
of tb and ib from W decay (Drell-Yan); (c) sum of tb and ib from PF-gluon fusion. See text for 
details. 



the lower energy of the SppS collider, top production through W decay dominates in the 
mass region 40-80 GeV/c 2 , (see Fig. F2I|). Once the experimental evidence started pointing 
towards higher top masses (see Section |VI|) , it became clear that top searches at the SppS 
were not competitive with those at the Tevatron, due to the lower it cross-section at \fs = 
630 GeV. 

It should be emphasized that top production is a very rare process in pp collisions. The 
total inelastic cross-section at the Tevatron is approximately 60 mb (F. Abe et al, 1994b), 
ten orders of magnitude higher than a(ti) for M top = 175 GeV/c 2 . Therefore, in trying to 
isolate a top signal, both excellent background rejection and high luminosities are critical. 

In the remainder of this Section and in most of this review we will concentrate on the 
pp — > ti reaction. We will however revisit the W^-gluon fusion process in Section [X|, since 
it is interesting in its own right and its study will become accessible in the not too distant 
future. 



B. Top quark hadronization 

Quarks are not observed as free particles but are confined to form hadronic bound states. 
The top quark however is unique in that its mass is high enough that it can decay before 
hadronization. According to the Standard Model, top quarks undergo the weak decay t — > 
Wb, where the W boson is real if M top > M w + M b , and virtual otherwise. Decay modes such 
as t — > Ws and t — > Wd are also allowed. They are suppressed by factors of | V^ s | 2 /| V t b\ 2 ~ 
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FIG. 23. Diagram contributing to the 0(a s ) corrections to the Drell-Yan process pp —> W —* tb. 

1CT 3 and |Vy 2 /l V t b\ 2 ~ 5 x 10~ 4 respectively, where Vy is a Cabibbo-Kobayashi-Maskawa 
(CKM) mixing- matrix element (Montanet et al, 1994). The expected width of the top 
quark, and hence the lifetime, as a function of its mass, is shown in Fig. [24]. 
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FIG. 24. The Standard Model width of the top quark as a function of its mass (Bigi, 1986). Note 
the transition between the region of virtual and real W decays which occurs at Mt op ~ Myy + 

The hadronization process, which is non-perturbative in nature, is not well understood. 
However, the formation of hadrons is estimated to take place in a time of order Aqq D rs 
O(100 MeV)- 1 « O(10~ 23 ) seconds (Bigi, 1986). As can be seen in Fig. |4[ the top lifetime 
becomes shorter than this characteristic time if the top mass is higher than approximately 
100 GeV/ c 2 . A more quantitative treatment is given by Orr (1991), and is briefly summarized 
here. In this model, the t and the t emerging from the hard scatter are linked by color strings 
to the remnants of the proton and anti-proton. When the separation between the outgoing 
quarks and the color-connected remnants exceeds a distance of order 1 fm, the stretched 
color string is expected to break, resulting in the creation of fragmentation particles from 
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the vacuum, and possibly the formation of a bound state top hadron. For M top > 165 
GeV/c 2 , the top production kinematics at the Tevatron are such that essentially all top 
quarks are expected to decay before having travelled that minimum distance. Conversely, 
for top masses below 120 GeV/c 2 , the overwhelming majority of top quarks will survive to 
a distance of 1 fm, and hadronization effects are expected to occur. Varying the assumption 
on the hadronization distance by a factor of two results in a mass-shift of order 20 GeV/c 2 
for the transition region between hadronization and free-quark decay. 

The exceedingly short top-quark lifetime is due not only to the very high mass, but also 
to the fact that M top > M b , so that the the top can decay into Wb, and this decay mode is not 
CKM suppressed. A heavy I 3 = —1/2 fourth generation quark (V) would decay into Wu', 
where vl here stands for an I3 = +1/2 up-type quark. The decay rate would be proportional 
to the square of the CKM mixing-matrix element that connects b' and vl . If the b' were lighter 
than the fourth-generation up-type quark, then only generation-changing, CKM suppressed, 
decays would be allowed. As a result, the lifetime of such a fourth generation quark could 
be considerably longer than that of a top-quark of the same mass. Fourth generation heavy 
hadrons, as well as new quarkonia states, could then still be allowed to form. 

Even if hadronization effects do occur in tt production, their effects, although potentially 
interesting, are not expected to be experimentally observable, at least in the forseeable 
future. The reason is that the fragmentation of heavy quarks is hard, i.e. the fractional 
energy loss of the top quark as it hadronizes is small (Peterson et al, 1983). Distortions 
to the kinematics of the top quark from the perturbative partonic calculation are minimal. 
Additional particles produced in the hadronization process have little effect on the overall 
event topology. If a top hadron is indeed produced, the kinematics of the top decay will 
not be very different from that of a free quark decay since the companion quark is so much 
lighter. The fragmentation of the 6-quark produced in top decay could potentially be more 
seriously affected. This is because the color string would link the 6-quark to the light quark 
produced in the top fragmentation rather than the proton or anti-proton remnant (Orr, 
1991). However, all top-quark experimental studies to date have not been precise enough to 
be sensitive to fragmentation assumptions. 

C. Underlying event 

After the hard collision, the remnants of the proton and anti-proton also hadronize. This 
process cannot be described within the framework of perturbative QCD, and is therefore 
poorly understood. The particles from the remnant hadronization form what is usually 
referred to as the underlying event. The structure of the underlying event is similar to that 
of the bulk of soft pp collisions (the so-called minimum bias events). 

Minimum bias events are events collected with a simple interaction trigger. This trigger 
usually consists of a coincidence between large banks of scintillator counters in the very 
forward and backward regions, and is highly efficient for all types of inelastic pp collisions, 
except for singly diffractive events. In minimum bias events, the average transverse mo- 
mentum of the hadrons is 500 MeV/c (Para, 1988). (Transverse momentum is the 
component of momentum perpendicular to the direction of the beams) . Most of the energy 
is carried away by particles which remain inside the beam pipe and are not seen in the 
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FIG. 25. The expected rapidity (y) and transverse momentum (Pt) distributions for top quarks 
at the Tevatron. The predictions are form the ISAJET Monte Carlo program. 

detector. The charged particle multiplicity in the central region per unit pseudo-rapidity 
(dN ch /di]) grows approximately logarithmically with the center-of-mass energy, and at = 
1.8 TeV is dN ch /dr] ~ 4 (F. Abe et al, 1990e). (NB: the rapidity y of a particle is defined in 
terms of the longitudinal Lorentz boost, with (3 = tanlny, to the frame in which the particle's 
momentum is purely transverse. Rapidity can be written as 
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1 E + P z 
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2 E-P 7 
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where E is the energy of the particle, M is its mass, and P z and Pp are the components of 
momenta parallel and transverse to the beam direction. Pseudo-rapidity {rf) is the rapidity 
calculated neglecting the particle's mass, rj = — In tan(6 l /2), where 9 is the polar angle with 
respect to the proton direction. Zero rapidity or pseudorapidity corresponds to particles 
moving at 90° from the beamline; high values of \y\ or \t]\ implies very forward or backward 
going particles). 



D. Modelling of top quark production 

The reliability of the modeling of tt production in is an important issue. Top production 
is usually modelled using a QCD shower Monte Carlo program, such as ISAJET (Paige 
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and Protopopescu, 1986), HERWIG (Marchesini and Webber, 1984 and 1988) or PYTHIA 
(Sjostrand and Bengtsson, 1987). These Monte Carlos are used by the experimental groups 
to calculate the tt acceptance and kinematics, and to model the resolution of the top quark 



mass measurement, see Section IX. 



In all these Monte Carlo programs, the initial hard scatter is generated from tree-level 
matrix elements convoluted with parametrizations of the parton distribution functions. Ini- 
tial and final state partons are then developed into a gluon and qq radiation cascade, with 
angular and energy spectra based on the QCD Altarelli-Parisi splitting function. The QCD 
shower is terminated when the virtual invariant mass of the parton in the cascade becomes 
smaller than a minimum value which is of order 1 GeV for HERWIG and PYTHIA and 
6 GeV for ISAJET, at which point perturbative QCD is expected to break down. Phe- 
nomenological models are then employed to combine the remaining partons into hadrons. 
The underlying event is also modelled in a phenomenological way, with a number of pa- 
rameters tuned to reproduce the hadron multiplicities and transverse momentum spectra 
measured in soft pp collisions (minimum bias events). Short lived particles are made to 
decay, with branching ratios and decay models based on the compilation from the Particle 
Data Group (Montanet et al, 1994). 

The main differences between these Monte Carlo event generators reside in the model- 
ing of the radiation processes. ISAJET employs an independent fragmentation model, i.e. 
radiation from each parton occurs independently from the structure of the rest of the event, 
whereas in both HERWIG and PYTHIA radiation is more realistically emitted taking into 
account color-correlations between all partons in the initial and final state. The output of 
these Monte Carlo event generators consists of a list of stable particles, which can then be 
fed to a detector simulation for detailed studies of the expected signature of a top event. 

In the pair production process, the t and t quarks are produced in the central rapidity 
region, with Pt of order M top /2, (see Fig. ^5|). These features of tt production can be simply 



understood from the properties of the i+j—*tt process. The cross-section for the qq — > tt 
subprocess, which dominates at high top mass, is given at lowest order by (see for example 
Nason, Dawson, and Ellis, 1988) : 
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The parton-parton cross-section as a function of s rises from zero at threshold (s = 4M^ op ), 
reaches a maximum at s — 5.6M^ op , and then falls off asymptotically as 1/s. When convo- 
luted with the falling qq luminosity, (see Fig. ^0[), the maximum of the qq — > tt cross-section 
is shifted down to s ~ A.5Mf op . Therefore, the most probable energy for a top quark is E 
~ V4.5M top /2 « l.lM top , and the most probable momentum is P ~ 0.4M top . 

From the definition of rapidity, it is clear that the maximum of \y\ occurs as Pt — > and 
at maximum E, which for pair-produced objects is E = 0.5i/i = 900 GeV at the Tevatron. 
For M top = 160 GeV/c 2 , the kinematic limit is then \y\ < 2.4; however, as can be seen from 



Fig. |25], most top quarks have \y\ < 1.5. High values of \y\ are suppressed because they 
require Pt — > 0, where the phase space factor also — > 0, and they require high values of E, 
i.e. high values of s. Both the parton-parton luminosities and parton-parton cross-section 
fall off with increasing s. 
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The lowest order diagrams (Fig. 16) lead to a back-to-back topology for the t and the t 
in the transverse plane, which is slightly modified by higher order corrections. Because the 
top quark momentum is not large compared to M top , the decay products are not significantly 
boosted along the original top quark flight path, leading to nearly spherical events. In the 
next Section we will turn to the discussion of top decays and signatures. 
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FIG. 26. Comparison of the transverse momentum distributions of the tt pair as predicted from 
the Herwig Monte Carlo and next-to-leading order (NLO) QCD. At Tevatron energies, and for 
M top = 176 GeV/c 2 . From Frixione et al, 1995. 

As was mentioned above, shower Monte Carlos are based on the LO matrix element for 
tt production and models of initial and final state radiation. It is interesting to compare 
these models with higher order QCD calculations. The earlier NLO order calculations of 
tt production (Nason, Dawson, and Ellis, 1988; Beenaker et al, 1991) are not sufficient, 
since these are calculations of single quark kinematic distributions such as Pt and rapidity, 
integrated over the whole phase-space for the other quark. More recently, a NLO calcula- 
tion of the doubly-inclusive cross-section for heavy quark production has become available 
(Mangano, Nason and Ridolfi, 1992). This calculation allows for comparisons not only of 
single t or i distributions, but also of correlated distributions, e.g. the tt invariant mass 
(M(tt)), the tt transverse momentum (Py(tt)), and the azimuthal separation between t and 
t (A0). Detailed comparisons between the Herwig model and NLO-QCD have been per- 
formed (Frixione et al, 1995). Excellent agreement is found in the shapes of distributions 
of quark rapidity, Pt, and M(tt), except for very large values of the latter two quantities. 
In this kinematic regime, multiple gluon emission from the final state top quarks becomes 
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important, and this process is not modelled by the NLO QCD calculation. Disagreement 
between NLO QCD and HERWIG is also observed in distributions of Pr(ti) and A0, (see 
Fig. |2(|). Note that at LO these distributions are delta-functions, with Fritt) = and A(f> 
= 180°, and deviations from the delta-function behavior are due entirely to higher order 
corrections. For small Pr(tt), multiple gluon emission is expected to dominate, and the 
HERWIG model is expected to be more realistic. 

Studies of the expected single top quark and rapidity distributions, calculated in 
NLO QCD and including the resummation of the leading soft gluon corrections, have also 
been performed (Kidonakis and Smith, 1995). The shapes of these distributions are found 
to be essentially identical to those calculated at NLO, which were shown to agree with the 
Herwig model. Comparisons of gluon emission in ti events from HERWIG and from a O(a^) 
matrix element calculation, including initial and final state gluon radiation have also been 
made (Orr, Stelzer, and Stirling, 1995). These authors find larger contributions of gluon 
radiation in HERWIG than in the matrix element calculation. This effect may also be due 
to the absence of multi-gluon emission in the calculation. 

As we will discuss in Section 1X1 understanding gluon radiation in ti events is crucial for 



a precise determination of the top mass. We expect that this subject will attract more and 
more attention in the next few years. 



V. TOP QUARK SIGNATURES 

Since the top quark decays with a very short lifetime, only its decay products can be de- 
tected. Therefore, to understand the experimental signature for a top event, we first discuss 
the decay modes of the top quark. In this Section we review the top quark decay properties 
and we present a general discussion of how the top quark can be observed, paying particular 
attention to the background sources. We will concentrate on the signature for pp — > ti, since 
this is the most important production mechanism at Tevatron Collider energies. 



A. Standard Model top quark decay modes 

As mentioned in the previous Section, according to the Standard Model the top quark 
decays as t — > Wb, where the W boson is real or virtual depending on the top mass. (Non- 
standard Model decay modes of the top quark will be reviewed in Section |V1.C| ). The W will 
subsequently decay into fermion pairs, either W —>■ Iv or W — > qq, where I denotes a charged 
lepton, and qq denotes a light-quark pair, ud or cs. At tree level, the W couples with equal 
strength to leptons and quarks, so each W decay mode occurs with equal probability. There 
are three leptonic channels (ev, \w, and rv), and six hadronic channels (ud and cs, with 
three possible color assignments), hence each decay mode has a branching ratio of 1/9. QCD 
corrections enhance the branching ratios of the hadronic modes by a factor of (1 + a s /ir) 
~ 1.05. Given the W branching fractions, it is a simple matter to list the ti decay modes, 



see Table 111. 
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Decay mode Branching ratio 

ti -> qq qq bb 36/81 

tt -> qq eu bb 12/81 

ti^qqiivbb 12/81 

tt^qqTvbb 12/81 

tt^ev^ivbb 2/81 

tt -> ez/ rzv 66 2/81 

tt^iivTubb 2/81 

ti^evevbb 1/81 

— > /iz/ /iz/ 66 1/81 

tt^TVTvbb 1/81 



TABLE III. Decay modes for a pair and their lowest order branching ratios assuming Standard 
Model decays. 



B. Detection of the top decay products 

The possible final states contain combinations of electrons, muons, taus, neutrinos, and 
quarks. Here we briefly illustrate techniques for detection of the top-quark decay products. 

General-purpose pp collider detectors are needed for top quark physics. These detectors 
are designed to cover as much as possible of the solid angle around the interaction point, 
and are composed of a number of sub-detectors optimized for study of different aspects of 
the event. A number of such detectors have been used (UAl and UA2 at the CERN SppS), 
are still in operation (CDF and DO at Fermilab's Tevatron), or are now being designed and 
constructed (CMS and ATLAS, at the proposed LHC pp collider). While the details of 
the design of these detectors are different, their overall structure is in general quite similar. 
The region immediately surrounding the interaction region is instrumented with detectors 
designed to measure the trajectories of charged particles. Except for UA2 and DO, the 
tracking volume is immersed in a magnetic field for momentum measurement. The tracking 
volume is surrounded by calorimeters, where measurements of the energy of electromagnetic 
and hadronic showers are performed. Calorimeters are segmented both longitudinally and 
transversely to the direction of flight of particles originating from the interaction point. 
Transverse segmentation is necessary to measure the position of the showers, while lon- 
gitudinal information is used to distinguish between electromagnetic (EM) and hadronic 
(HAD) showers. Calorimeters cover most of the Att solid angle around the interaction re- 
gion, however the very forward and backward regions must be left uninstrumented to allow 
for the passage of the beam pipe. Muon detectors consisting of additional tracking devices, 
hadron absorbers, and possibly magnets for momentum measurement, are placed outside 
the calorimeter. Drawings of the two collider detectors at Fermilab's Tevatron are shown in 



Fig. 27 and Fig. 28. 
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FIG. 27. The DO detector at the Tevatron. EC and CC are liquid-argon-uranium calorimeters. 
The central detector provides tracking information. Muons are detected using the five toroids (CF, 
EF, SAMUS) and the proportional drift tube (PDT) systems, CM and EM. From Snyder, 1995a. 
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FIG. 28. A side-view cross-section of one quadrant of the CDF detector at the Tevatron. The 
detector is forward-backward symmetric about the interaction region, which is at the lower-right 
corner of the figure. SVX, VTX, CTC, and CDT are tracking detectors. CEM, CHA, WHA, PEM, 
PHA, FEM, and FHA are calorimeters. CMU, CMP, CMX, and FMU are muon detectors. BBC 
is a bank of scintillators. CPR and CES are multiwire proportional chambers placed in front and 
in the middle of the central electromagnetic calorimeter (CEM). From F. Abe et al, 1994a. 



1. Detection of electrons and muons 

Electrons are identified as highly electromagnetic showers in the calorimeters. If momen- 
tum information from the tracking system is available, consistency between the measured 
momentum of the electron candidate and the energy of the corresponding EM-shower pro- 
vides a powerful handle for rejection of backgrounds from e.g. hadronic shower fluctuations 
and overlaps between hadron tracks and photons from tt° decay. Information from the trans- 
verse and longitudinal shapes of the shower, from ionization measurements in the tracking 
chamber (dE/dX), and from the response of transition radiation and preshower detectors, 
are also used for electron identification. 

Electrons from W decays in top events have high transverse momentum, (P^, see Fig. 
and are expected to be isolated, i.e. well separated from the other decay products of the two 
top quarks in a ti event. These electrons can be identified with high efficiency, and their 
energy can be measured very precisely in the calorimeter (see Table fV\j . On the other hand, 
identification of electrons from b —>■ cev in a top event is much more problematic. These 
electrons have lower transverse momentum than electrons from W decays and, since the b is 
highly boosted, the nearby hadrons from the b fragmentation and b or c decay may deposit 
their energy in the same calorimeter cells as these electrons. 

Muons can also be reliably identified as charged particles that penetrate the calorimeter 
and reach the outside muon detectors. Backgrounds to the muon signal arise from decays 
in flight of pions and kaons, and from hadrons that traverse the calorimeter and hadron 
absorbers without interacting (punchthrough) . If there is a magnetic field in the inner 
tracking system (before the calorimeter), then the muon momentum is precisely measured; 
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Detector 


Electron energy resolution 


Muon momentum resolution 


UAl 
UA2 
CDF 
DO 


<t(E)/E = (0.15-0.21)/^ .03 
<r(E)/E = 0.17/VE © .02 
cr(E T )/E r = 0.14/VW © -02 

<t(e)/e = o.i5/Ve © .01 


<r(P)/P = 0.005P 
(t(P t )/P t = 0.0009P r 0.0066 

<t(p)/p = o.oip e .2 





TABLE IV. Electron energy (GeV) and muon momentum (GeV/c) resolutions in the central region 
for the UAl (Albajar et al, 1989), UA2 (Alitti et al, 1992a), CDF (F. Abe et al., 1994a), and DO 
(Abachi et al, 1994) detectors. The UAl muon momentum resolution is for measurements in the 
central-detector, for muons at 90° from the direction of the dipole field. The UAl electron energy 
resolution changed between 1983 and 1985 due to radiation damage of the scintillator. The symbol 
© indicates that the two terms are added in quadrature. The subscript T refers to components 
transverse to the beam direction. 



otherwise, (e.g. in the DO experiment) it is measured, with worse resolution, in the outer 
muon detector (see Table [TV]) . 



2. Detection of quarks 

Quarks hadronize and are detected as collimated jets of particles, (see Fig. ^). Jets in 
pp collisions are reconstructed by summing up the energy deposited in the calorimeter cells 
within a fixed cone in 77 — space, where 77 is the pseudorapidity, and <fi is the azimuthal angle 
around the beamline. The fixed cone algorithm is used because jets are are approximately 
circular in 77 — space, because the 77 — size of a jet of a given P T is independent of the 
rapidity of a jet, and because this size is only weakly dependent on the transverse momentum 
of the jet, as we will briefly discuss below. 

If the typical longitudinal and transverse momentum components of a fragmentation 
particle with respect to the jet axis are qr and qi, then the typical spread of the jet will be 
A9 « qr/qi and A0 ps qr/ili sin9), for q T « q^ Then, A77 = (dr}/d6) A9 ps — <Zt/(<7/ sin#) 
= —A<f), i.e. jets are approximately circular in 77 — <fi space. Since the rapidity of a massless 
particle under a longitudinal boost changes as y — > y + Ay, where Ay only depends on the 
boost, in the limit that the mass of the fragmentation hadrons is small, the 77 size of a jet of 
a given is invariant under longitudinal boosts, i.e. independent of the 77 of the jet itself. 
The size of a jet does vary slightly with its transverse momentum; for example a simple 
model of jet fragmentation uniform in rapidity along the jet axis predicts that the angular 
size of the cone containing half of the particles in the jet varies as 1/yE, where E is the 
energy of the jet. 

The size of the cone used in jet reconstruction must be matched to the size of a jet. 
On average, of order 70% of the jet energy is contained within a cone of radius AR = 



VA77 2 + A(f)' 2 = 0.4 (F. Abe et al, 1991c and 1993a; Linneman, 1995.). See also Fig. 131 



NB: AR< 0.4 translates into A0 < 23° in Fig. gl. 

The energy of a jet is defined as the energy of the corresponding calorimeter cluster. The 
resolution in the measurement is typically only of order ct(Et)/Et ~ 1.0/ \[Et (Et in GeV), 
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FIG. 29. The expected lepton transverse momentum from t — ► — > /i^ from the ISAJET Monte 
Carlo generator. This is for pp — > ti at -^s = 1.8 TeV. Detector effects are not included. 



(see Fig. |^). This poor resolution is due to (i) the intrinsic large fluctuations in the response 
of calorimeters to hadronic showers, (ii) differences in the calorimeter response between 
charged hadrons and electrons or photons, (iii) energy loss in uninstrumented calorimeter 
regions, e.g. in the vicinity of boundaries between calorimeter modules, (iv) energy loss 
due to the use of a finite cone-size in jet reconstruction, and (v) overlaps between the jet 
and hadrons from the underlying event. The direction of a jet is measured by linking the 
position of the energy cluster in the calorimeter with the position of the interaction point. 
The resolution on the angular measurement is a few degrees. 

The number of detected jets for a given decay mode in a tt event is not expected to 
correspond to the number of quarks in the final states listed in Table III. There are a 
number of reasons for this. First, as the of the parton becomes small, identification 
of the corresponding jet becomes more and more problematic, as it tends to blend with 
the underlying event. In practice, one imposes a minimum cutoff on the jet transverse 
momentum which is set to a value at least of order 10 GeV/c. Furthermore, as will be 
discussed further in this Section, backgrounds to the top signal consist mainly of events 
with low jets. Therefore, to achieve the needed background rejection, the minimum 
jet Pt requirement is often chosen to be higher than 10 GeV/c. A second source of jet 
reconstruction inefficiency is jet merging. Nearby jets can be resolved only if their separation 
in rj — <p space is larger than a minimum distance of the order of the clustering radius used 
in jet reconstruction. Top events have a large number of partons in the final state, and 
the probability that at least two of them will be too close to be separately identified is 
substantial. In those cases, the two nearby jets are merged and are reconstructed as a single 
jet. 
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FIG. 30. A pp — > jet-jet event in CDF. Here we show the reconstructed tracks in the transverse 
plane. The two-jet structure is apparent. 



To illustrate some of these effects, we show in Fig. |33| the expected transverse momenta 
and separation (AR) in rj — space for quarks in tt — > qq Iv bb events for M top = 170 
GeV/c 2 . Different experiments and different analyses use different cone clustering radii, 
typically between 0.3 and 1.0. In a significant fraction of events, the minimum AR between 
quarks is small enough that at least two of the quark-jets are expected to be merged. Even 
for a high top mass, the fraction of events with at least one relatively soft jet is substantial 
(again, see Fig. |33|). For lower top masses of course the transverse momenta will be even 
lower. The situation for M top close to Mw is particularly difficult. In that case the kinetic 
energy liberated in the t — > Wb decay is low and the b momentum in the top rest frame 
is small. Even after boosting to the laboratory frame, the b momentum remains soft, (see 
Fig. R§. 
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FIG. 31. The relative distribution in calorimeter cells at an angle <p with respect to the trans- 
verse jet thrust axis for CDF di-jet data with a 30 GeV Ey jet trigger threshold. Note that in 
these di-jet events one jet is at <j> = 0°, and the other jet is at = 180°. From F. Abe et. al, 1991c. 
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FIG. 32. From the DO experiment, Abachi 1995d; jet energy resolution as a function of jet trans- 
verse energy (Ey) as computed from di-jet and photon-jet events in four pseudo-rapidity regions. 
The fits are of the form (a E /E) 2 = (N/E) 2 + (S/^E) 2 + C 2 . The fitted values of TV, S, and C 
are : (a) N = 7.07, S = 0.81, C = 0.0; (b) N = 6.92, S = 0.91, C = 0.0; (c) N = 0.0, S = 1.45, 
C = 0.052; (d) N = 8.15, S = 0.48, C = 0.0. Jets are reconstructed using a cone-size of 0.5. 
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FIG. 33. (a) The expected transverse momentum of the lowest (solid) and highest (dashed) 
quark in lepton + jets events; (b) AR between the closest two quarks in pp — > ti — > qq Iv bb. 
Results are from the ISAJET Monte Carlo event generator at y/s = 1.8 TeV. Detector effects are 
not included. 
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FIG. 34. The expected transverse momentum distribution of 6-quarks from top decay for Mt op = 
170 GeV/c 2 (solid) and M top = 90 GeV/c 2 (dashed). Results are from the ISAJET Monte Carlo 
event generator for the process pp — > tt at ^/s = 1.8 TeV. Detector effects are not included. 
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The situation gets even more complicated when initial and final state radiation are taken 
into account. Initial state gluon radiations gives a net to the ti system, and therefore 
alters the jet spectrum calculated at tree-level; the radiated gluons can be detected as 
additional jets in the final state; large angle radiation from the final state quarks softens the 
spectrum of reconstructed jets, and can also result in additional jets. 

All of these jet reconstruction effects are important, and must be studied using QCD- 
inspired Monte Carlo event generators, in conjunction with a simulation of the detector 
response, see for example Fig. They constitute one of the major systematic uncertainties 
in the calculation of the ti acceptance and, more importantly, in the determination of the 
top mass, see Section fX. 
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FIG. 35. Jet multiplicity in tt — ► lubqqb (lepton + jets) events from the ISAJET Monte Carlo and 
the CDF detector simulation. Solid line : M top = 200 GeV/c 2 ; dashed line : M top = 120 GeV/c 2 . 
Jets are reconstructed using a cone-size AR=0.4, and must have \r/\ < 2. The jet transverse energy 
threshold is 15 GeV, without application of jet energy corrections. (The jet energy corrections will 
be described in Section IX. A. 3; a 15 GeV jet in CDF is corrected on average to 23 GeV). Note 
that in the absence of gluon radiation these events should have Nj e t s < 4; the significant fraction of 
events with additional jets in the final state is an indication of the importance of gluon radiation. 



3. Detection of neutrinos 

Neutrinos are detected by missing momentum techniques: since the initial center-of- 
mass momentum is zero, the vector sum of the momenta of all of the neutrinos in the event 
is inferred as the negative of the vector sum of the momenta of all the detected particles. 
However, because the very forward and backward detector regions are uninstrumented, longi- 
tudinal information is lost, and only the transverse components of the momenta of neutrinos 
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Detector 




Reference 


UAl 
UA2 
CDF 
DO 


0.7 ^E T (in GeV) 
0.8 (Y,E T ) 0A (in GeV) 
0.7 y/EE T (in GeV) 
1.08 GeV + 0.019 ZE T 


Albajar et al., 1989 
Alitti et al, 1990b 
F. Abe et al, 1994a 
Abachi et al, 1995d 





TABLE V. Missing transverse energy resolution for minimum bias events. 

can be measured. In practice, what is measured is not the momentum of all of the particles, 
but rather the energy deposited in the calorimeter. The missing transverse energy vector 

(JZl T ) is defined as % T = —T^E^, where the sum is over all calorimeter cellss, and E % T is a 
vector whose direction points to the ith cell and whose magnitude is equal to the transverse 
energy deposited in the ith tower. The missing transverse energy vector must be corrected 
for detected muons which only loose a minimal amount of energy in the calorimeter, and is 
then associated with the neutrino transverse momentum. 

The resolution on the neutrino energy measurement is very much dependent on the event 
topology, since it depends directly on the resolution in the measurements of all the leptons 
and jets in the event. Since leptons are in general well measured, the uncertainty in the 
measurement of JZl T arises mostly from errors in the measurements of jet energies. It is 
customary to parametrize the resolution in JZl T , c(E't), as a function of the total transverse 
energy in the event, T>E T , see Table [V]. 

4. Detection of tau leptons 

Taus are very hard to identify. Approximately 36% of the time, a tau lepton will decay 
into a muon or an electron. The signature for an event with a t — > W — > tv — > Ivv, where 
I = e or /i, is very similar to that of an event with a t — > W — > Iv decay, except that 
the final state lepton will in general have lower momentum. Taus that decay hadronically 
are detected as jets. Separation between jets from hadronic decays of taus and quarks or 
gluon jets in pp collisions has been achieved in e.g. measurements of the pp —> W — > tv 
cross-section (Albajar et al, 1989; Alitti et al, 1991b; F. Abe et al, 1992b), and in searches 
for non Standard- Model top quark decays, see Section |VT.C[ The separation is based on the 
distinctive narrowness of a jet from the hadronic decay of a high tau-lepton, and/or the 
characteristic one- and three-prong track multiplicities. However, the efficiency for detecting 
hadronic taus is so low, and the backgrounds from jet fluctuations are so high, that these 
techniques are only now just beginning to be applied successfully in the context of a Standard 
Model top search. 

For the remainder of this article, we will refer to ti final states with zero, one, or two 
leptons (e or /i) from W decay as all hadronic, lepton + jets, and dilepton respectively. Sig- 
natures that include the explicit identification of hadronic tau-decays will not be considered 
in this review. We now turn to a discussion of the ti signatures in these three channels. 
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C. All hadronic mode 

The all hadronic final state, (ti — > qq qq bb, see Table |111| ) is the most common, but it 
competes with very high backgrounds from pp — > 6 jets (see for example Benlloch, Wainer, 
and Giele, 1993). The cross-sections for this QCD process at the Tevatron is higher than 
the top cross-section by approximately three orders of magnitude. Despite the extremely 
high background levels, it may be possible, with sufficient statistics, to isolate a top signal 
in this mode by applying further kinematic cuts, and by identifying the fr-quark(s) in the 
final state. These issues are being carefully studied by the experimenters (Castro, 1994; 
Tartarelli, 1996; Narain 1996). In this review we will concentrate on the dilepton and lepton 
+ jets modes. 



D. Dilepton mode 



The signature for the dilepton final state (ti — ► Iv Iv bb, see Table p|) consists of two 
leptons, two 6-jets, and % T from the two neutrinos. Since the leptons originate from W decay, 
they tend to be isolated and to have high transverse momenta. The JZl T is also expected to 
be high, (see Fig. ^). Typical minimum lepton transverse momentum or JZl T requirements 
are set around 20 GeV/c. Despite the low branching ratio, this mode turns out to be very 
important because background levels are very low. 
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FIG. 36. Expected sum of neutrinos' transverse momenta in the dilepton channel (pp — > tt — > 
Ivblvb) at y/s = 1.8 TeV. From the ISAJET Monte Carlo event generator, for M top = 160 GeV/c 2 . 



The most probable way to obtain two isolated leptons in pp collisions is through the 
Drell-Yan process, see Fig. 37. This mechanism yields e + e~ and pairs, but not e/i 
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FIG. 37. Lowest order Feynman diagram for Drell-Yan production of lepton pairs. 



pairs, except through rr production, with both taus decaying leptonically. Here we begin 
by addressing backgrounds from direct Drell-Yan production of ee and \x\i pairs. 

The dominant Drell-Yan pp — > Z — > resonance can be easily eliminated by a 
invariant mass cut, with a modest (~ 25 %) loss in top acceptance. After the Z removal, 
the rate of high transverse momentum Drell-Yan pairs is still approximately two orders of 
magnitude higher than the tt dilepton rate (for M top « 150 GeV/c 2 ). Additional background 
rejection can be obtained because (i) in Drell-Yan events there are no additional emitted 
jets at lowest order, and (ii) there are no neutrinos, and hence zero % T except for resolution 
effects. Higher order QCD corrections to the diagram shown in Fig. 37 give rise to final 
state jets, for example from gluons radiated off the q and q lines. For each additional jet, 
the rate is reduced by a factor of 0(a s ) ~ 0.15. Since there are two 6-quark jets in tt 
dilepton events, one can achieve a significant background rejection factor, while maintaining 
efficiency for top, by demanding that at least two jets be detected in addition to the l + l~ 
pair. In conjunction with a E T requirement, the Drell-Yan background can then be reduced 
to a tenth or less of the expected top signal for top masses as high as 200 GeV. 

As was mentioned above, Drell-Yan rr production, followed by leptonic decays of both 
taus (r — > luu), constitutes an additional source of lepton pairs. Events from the Z — > rr 
resonance cannot be easily removed because the invariant mass information is lost due to 
the presence of four neutrinos in the final state. However, the transverse momenta of the 
leptons and the JZ! T for these events are significantly lower than in tt dilepton events, (see 
, and |36|) . By requiring high transverse momentum leptons, high Er , and two 
jets, this background can be reduced to approximately the same level as the ee and \x\i 
Drell-Yan background. We note here that, unlike in the case of direct Drell-Yan production 
of ee and pairs, this process can result in e\i final states. 

Diboson production, (see Fig. 39), constitutes an additional source of high dileptons 
and Er . These are exceedingly rare processes, which are extremely interesting in their own 
right. The cross-section for WW production has been calculated to next-to-leading order, 
and at <Js — 1.8 TeV it is estimated to be a(WW) « 10 pb (Ohnemus, 1991a). This is 
the same as a(tt), for M top « 160 GeV/c 2 . When both PF-bosons decay leptonically, the 
kinematics for the leptons and the neutrinos are very similar as those expected from ti, 
which also results in a WW pair in the final state. The most efficient method that can be 
used suppress this background is to require that there be jets in the event. Just as in the 
Drell-Yan process, there are no jets at leading order in WW events. By demanding that 
there be two jets, the background is reduced by a factor of order a 2 s ~ 0.02. Backgrounds 
from WZ and ZZ production are smaller by over one order of magnitude because (i) the 
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38. The expected lepton and neutrino transverse momenta in pp — 
► l~vv at y/s = 1.8 TeV. From the ISAJET Monte Carlo generator. 
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FIG. 39. Lowest order Feynman diagram for production of WW pairs in pp collisions; similar 
diagrams lead to WZ and ZZ production. 



WZ and ZZ production cross-sections are significantly lower than that of WW (Ohnemus, 
1991b; Ohnemus and Owens, 1991), (ii) the leptonic branching ratios of the Z are a factor 
of three smaller than those of the W, and (iii) l + l~ pairs from Z decays can be eliminated 
with an invariant mass cut. 

Additional backgrounds to the dilepton signal from fake leptons as well as doubly- 
semileptonic decays of bb pairs also have to be considered, but are generally found to be 
small. In a given analysis the signal-to-background level can be tuned by the choice of re- 
quirements. In general, raising the minimum cut on the jets eliminates more background 
events than signal events. The reason for this is that for sufficiently high M top the spec- 
trum of 6-jets in dilepton top events is harder than the brehmstrahlung-like spectrum of jets 
in all the processes listed above. Higher JZl T or lepton P T requirements would considerably 
lower all backgrounds except the diboson background; the Z — > rr background could be 
entirely eliminated by requiring the invariant mass of the l + l~ pair to be higher than the Z 
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FIG. 40. Feynman diagram for W+ 4 jet production in pp collisions. Many other diagrams also 
contribute. 



mass; requiring that jets be 6-tagged would reduce all backgrounds by about two orders of 
magnitude. With sufficient luminosity, it should be possible to obtain very pure ti samples 
in the dilepton mode. 



E. Lepton + jets mode 

The branching ratio for this mode (ti — > qqlvbb) is quite large, 24/81, see Table |TTT| . The 
signature consists of one isolated, high lepton (e or /i), E^ from the neutrino, two light 
quark jets (u, d, c, or s), and two 6-quark jets. 

As was discussed in Section |V.B.2| , in practice the number of detected jets is not expected 
to be always four. In order to maintain high efficiency, in most analyses the number of 
jets requirement is usually relaxed to > 2 or > 3, except for very high top mass, or where 



detection of a fourth jet is essential, e.g. for the determination of the top mass, see Section IX 



1. W+ jets background 

For sufficiently high top mass, the dominant background to ti in the lepton + jets 
channels is due to W + multijets production, (see Fig. The inclusive W production 
cross-section at the Tevatron is <r(pp — > W) ~ 20 nb, over three orders of magnitude higher 
than the ti cross-section for M iop > 150 GeV/c 2 . The cross-section for W+ N jets is 
suppressed by factors of order af. 

Other backgrounds, such as pp — > bb + jets followed by b — > civ, and backgrounds from 
fake leptons also need to be carefully evaluated. For top masses above ~ 40 GeV/c 2 , where 
an efficient lepton + jets top selection can be devised based on both a high lepton and 
high missing transverse energy, these backgrounds are in general found to be much smaller. 

The importance of the W + multijets process as a background to top production, as well 
as to more exotic phenomena (e.g. supersymmetry), was noted soon after the first results 
from the SppS collider became available. As a result, a large theoretical effort was directed 
towards the calculation of the cross-sections and kinematic properties for pp — > W (or Z) 
+ N jets. The first calculations, for N=l or 2, were performed in the mid-eighties (Ellis 
and Gonsalves, 1985; Ellis, Kleiss, and Stirling, 1985; Kleiss and Stirling, 1985; Gunion and 
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Kunstz, 1985). These were then extended to final states with 3 (Berends, Giele, and Kuijf, 
1989; Hagiwara and Zeppenfeld, 1989) and 4 jets (Berends et al, 1991). 

These calculations are performed at tree-level, and therefore they diverge as the angu- 
lar separation between outgoing partons becomes very small or as their transverse momenta 
tend to zero (collinear and infrared divergences). However, far enough away from the regions 
of divergence, calculations are expected to be quite reliable, within the estimated theoret- 
ical uncertainties due to the missing higher order terms. Agreement is found between the 
measured W + jets cross-section and the theoretical prediction, (see Fig. [41]). The degree of 
confidence in these LO calculations is such that the theoretical predictions have been used 
to extract the value of a s from the W+ jet data (Lindgren, 1992; Alitti et al, 1991c). More 
recent measurements of a s (Abachi et al., 1995c) have been based on the full NLO (order 
a s ) calculation for pp — > W (Giele, Glover, and Kosower, 1993). 
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FIG. 41. Product of W cross-section (aw) times leptonic branching ratio as a function of jet 
multiplicity in pp collisions at \J s = 1.8 TeV. The LO QCD predictions are shown for two different 
choices of the renormalization and factorization scale \i. The jet threshold was set at 15 GeV. 
From F. Abe et al, 1993b. 



The LO theoretical calculation for pp — > W (or Z) + N jets, for N up to 4, is implemented 
in the VECBOS Monte Carlo event generator (Berends et al, 1991), which is extensively used 
by the experimenters to model the W + jets background. VECBOS is a parton-level Monte 
Carlo generator, i.e. its output consists of parton four-vectors only. In order to properly 
simulate the response of the detector, hadronization effects are included by interfacing the 
VECBOS event generator with hadronization models based e.g., on independent parton- 
fragmentation (Field and Feynman, 1978), or the HERWIG model (Marchesini and Webber, 
1984 and 1988). A model of the underlying event also needs to be included. 

In what follows we will show several comparisons of theoretical expectations for tt and 
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W+ jets at Tevatron energies. These comparisons are performed at the parton-level only, 
i.e. hadronization, as well as detector effects, such as resolution smearing, efficiencies etc., 
are not included. Therefore, the discussion presented here is intended only as a general 
illustration of the issues involved. 

The two processes are modeled with the ISAJET (for tt) and VECBOS (for W+ jets) 
event generators. The ISAJET model employed here does not include initial and final 
state gluon radiation effects; the tt rates are normalized to a recent pp — > tt cross-section 
calculation (Laenen, Smith, and van Neerven, 1994). In order to mimic actual experimental 
conditions, and to avoid the infrared and collinear divergencies in the W+ jets calculation, 
we impose the following requirements : 

• Pt of partons (quarks or gluons) > 15 GeV/c 

• |^| of jets < 2 

• Pt of leptons (electrons, muons and neutrinos) > 20 GeV/c 

• 1 77 1 of electron or muon < 1 

• AR between jets > 0.5 
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FIG. 42. Comparison of expected lepton + jets rates at the Tevatron for tt and W+ jets as a 
function of jet multiplicity and top mass. Detector and tt gluon radiation effects are not included. 
The W+ jets theoretical predictions are given as bands which reflect the effects of reasonable 
variations in the fi 2 scale. See text for details. 



In Fig. ^2] we show the expected rates of lepton + jets events from tt and W+ jets. (To 
get a feeling for the effects of gluon radiation, which are not included here, see for example 
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Fig. |35|). Note the drop in efficiency for detecting a fourth jet in tt events, even at top 
masses of 200 GeV. As discussed in Section V.B.2| , this is due to (i) final state quarks having 
Pt below the threshold (15 GeV/c in this case), and (ii) to the effect of jet merging. The 
expected signal-to-background, for a lepton + > 3 jets selection, with a 15 GeV jet Pt 
threshold, varies between about 1/1 at M top = 100 GeV/c 2 and 1/10 at M top = 200 GeV/c 2 . 
By requiring > 4 jets, the signal-to-background for high top quark mass is significantly 
improved (by approximately a factor of 3 for M top = 200 GeV/c 2 ). As we will show, further 
improvements in signal-to-background can be achieved by raising the Pt threshold on the 
jets. For top masses in the neighborhood of the W-m&ss, 6-quarks are expected to have a 
soft transverse momentum spectrum, (see Fig. |34|). As a result, the probability of detecting 
more than two jets in this mass region is particularly low. The required number of detected 
jets in a lepton + jets top search is therefore in general dependent on the top mass region 
that is being explored. 

The cross-section for W+ N jets is proportional to af, therefore for each additional jet 
the W cross-section drops by a factor of order a s . The W+ jets predictions in Fig. [42] are 
derived from VECBOS, which is based on a tree-level calculation with significant uncertain- 
ties. These uncertainties can be partially quantified by the stability of the calculation under 
changes in the factorization and renormalization scale fi. Here we present these predictions 
as bands that reflect the variation between the choices fi 2 = M 2 ^ and /i 2 =< Pt > 2 , where 
< Pt > is the average transverse momentum of the partons in the event. 

The relative uncertainty on the W+ N jets cross-section due to the choice of fi grows with 
the number of jets. This is because a(W+ N jets) is proportional to a^(fi 2 ), so that when 
choosing a different fi 2 , 5a jo becomes proportional to N 5a s /a s (where we have ignored 
the fi 2 dependence of the parton distribution functions that need to be convoluted with the 
partonic cross-section). 

Because of these theoretical uncertainties, the existence of the top quark cannot be 
firmly established based on just the observation of an excess of W+ jets events. Additional 
information need to be employed to isolate a potential top signal. A number of possibilities 
will be discussed next. 



2. Separation of W+ jets and tt for M top < M w + M b 

As discussed previously, if M top < M w + M b , the W from the decay t — ► Wb will be 
virtual (W*). In this case, the invariant mass of the Iv pair in tt events from W* — > Iv 
will be smaller than both Mv^ and M top . In contrast, the Iv pair in W+ jets events, which 
originates from real W decays, has invariant mass equal to the W mass. However, as was 
mentioned in Section |V.B.3| , the longitudinal component of the neutrino momentum cannot 
be measured, so the Iv invariant mass cannot be calculated. Fortunately, the transverse 
mass (Mt) of the Iv pair still provides significant discrimination between real and virtual 
H^-decays (Rosner, 1989). The transverse mass is defined as the pseudo-invariant mass of 
the lepton and the neutrino constructed from the transverse components only : 

M 2 T = (P T u + Pti) 2 - (Pt, + Pti? = 2P Tv P Tl {l - cosA0) 
where Pjv and Pti are the neutrino and lepton transverse momenta respectively, is the 
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FIG. 43. Transverse mass distribution in W — > ev decays. From the CDF collaboration, Frisch, 
1995. This distribution is used to measure the W mass by fitting to Monte Carlo expectations; the 
arrows delimit the range of the fit. 

angle between the two transverse momentum vectors, and where we have ignored the mass 
of the lepton compared to its momentum. Transverse mass distributions, (see Fig. ^3|), have 
the following properties : (i) My is less than or equal the invariant mass of the lepton- 
neutrino pair; (ii) My distributions are invariant under longitudinal boosts, hence they are 
independent of the longitudinal momentum of the pair; (iii) they result in Jacobian peaks 
at the original Iv invariant mass; (iv) they are fairly insensitive to the total transverse 
momentum of the pair. (The transverse mass distribution of a pair of invariant mass M 
and transverse momentum differs from that of a = pair by corrections of order w 
(Pt/M) 2 , see for example Barger and Phillips, 1987). 

In Fig. [H] we show expected transverse mass distributions for W+ jets and ti, with 
Mtop = 70 GeV/c 2 . A top signal would result in a significant distortion of the transverse 
mass spectrum of W+ 2 jets events, even after accounting for smearing due to resolution 
effects (see Section |VT| ). By concentrating on just the shape of the transverse mass distri- 
bution, uncertainties due to the theoretical expectation of the W+ jets rate do not enter 
in the analysis. Furthermore, the shape of the transverse mass distribution for the W+ 
jets background is dominated by the kinematics of the W — > Iv decay, and depends only 
weakly on the modeling of the W production properties. We note that the transverse mass 
method can be used to separate a top signal from the W+ jets background for both the ti 
and W — » tb production mechanisms. 

For higher top masses, the Iv transverse mass distributions in ti and W+ jets events 
become indistinguishable, since top quarks will decay into real W bosons. In order to 
separate signal from background in the lepton + jets mode, one then has to rely on the 
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Transverse Mass (GeV/c 2 ) 

FIG. 44. Expected transverse mass distributions for W+ 2 jets (solid) and pp — ► ti — > lepton + 
2 jets, for M top = 70 GeV/c 2 (dashed). This is at Tevatron energies (1.8 TeV). The W+ 2 jets 
calculation has been performed with a scale fi 2 =< Pt > 2 . Detector and it gluon radiation effects 
are not included. Relative normalizations are from Monte Carlo. See text for details. 

kinematic differences between the two processes, and/or the fact that ti events always contain 
two 6-quarks in the final state. 

3. Kinematic differences between ti and W+ jets 

Several possible ways of extracting a top signal from the W+ jets background using 
kinematic signatures have been suggested in the literature (Berends et al, 1989; Baer, Barger 
and Phillips, 1989; Agrawal and Ellis, 1989; Giele and Stirling, 1990; Berends, Tausk, and 
Giele, 1993; Barger, Ohnemus, and Phillips, 1993; Benlloch, Sumorok, and Giele, 1994; 
Cobal, Grassmann, and Leone, 1994; Barger et al, 1995). Briefly, the differences between 
the two processes are the following : 

• In ti events, the invariant mass of three of the jets should reconstruct to the top mass, 
and two out of these three jets should have invariant mass equal to the W mass. This 
is the consequence of the decay chain t — > Wb followed by W — > qq. Of course no such 
invariant mass enhancements occurs for the W+ jets background. 

• Jets in W+ jets events tend to have lower transverse momentum than jets in ti events. 
This is due to the fact that jets in W events arise from a brehmstrahlung-like process. 
For the same reason, these jets also tend to be emitted more in the forward direction 
than jets from the decay of centrally produced top quarks. 
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• Top events tend to be more spherical and aplanar than W+ jets events. The reason for 
this is that the QCD W+ jets matrix element introduces significant spatial correlations 
between the jets, e.g. gluon brehmstrahlung (q — > qg) and gluon splitting (g — > qq) 
favor small opening angles between partons in the final state. 



At first glance two-jet and three-jet invariant masses appear to be the most attractive 
discriminators. Unfortunately, because of the poor jet energy resolution and the number of 
possible jet combinations that may be present in a given event, this method turns out to 
be useful only if very large statistics data sets are available. We will not further discuss jet 
invariant masses in this Section. We will however revisit the issue in Section IX, where the 
measurements of the top mass will be reviewed. 

As an illustration of the kinematic differences between W and top events at the Tevatron, 
we show the Monte Carlo transverse momentum and pseudorapidity distributions of jets in 
lepton + 4 jets events for W+ 4 jets and ti (see Figs. [45], |46|, and |47|). As anticipated, jets 
in top events tend to be more central and to have higher transverse momenta. The scalar 
sum of the transverse momenta of all of the jets in the event is a simple global variable that 
is expected to provide good discrimination between signal and background, (see Fig. flop . 

Note that the choice of scale in the W+ jets calculation affects not only the expected 
rate, but also the shape of kinematic distributions, particularly the jet transverse momentum 
spectrum. Comparing the distributions for the two equally arbitrary, and a-priori equally 
reasonable choices, /i 2 =< Pt > 2 and /i 2 = M 2 ^, it is found that the /i 2 = M 2 ^ choice results 
in a harder jet Pt spectrum, see Figs. ^5| and |6. This can simply be understood as follows. 
The cross-section for a W+ N jets event is just given by the convolution of the relevant 
matrix element with the parton distribution functions. Neglecting the /^-dependence of the 
parton distribution functions, the only scale dependence is due to the factor a^(fi 2 ) which 
appears in the W+ N jets matrix element. Because of the running of the strong coupling 
constant, the choice of an event-by-event scale such as /i 2 = < Pt > 2 results in a higher 
probability for events with low Pt jets as compared to what one would obtain by choosing 
a global scale like fi 2 = M 2 ^. 
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FIG. 45. Expected transverse momentum (Pt) distributions of the first, second, third, and fourth 
highest Pt jet in W+ 4 jets and tt events. For pp collisions at s = 1.8 TeV. Solid line : W+ 



4 jets, scale /i 2 =< Pt > 2 



dashed line 



W+ 4 jets, scale /x 2 = M 2 ^ 



dots : tt, M top = 170 



GeV/c 2 . Note the expanded horizontal scale for the third and fourth jets and that all jets have 
been required to have Ey > 15 GeV. Detector and tt gluon radiation effects are not included. See 
text for details. 
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FIG. 46. Expected scalar sum of transverse momenta of 4 jets in W+ 4 jets and ti events. For pp 
collisions at y/s = 1.8 TeV. Solid line : W+ 4 jets, scale /x 2 =< Pt > 2 ; dashed line : W+ 4 jets, 
scale fi 2 = Myy ; dots : ti, M top = 170 GeV/c 2 . Detector and ti gluon radiation effects are not 
included. See text for details. 
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FIG. 47. Pseudorapidity distribution of jets in W+ 4 jets and tt events. For pp collisions at ^/s = 
1.8 TeV. Solid line : W+ 4 jets, scale p 2 =< P T > 2 ; dashed line : ti, M top = 170 GeV/c 2 . 
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Another discriminant that can be employed is aplanarity. Aplanarity is defined as A 
= 3/2Ai, where Ai is the smallest eigenvalue of the matrix M ab = J^PaPb/ where 
Pi are the cartesian components of momentum of the parton, P is the magnitude of the 
three-momentum, and the sum is over all final state objects: jets, electrons, muons, and 
neutrinos. Zero aplanarity corresponds to planar events. W+ jets events are expected to be 
more planar than ti events, (see Fig. [18|). 
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FIG. 48. Expected aplanarity for W+ 4 jets and ti events. For pp collisions at y/s = 1.8 TeV. 
Solid line : W+ 4 jets, scale /j, 2 =< Pt > 2 ; dashed line : W + 4 jets, scale [i 2 = Mf^ ; dots : ti, 
Mtop = 170 GeV/c 2 . Aplanarity here is calculated in the laboratory frame from the four jets, the 
lepton, and the neutrino. Detector and ti gluon radiation effects are not included. See text for 
details. 

It is clear from this discussion that in order to extract a top signal from kinematic 
distributions, special attention must be payed to the systematic uncertainties associated 
with the theoretical modeling of the W+ jets backgrounds. The theoretical uncertainties 
can be bounded by comparing the data with theoretical predictions in kinematic regions 
where the top signal is small, e.g. in samples of events with low jet multiplicity and/or low 
transverse momentum jets. Deviations from the W+ jets expectations in kinematic regions 
where the top quark is expected to contribute would then signal the presence of ti events 
in the sample (or possibly of some other source of W+ jets events beyond standard QCD 
production). 

Samples of Z+ jets events in principle provide the ideal testing ground for the W+ jets 
calculation. Unfortunately, at y/s = 1.8 TeV the cross-section for pp — > Z is a factor of 3.3 
smaller than that for pp — > W. Furthermore, the leptonic branching ratio of the Z is also 
a factor of three smaller than that of the W. As a result, the number of reconstructed Z 
events is approximately one order of magnitude smaller than that of W events. Because 
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of the limited statistics, Z events, although useful as a first order check, do not provide 
stringent bounds on the modeling of vector boson + > 3 jets production. We will discuss 
these issues in more detail in Section IVTTTl where the experimental results will be reviewed. 



4. 6-quark tagging 

An alternative method that can be employed to separate the ti signal from the W+ 
jets background is to tag the 6-quarks in top events. Each ti event contains one b- and one 
6-quark from the decays t — > Wb and i —>■ Wb, whereas the jets in W events arise mostly 
from the fragmentation of gluons and light quarks. 

There are two ways to detect the presence of 6-quarks. The first method is based on 
the detection of additional leptons from the semileptonic decays b — > civ or b — > c —>■ slv. 
The semileptonic branching ratios of bottom and charm quarks are approximately 10% per 
lepton species (Montanet et ai, 1994). There is on average about one lepton (electron + 
muon) from b or c decay in each top event, where we have also included the contributions 
from c-quarks from W — > cs, which occurs in one half of all hadronic W decays. From 
an experimental point of view, detection of these leptons is more difficult than detection 
of leptons from W decays, because these leptons tend to have a much lower transverse 
momentum (compare Figs. |49| and ^9]). Furthermore, these leptons are not isolated but are 



accompanied by nearby hadrons from the 6-quark fragmentation and the 6-hadron decay. 
This makes efficient detection of electrons particularly challenging. 
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FIG. 49. (a) Expected transverse momentum of leptons from t — > b — > civ Solid line : M.t op = 110 
GeV/c 2 ; dashed line : M t0 p = 170 GeV/c 2 . (b) Same as (a), but for t — > b — > c — > slv. For pp — > tt 
at -y/s = 1.8 TeV, from the ISAJET Monte Carlo generator. Note the different horizontal scales in 
(a) and (b). 
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Bottom and charm quarks can also be tagged by exploiting the long lifetime of b- and 
c-hadrons. The recent compilation from the Particle Data Group (Montanet et al.) reports 
a lifetime of 1.537 ± 0.021 ps for 6-hadrons, and 0.415 ± 0.004 ps and 1.057 ± 0.015 ps 
for neutral and charged D mesons respectively. As a consequence of the long lifetime, b- 
hadrons in a top event are expected to travel several mm before decaying, (see Fig. |50"| ). 
With the advent of silicon microstrip-vertex detectors, the position of decay vertices can 
be measured with resolutions of order 100-150 /im. It then becomes possible to separate 
with good efficiency the secondary vertex where the 6-decay occurs from the primary pp 
interaction point. 

There are two main sources of background to the lepton- or vertex-tagged ti signal. 
The first one is instrumental. Hadrons originating from the fragmentation of gluons and 
light quarks in W+ jets events can be misidentified as muons or electrons. This happens for 
example when hadronic showers in the calorimeter fluctuate to mimic the electron signature, 
or when kaons and pions decay to muons. Also, track mismeasurements or decays of other 
long-lived particles such as A and K s can result in the reconstruction of spurious detached 
vertices. Needless to say, these effects have to be considered very carefully. As we will discuss 
in Section |VII1| , sufficiently good instrumental background rejection has been achieved, and 
methods to precisely estimate the remaining background have been developed. 
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FIG. 50. Distance travelled by 6-hadrons in top events before decaying. Solid line : M top = 
110 GeV/c 2 ; dashed line : M top = 170 GeV/c 2 . This distance is calculated by convoluting the 
momentum spectrum of 6-hadrons from the ISAJET pp — ► ti (y/s = 1.8 TeV) Monte Carlo with 
their lifetime. 

The second background source stems from the fact that a small fraction of jets in W+ jets 
events will contain heavy-quarks (b or c). In the absence of additional kinematic information, 
these events constitute an irreducible physical background to the lepton + jets + fe-tag 
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FIG. 51. Lowest order Feynman diagram for production of W+ charm in pp collisions. At higher 
order, gg — > Wcs diagrams also contribute. 
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FIG. 52. Lowest order Feynman diagram for production of WQQ in pp collisions. Q = c or b. 



ti signature which needs to be carefully evaluated. A comprehensive discussion of these 
backgrounds, in the context of the CDF experiment, is given by F. Abe et ai, 1994a, and 
will be summarized below. The discussion here applies equally well to lepton or vertex 
6-tagging methods. 

Heavy-quarks in W+ jets event can be produced singly, in the process sg — > Wc or 
dg — > Wc, (see Fig. 51), or in pairs, when a gluon in a W+ jets event splits into a cc or bb 
pair, (see Fig. B2t). The gluon splitting probability, g — > QQ, Q = c or b, is estimated to be 



of order a few percent, with significant theoretical uncertainties (Mangano and Nason, 1992; 
Muller and Nason, 1985 and 1986). Experimental extractions of the g — > cc probability in 
both pp (Ikeda, 1990, F. Abe et ai, 1990d) and e + e _ collisions (Akers, 1995) are found to be 
in agreement with the results of the theoretical calculations, within the large experimental 
and theoretical uncertainties. In what follows we will refer to these two processes as Wc 
and WQQ respectively. 

The Wc background estimation begins by computing the fraction of W+ jets events that 



contain a single c-quark, from diagrams like the one shown in Fig. |51], using the VECBOS and 
HERWIG event generators. This fraction is found to be of order 8%, with small variations 
depending on the jet multiplicity and the choice of the input parton distribution function 
for strange quarks in the proton. To obtain an absolute background prediction, this fraction 
is then multiplied by the number of observed W+ jets event and the tagging efficiency for 
c-quarks. The tagging efficiency for these events clearly depends on the details of the tagging 
algorithm. However, it is in general much smaller (typically by a factor of order 5) than 
the tagging efficiency for ti since (i) the tagging efficiency for c-quarks is lower than that of 
6-quarks because e.g. there are fewer tracks in a c decay than in a b decay, and (ii) there are 
multiple b- and c-quarks in a top event that can potentially be tagged. Therefore, as a result 
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FIG. 53. Two of the possible LO Feynman diagrams for jet production in pp collisions: gg — > gg 
and gg -> gg. 

of the small probability for a Ifc event and its small tagging efficiency, the Wc background 
is much smaller than the expected tagged ti rate. The highly uncertain overall normalization 
of the W+ jets theoretical calculation does not enter in the background estimate, since only 
the fraction of Wc events is taken from theory. 

The WQQ background is more important because the tagging efficiency for Wbb is 
comparable to that of ti. This background can be estimated in two ways. The first method 
(Method T) requires a minimum of theoretical input, and is expected to yield an overestimate 
of the WQQ background. The alternative method, (Method IT), is based on the state-of- 
the-art theoretical understanding of heavy quark production, 

At the heart of the Method I background calculation is the assumption that the heavy 
flavor content (b or c) of jets in pp — > jets (generic-jets, see Fig. |53|) is the same or larger 
than the heavy flavor content of jets in W+ jets events. Accepting this assumption for now, 
we proceed to describe the Method I background calculation. 

The generic-jet sample includes gluon and light-quark jets, as well as a small fraction of 
heavy-quark jets. A generic-jet can be tagged due to instrumental effects resulting in a false 
tag of a light-quark or gluon jet, or due to the b- or c-quark contribution. Operationally, a 
tag-rate is measured for generic-jets as the probability of tagging a generic-jet as a function 
of several relevant variables (e.g. jet Pt, track multiplicity). The generic-jet tag-rate is then 
applied to the sample of jets in W+ jets events, to predict an upper limit for the sum of 
the instrumental and WQQ backgrounds to the ti — ► W+ jets + 6-tag signature. From an 
experimental point of view this has the advantage that both the instrumental and WQQ 
backgrounds are estimated simultaneously directly from the data. No a-priori knowledge of 
the tagging efficiency, or the WQQ content of the sample, is needed. 

We now turn to a discussion of the theoretical assumption on which the Method I WQQ 
background calculation is based. At Tevatron energies, generic-jets in the relevant range 
(20-150 GeV/c) consist predominantly of gluon-jets. There will also of course be a con- 
tribution from light-quark jets as well as pp — ► QQ. The lowest order Feynman diagrams 
for direct production of bb and cc pairs in pp collisions are identical to the ones for ti pro- 
duction shown in Fig. 16. For the b- and c-quarks however, it is found that higher order 
diagrams such as the ones displayed in Fig. |53| contribute a very significant amount to the 
pp — > QQ + X cross-section. For instance, the gluon splitting diagram is believed to account 
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FIG. 54. Higher order Feynman diagrams for QQ production in pp collisions: (a) gluon splitting; 
(b) flavor excitation. 

for of order 70% of the bb and cc production rate (Mangano, Nason, and Ridolfi, 1992). 

The heavy-flavor content of generic-jets is expected to be higher than that of jets in 
W+ jets events for two reasons : (i) in W+ jets events there are no contributions from 
lowest order, direct, QQ production (Fig. 16) and next-to-leading order flavor excitation 
(Fig. |54|b), and (ii) generic-jets consist mostly of gluon jets which can result in QQ pairs 
via gluon splitting, whereas the jets in W+ jets events consist of an approximately equal 
mixture of gluon and light-quark jets. These very simple qualitative arguments are borne 
out by a more quantitative calculation, as we will illustrate below. 

The second way of estimating the WQQ background (Method II) is based on an explicit 
calculation of the WQQ process. Just as in the Wc case, uncertainties on the overall 
normalization of the W+ jets calculation are minimized by using as the theoretical input the 
fraction of W+ jets events that contain a QQ pair, rather than the absolute rate prediction. 
The absolute background level is then estimated by multiplying this fraction by the number 
of observed W+ jets events and the tagging efficiencies for Wcc and Wbb. 

The WQQ fraction is estimated using a combination of the HERWIG Monte Carlo event 



generator and the lowest order WQQ matrix element (Fig. [52]) calculation (Mangano, 1993). 
This calculation differs from the W+ jets calculation in that it includes all mass effects and is 
therefore free from collinear and infrared divergencies. Since the calculation does not include 
higher order terms, it cannot be used directly to estimate the WQQ rate in events with more 
than two jets, which are the most relevant for the top search. Instead, the results of the 
exact calculation are compared to those of HERWIG, where the QQ pairs are produced by 
gluon splitting from initial- and final-state parton evolution. The HERWIG results are found 
to be in good agreement with those of the WQQ matrix element calculation. The gluon 
splitting process in HERWIG also provides a good description of the measured tagging rate 
in generic-jet events, giving further evidence for the validity of the HERWIG model. Thus, 
this model is used to predict the fraction of WQQ events as a function of jet multiplicity. 
For the sample of W+ > 3 jets, it is estimated that the fractions of W events containing 
a bb or cc pair are approximately 3% and 5% respectively. The systematic uncertainties 
associated with these predictions are clearly significant and are estimated to be at the level 
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of 80%. The WQQ background estimated with this method turns out to be approximately 
a factor of three lower than the conservative estimate made under the assumption that the 
heavy flavor content of jets in W+ jets events is the same as that of generic-jets. 

To assess the impact of the WQQ background on the top search, we re-examine the 
predicted rates of lepton + jets from ti and W+ jets. For example, from Fig. [E^, the signal- 
to-background for a W+ 4 jets selection, with jet Pt threshold set at 15 GeV/c, is expected 
to vary from approximately 3/1 at M top = 120 GeV/c 2 to approximately 1/3 at M top = 
200 GeV/c 2 . With the fraction of W+ jets events containing a QQ pair given above, the 
signal-to-background, (ti vs. WQQ) will be at least of the order of 7/1. Thus, provided the 
instrumental backgrounds can be adequately controlled, and enough luminosity is available 
to the experimenters, the top signal can be separated from the W+ jets background with 
6-tagging methods. 

VI. EARLY SEARCHES FOR THE TOP QUARK 

Experimental searches for the top quark began immediately following the discovery of 
the 6-quark. In this Section we will review these searches up to approximately 1990. 

A. Searches in e + e~ collisions 

In e + e~ collisions, top quarks would be produced in pairs through e + e~ annihilation into 
a photon or a Z. Since at leading order this is a purely electroweak process, the production 
cross-section can be accurately calculated. At center-of-mass energies well below the Z mass, 
where annihilation of the e + e~ pair into a photon dominates, ti production would manifest 
itself as an increase by an amount <5R ~ 3Q 2 op = 4/3 in the ratio R = a(e + e~ —>■ hadrons) 
/cr(e + e~ — > well above the energy threshold for the production of a ti pair. This 

expected increase in R is largely independent of the top quark decay mode, as long as the 
top-quark decays into final states containing hadrons. 

Between 1979 and 1984, measurements of R were performed at the PETRA e + e~ collider 
in the center-of-mass energy range between 12 and 46.8 GeV (Barber et al, 1979 and 1980; 
Berger et a/., 1979 and 1980; Bartel et al, 1979a, 1979b, and 1981; Brandelik et al, 1982; 
Adeva et al, 1983a, 1983b, 1985, and 1986; Behrend et al, 1984; Althoff et al, 1984a 
and 1984b). The value of R was found to be consistent with Standard Model expectations 
without a top-quark contribution. Event topology studies gave no evidence for excesses 
of spherical, aplanar, or low thrust events that could be attributed to ti production. The 
measured rate of prompt muons was also found to be in agreement with expectations from 
models of c- and 6-quark production and decay, and could not accommodate a contribution 
from semileptonic decays of top quarks. Existence of the top quark with mass below 23.3 
GeV/c 2 was ruled out at the 95% confidence level (C.L.). 

Similar searches were later performed at the TRISTAN collider, which reached a center- 
of-mass energy of 61.4 GeV (Yoshida et al, 1987; Sagawa et al, 1988; Adachi et al, 1988; 
Igarashi et a/., 1988; K. Abe et a/., 1990). No evidence for top quark production was 
reported, resulting in a lower limit on the top mass of 30.2 GeV/c 2 . 
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In 1989-90, the SLC and LEP e + e~ colliders, with y^s pa M z became operational. Studies 
of event topologies and measurements of the Z width were found to be inconsistent with 
a Z — > ti contribution, and resulted in lower limits on the top quark mass as high as 45.8 
GeV/c 2 (Abrams et al, 1989; Akrawi et al, 1990a; Decamp et al, 1990; Abreu et al, 1990b 
and 1991; Adriani et al, 1993). 



B. Early searches in pp collisions assuming Standard Model top-quark decay 

With the emergence in the 1980's of pp colliders, first at CERN and then at Fermilab, 
and with evidence from e + e~ experiments pointing towards a very high mass for the top 
quark, focus in the search for top rapidly shifted to hadron colliders. The obvious advantage 
of hadron colliders for top physics is the high center of mass energy, which enables the 
exploration of higher mass regions. However, unlike in e + e~ collisions, the large backgrounds 
make it impossible to directly search for the top quark in a model-independent way. It is 
necessary to concentrate on particular signatures, based for example on the Standard Model 
decay modes of the top quark discussed in the previous Section. 

In this Section we will discuss searches for the top quark in pp collisions, assuming 
Standard Model top-quark decay, that were carried out in the mid to late 1980 's at CERN 
and Fermilab. No evidence for top quark production was uncovered, leading to lower limits 
on the top quark mass as high as 91 GeV/c 2 , at the 95% C.L.. These searches do not 
result directly in limits on the top mass, but rather in upper limits on the product of top 
production cross-section and branching ratio for top quark decay. To turn these limits into 
mass limits for the top quark, it is necessary to (i) assume that the top quark decays as 
t — > Wb, as prescribed by the Standard Model, and (ii) use theoretical expectations for 
the production cross-section. As discussed in Section |TV], the main production mechanisms 
are W — > tb decays at SppS energies, and ti pair production at Tevatron energies. The 
expected cross-section for pp — > W — *> tb can be reliably predicted from direct measurements 
of pp — > W — > lu; on the other hand, there are significant theoretical uncertainties on the 
predicted top-pair production cross-section as a function of the top mass, see Section [TV], 



Figs. ^T] and |22j. Since these theoretical uncertainties are difficult to quantify, limits on 
M top in the absence of a top signal in the data are placed based on the lower range of the 
calculation of the top production cross-section. This results in conservative 95% CL. lower 
limits on M fop . 

Initial results reported by the UAl collaboration (Arnison et al, 1984; Revol, 1985) at 
the CERN SppS collider (-y/s = 630 GeV/c), seemed to be consistent with production of 
a top quark of mass 40 ± 10 GeV/c 2 . These results were based on the observation of 12 
isolated lepton + 2 jets events, with an expected background of approximately 3.5 events 
in an exposure with an integrated luminosity of 200 nb _1 (Revol, 1985). In these events, 
the invariant mass of the lepton, neutrino and one of the jets was found to cluster around 
a common value of approximately 40 GeV/c 2 , while the invariant mass of the two jets, the 
lepton, and the neutrino was found to be consistent with the W mass. This is the expected 
signature for the process pp — ► W — > tb, followed by t — > blv. The excess of events over the 
background prediction was also consistent with the expected top production cross-section, 
providing further evidence for a top quark with M top pa 40 GeV/c 2 . 
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These first results were however not supported by a subsequent UAl analysis (Albajar et 
al, 1988), with a higher statistics data sample, as well as a more complete evaluation of the 
backgrounds. This analysis was based on samples of events with one isolated muon + > 2 
jets or one isolated electron + > 1 jet. The integrated luminosity was 700 nb _1 . Given this 
integrated luminosity, and the expected top production rate, this search was sensitive to a 
top quark with M top < 55 GeV/c 2 . A maximum lepton- neutrino transverse mass requirement 
of 40 (45) GeV/c 2 was imposed in the muon (electron) sample to substantially reduce the 
W — > Iv background, while maintaining good efficiency for top quarks in the relevant mass 
range. No missing transverse energy (JZl T ) requirement was imposed, since for low top mass 
the probability for both the lepton and the neutrino in t — > blu to have high transverse 
momentum is low. 

With no Ex requirement, the backgrounds from fake leptons, Drell-Yan l + l~ pairs, as 
well as bb and cc production are significant. The number of observed events was found to be 
fully consistent with the non-top background contribution only. Furthermore, it was shown 
that the event selection requirements resulted in background invariant mass distributions 
that were similar to those expected from top-quark production and decay (for M top ~ 40 
GeV/c 2 ). Thus, this property of the background accounted for the features of the invariant 
mass distributions observed in the previous analysis. Dilepton events were also studied, and 
their properties were found to be in agreement with expectations from semileptonic decays 
in bb and cc events. 

An upper limit on the combined W — * tb and ti cross-sections was extracted from the 
background-subtracted lepton + jets event rates, as well as kinematic distributions such as 
those of jet Et, $ t , and lepton isolation. At the time that these results were obtained, only 
a LO calculation for pp — > ti was available. Using a conservative value for this cross-section 
the limit was inferred to be M top > 44 GeV/c 2 at the 95% confidence level. A subsequent 
re-evaluation (Altarelli et al, 1988), based on the NLO calculation of heavy-quark pair- 
production (Nason, Dawson, and Ellis, 1988) resulted in a slightly modified lower limit 
> 41 GeV/c 2 . This limit is obtained from the ti cross-section corresponding to the 
lower range of the theoretical prediction. 

The lesson to be learned from the UAl experience is that the reliability of background 
estimates is of paramount importance. Their first analysis did not include the J/^, T, and 
Drell-Yan backgrounds, and the bb + cc backgrounds were underestimated by a factor of four. 
The top signature at hadron colliders is complicated, and involves comparing the number 
of observed events and/or kinematic distributions with background expectations. If at all 
possible, dependences on uncertain theoretical models of background processes should be 
minimized. 

More sensitive searches for the top quark were performed in the period 1988-89. At 
CERN, the new Antiproton Accumulator Complex (AAC) was commissioned, resulting in 
luminosities as high as 3 x 10 30 cm~ 2 s" 1 for the SppS. The UAl electromagnetic calorimeter 
was removed to allow for a replacement based on Uranium-TMP (tetramethylpentane) tech- 
nology. This new calorimeter was unfortunately not ready to be installed, resulting in the 
loss of electron identification in UAl. The UA2 detector at the CERN SppS was significantly 
upgraded, with improved calorimeter coverage, and enhanced electron detection capabilities. 
In the US, the first high-luminosity (up to 2 x 10 30 cm~ 2 s^ 1 ) run of the Tevatron Collider 
also started in 1988, with the CDF detector ready for data-taking. 
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The 1988-89 UA1 top search (Albajar et ai, 1990) was based on an integrated luminosity 
of 4.7 pb _1 , a factor of seven higher than had been previously available. Because of the 
missing EM calorimeter, only fi + > 2 jets and /x/i final states were considered. The /i + 
jets analysis was essentially an extension of the previous UAl search (Albajar et ai, 1988), 
with the maximum transverse mass requirement raised to 60 GeV/c 2 . For each event, a 
likelihood-like variable (C\) was defined to discriminate between top and background on 
an average basis. This likelihood was based on the Et, the isolation of the muon, and the 
opening angle between the muon and the highest E^ jet in the event. 
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FIG. 55. The ln(£i) distribution compared with the expected background and top contributions. 
The shaded histogram is for the K and ir — > [iv background. The other major background con- 
tribution is from bb and cc events. The expected top contribution, scaled up by a factor of 10, is 
shown in the dashed histogram. From the UAl collaboration, Albajar et al, 1990. 



The C\ distributions of data, expected background, and expected signal, are shown in 
Fig. [5^. The data are well- described by the background contribution only. Based on the 
number of observed events with ln(£i) > 4, and the expected top production cross-section, 
the fi+ jets data from UAl results in a lower limit on the mass of the top quark of 52 GeV/c 2 
at the 95% confidence level. 

A similar likelihood variable was defined for /z/i events, based on the transverse momen- 
tum and isolation of the highest muon, and the azimuthal opening angle between the 
two muons. Recall that at yfs = 0.63 TeV for M top > 40 GeV/c 2 the major source of top 
events is W decays, W — > tb, see Section [TV|. Hence, the dilepton final state arises from 
semileptonic decays of both the t- and 6-quarks, and the muon from b —> c\iv is not expected 
to be isolated. Studies of the likelihood distribution for /x/i events also gave a null result. 
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Results from the /i+ jets and the data were combined with the lower statistics earlier 
UA1 results, yielding a lower limit on the top mass of 60 GeV/c 2 (95% C.L.). 

A better limit on the top-quark mass was also obtained at the same time by the UA2 
collaboration, based on an integrated luminosity of 7.5 pb -1 at the SppS (Akesson et al, 
1990). Without muon detection capabilities, the UA2 results were based entirely on the 
electron + > 1 jet channel. The UA2 search strategy differed considerably from that of 
UAl, with a Fj t > 15 GeV and an electron isolation requirement imposed in the event 
selection. As a result, backgrounds from bb, cc, and fake electrons were highly suppressed, 
and contributed only of order 10% to the data sample. The bulk of the background was due 
to W — > ez/+ jets events. As was discussed in Section |V.E.2| , for M top < M w + M b , the W in 
the decay t — > Wb is virtual. The differences between the lepton-neutrino transverse mass 
distributions of real and virtual W — > Iv can then be exploited to separate a top signal from 
the background. 

In Fig. we show the transverse mass distribution for electron + jet + FJ T events in 
UA2. No evidence for an excess of low transverse mass events was found, resulting in a 
lower limit of M top > 69 GeV/c 2 . 
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FIG. 56. The electron-neutrino transverse mass distribution in electron + jet + E!t events. The 
data are inconsistent with a top contribution. From the UA2 collaboration, Akesson et al, 1990. 

The first high statistics run of the Tevatron Collider also took place in 1988-89, with an 
integrated luminosity of 4 pb _1 recorded by the CDF collaboration. As was discussed in 
Section [TV|, at the center-of-mass energy of the Tevatron (-^/s = 1800 GeV), top production 
in the relevant M top range is dominated by the pp — > tt process. The CDF collaboration 
searched for top in both the lepton + jets (see Section |V.E| ) and dilepton (see Section |V.D|) 
modes. 
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The CDF top search in the isolated electron + > 2 jets channel (F. Abe et ai, 1990a and 
1991a) was qualitatively similar to the UA2 search. An explicit FJ T requirement was imposed, 
yielding a data sample containing W — > ev + 2 jets events, with a small contamination from 
semileptonic b- and c-quark decays, as well as fake electrons. The resulting transverse mass 
distribution was found to be consistent with no top contribution, (see Fig. pTP - The existence 
of a Standard Model top quark with 40 GeV/c 2 < M top < 77 GeV/c 2 was ruled out. (The 
limit was not extended below 40 GeV/c 2 because of poor acceptance at low top quark mass). 
Consistent results were also found in a subsequent study of fj,+ jets events (Demortier, 1991). 
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FIG. 57. The transverse mass distribution of electron + two or more jets + Ft events. From the 
CDF collaboration, (F. Abe et ai, 1990a). The solid and dashed lines represent expectations from 
the W + 2 jets and tt (M top = 70 GeV/c 2 ) Monte Carlo calculations. 



Both the UA2 and CDF results were based on studies of the shape of the transverse mass 
distribution in lepton + jets events and were therefore independent of uncertainties in the 
theoretical predictions for the pp — > W+ jets cross-section. Furthermore, the shape of the 
transverse mass distribution depends only weakly on the details of W+ jets production. It 
is sensitive mostly to the kinematics of the W decay, and the % T resolution in the detector, 
and it is fairly insensitive to the transverse momentum of the W, see the discussion in 
Section [V.E.2| . Therefore, the influence of the theoretical modelling of the background was 
minimized in the CDF and UA2 lepton + jets top searches. 

The CDF collaboration also searched for tt production in the dilepton channel. Initially, 
only the e// channel was considered (F. Abe et ai, 1990b). As discussed in Section |V.D| , the 
tt dilepton signature consists in principle of two isolated high P T leptons, FJ T and two jets. 
Since the backgrounds in this channel are small, in order to maximize the top acceptance, no 
J£ T , isolation, or jet requirements were imposed. One event was observed, with an expected 
background of 1.4 events, mostly from Z — > tt followed by leptonic decays of both taus. An 
upper limit on the tt production cross-sections was then obtained based on the observation 
of one event, under the conservative assumption that this one event was due to tt production 
and decay. The mass region 28 < M top < 72 GeV/c 2 was excluded at the 95% C.L. 

A subsequent CDF search (F. Abe et ai, 1992a) in the dilepton channel, based on the 
same integrated luminosity of 4 pb _1 , also included the ee and \i\x channels, resulting in an 
increase in the top acceptance by a factor of two. As discussed in Section |V.D| , backgrounds 
in the ee and \i\x channels are in general higher than in the efi channel. Z —>■ ee and 
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Experiment Integ. luminosity 


Mode 


Mass limit (95% C.L.) 


Ref. 



UA1 


5.4 pb- 1 


H + jets, 


> 52 GeV/c 2 


UA1 


5.4 pb- 1 


fj, + jets, nfj, 


> 60 GeV/c 2 


UA2 


7.5 pb- 1 


e + jets 


> 69 GeV/c 2 


CDF 


4 pb" 1 


e/i 


> 72 GeV/c 2 


CDF 


4 pb- 1 


e+ jets 


> 77 GeV/c 2 


CDF 


4 pb" 1 


dileptons (ee, mm, efi) 


> 85 GeV/c 2 


CDF 


4 pb" 1 


dileptons and 1 + jets +6— tag 


> 91 GeV/c 2 



Albajar et al, 1988 
Albajar et al, 1988 
Akesson ei a/., 1990 
F. Abe et al, 1990b 
F. Abe et al, 1990a 
F. Abe et al, 1992a 
F. Abe et al, 1992a 



TABLE VI. Summary of lower limits on the top quark mass from pp collisions, circa 1992. See 
text for details. 



Z — > /i/i decays were eliminated by an invariant mass cut, leaving a large number of events 
from off-shell Drell-Yan production of ee and \x\x pairs. In order to control this background, 
additional requirements were imposed on the JZl T and the azimuthal opening angle between 
the two leptons. 

No events consistent with ti were found in the ee and /i/i channels. The resulting limit 
from the dilepton (ee + /i/i + e/i) channel was M top > 85 GeV/c 2 , at the 95% confidence 
level. 

The CDF experiment also searched for top-quarks in the electron or muon + > 2 jets 
channel, where the dominant W+ jets background was reduced by attempting to tag b- 
quarks through their semileptonic decay into muons, b — > /j or b — > c — > /j (F. Abe et a/., 
1992a). No events consistent with the top hypothesis were found. Despite the high branching 
ratio for the lepton + jets mode, the acceptance in this search was approximately a factor 
of 3 smaller than that of the search in the dilepton channel for M top ~ 90 GeV/c 2 . The low 
acceptance was the consequence of the low (~ 4.5%) muon tagging efficiency for lepton + 
jets top events. This low tagging efficiency was due to (i) the semileptonic 5-quark branching 
ratio, (ii) the limited (\rj\ < 0.6) muon coverage of the CDF detector as configured for the 
1988-89 run, and (iii) the low transverse momentum of 6-quarks in top decays for M top ~ 90 
GeV/c 2 . Combining results from this search and the dilepton search, resulted in a 95% C.L. 
lower limit on the top quark of 91 GeV/c 2 , (see Fig. ^). This limit was extracted using the 
NLO calculation of the ti cross-section (Ellis, 1991). The 91 GeV/c 2 limit corresponds to 
the point where the cross-section limit curve crosses the lower (i.e. more pessimistic) bound 
of the theoretical prediction. Using a more up-to-date calculation of a(ti) (Laenen, Smith, 
and van Neerven, 1994), would result in a limit of 95 GeV/c 2 . 

The lower limits on the top mass from the UAl, UA2, and CDF experiments are sum- 
marized in Table 0. With the top quark being so massive, hopes to observe the top quark 
at the CERN SppS were abandoned, because of the small top production cross-section at 
■s/s = 630, (see Fig. ^). In Section |VTII| we will discuss the most recent, higher statistics, 
searches for the top quark at the Tevatron Collider, which finally resulted in the discovery 
of the top quark. 
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FIG. 58. The 95% C.L. upper limits on the tt cross-section from F. Abe et al, 1992a. Circles : e/i 
channel only. Triangles : dilepton channel (e/i, ee, and [ifj,). Crosses : dilepton channel + lepton + 
jets with a b — > i-i tag. The band represents the NLO theoretical prediction for the tt cross-section 
from Ellis, 1991. The integrated luminosity was 4 pb _1 . The cross-section limits are a function of 
the top mass because the tt acceptance depends on the top mass. 



C. Searches for non-Standard Model top quark decay modes 



The lower limits on M top from pp collisions discussed in Section [VLB are not valid if the 
top quark does not decay as t — > Wb. If M iop + M b < M w , the decay W — > tb is possible 
and contributes to the width of the W boson, T(W). 

The width of the W is difficult to measure directly since the mass signature of the 
W is a Jacobian and not a Breit-Wigner (see Fig. [43]), and because the Ex resolution is 
much larger than the natural width of the W. However, constraints on the top quark mass 
independent of decay mode can be obtained from measurements of the leptonic W branching 
ratio Br(W — > Iv) = T(W —>■ lv)/T(W), based on theoretical expectations for T(W — ► Iv). 
At lowest order the branching ratio varies between Br(W — > Iv) = 1/9 for Mw < M top + 
and Br(W — > Iv) = 1/12 for for very light M top - Experimentally, a value for this branching 
ratio can be extracted from measurements of the ratio 

a(pp -> W) Br(W -> Iv) 



a(pp -> Z) Br(Z -> l+t 



(Cabibbo, 1983; Halzen and Marsula, 1983; Hikasa, 1984; Deshpande et al., 1985; Martin, 
Roberts, and Stirling, 1987; Berger et al., 1989). The value of the Z leptonic branching 
ratio is taken from the very precise measurements in e + e~ annhilations at LEP and SLC 
(Montanet et al., 1994). The values of a(pp —>■ W) and a(pp Z) can be calculated in 
QCD; note that theoretical uncertainties in these calculations largely cancel in the ratio, 
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and the uncertainty in a(pp — *> W)/a(pp — > Z) is only of order 1%. (Martin, Roberts, and 
Stirling, 1989; Hamberg, van Neerven and Matsuura, 1991; van Neerven and Zijlstra, 1992). 

Early measurements of this ratio from the SppS were found to be more consistent with 
a small top quark mass (Albajar et ai, 1987; Ansari et al, 1987). These results were used 
by several authors to set upper limits on the top quark mass, e.g. M top < 63 GeV/c 2 at 90% 
C.L. (Martin, Roberts, and Stirling, 1987); M top < 60 GeV/c 2 at 95% C.L. (Halzen, Kim and 
Willenbrock, 1988); M top < 70 ± 5 GeV/c 2 at 90% C.L. (Colas, Denegri and Stubenrauch, 
1988). However, subsequent higher statistics studies at both the SppS (Albajar et ai, 
1991a; Alitti et al, 1992b) and the Tevatron (F. Abe et al, 1990c, 1992c, 1994d; Abachi et 
al, 1995e), found results consistent with expectations for a high top quark mass. A lower 
limit M top > 65 GeV/c 2 at 95% C.L. can be extracted from the world average value of 
Br(W — > lv), (see Fig. ~ 
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FIG. 59. World average value for the inverse W —* Iv branching ratio (Abachi et al, 1995e), 
compared with Standard Model expectations as a function of the top quark mass. Also shown is 
the 95% lower limit on Mj op , independent of top quark decay mode, which can be extracted from 
the measurement. 



Direct searches for top quarks decaying into final states different from Wb have also been 
carried out at pp colliders. The decay t — > H + b is allowed in several models in which the 
Higgs sector is extended to include charged Higgs scalars {H + ). This includes non-minimal 
Standard Models, such as supersymmetry (Gunion et al., 1990). In the simplest version of 
these models, there are two Higgs doublets, and the decay t — > H + b dominates for M top < 
M^/+ Mb. At higher top masses, both t — > H + b and t — > Wb are allowed. The branching 
ratios for the two modes are functions of M top , M#+, and a theoretically unconstrained 
parameter tan/5, which is the ratio of vacuum expectation values for the two Higgs doublets. 
The charged Higgs scalar from top decay would then decay into the heaviest lepton or quark 
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pair (tv or cs), with branching fractions which also depend on the value of tan/3 (Glashow 
and Jenkins, 1987; Barger, Hewett, and Phillips, 1990; Dress and Roy, 1991). 

Searches for t — > H + b, H + — > tv based on the identification of hadronic tau decays have 
been carried out by the UAl (Albajar et al, 1991b), UA2 (Alitti et al, 1992c), and CDF (F. 
Abe et al, 1994c) collaborations. No evidence for this process was found, and limits were 
placed as a function of M top , M#+, and the H + — ► tv branching ratio. The most stringent 
limits are the result of a higher statistics analysis from the CDF collaboration (F. Abe et 
al, 1994e) based on the r — > e or fi signature, (see Fig. |60|) . 
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FIG. 60. Regions of the (M top — M^+) plane excluded at 95% C.L. for different values of the 
branching ratio for H + — > tv. From the CDF collaboration, F. Abe et al., 1994e. The vertical 
line reflects the lower limit M H + > 45 GeV/c 2 from the LEP experiments, (Buskulic et al., 1992; 
Abreu et al, 1990a; Adriani et al, 1992; Akrawi et al, 1990b). The horizontal line reflects the 
lower limit M top > 62 GeV/c 2 from F. Abe et al, 1994d, based on Br(W — > Iv). The integrated 
luminosity is 19 pb~k 
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VII. THE FERMILAB PROTON ANTI-PROTON COLLIDER 



The Fermilab pp Collider in Batavia, Illinois consists of seven accelerating structures, 
see Fig. 61]. It is the highest energy particle accelerator in the world, with a center-of-mass 
energy of 1800 GeV (Fermilab, 1984). 



A. Linear accelerators and synchrotrons 

Negatively charged hydrogen ions (H~) are initially accelerated by a Cockroft-Walton 
pre-accelerator to 750 keV. These ions are then accelerated in a 145 meter linear acceler- 
ator to 400 MeV. They are stripped of their two electrons as they are injected into the 
first proton synchrotron called the Booster. The Booster is a ring of 150 meter diameter, 
where 18 accelerating cavities are used to accelerate the protons to 8 GeV. The protons are 
then extracted and injected into the Main Ring, which is the original 2 kilometer diameter 
Fermilab proton synchrotron comissioned in 1972. 




DO detector 



FIG. 61. The Fermilab accelerator complex. From Thompson, 1994. 

The Main Ring was built with conventional electromagnets with a maximum dipole 
bending field of .65 Tesla, with 18 accelerating cavities operating at about 53 MHz, giving 
a boost of 2.5 MeV per turn. The Main Ring has accelerated protons up to 500 GeV. It is 
now routinely used to accelerate protons to 120-150 GeV which are then injected into the 
final synchrotron, the Tevatron, or into the anti-proton cooling ring. The Tevatron built in 
1983 uses superconducting magnets to bend and focus the beams, and is situated a meter 
below the Main Ring. The dipoles have a maximum dipole bending field of 4.4 Tesla and the 
focussing quadrupoles have a maximum field gradient of 67 Tesla/m. With 8 RF accelerating 
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cavities the Tevatron gives the protons a boost of 1 MeV per turn and can accelerate the 
protons up to 900 GeV with an energy resolution of ± 0.9 GeV. 

B. Anti-proton source 

As was mentioned above, the Main Ring also acts as an injector to the p source (Church 
and Marriner, 1993). The protons are extracted from the Main Ring at 120 Gev and delivered 
to a target where anti-protons are produced with peak momentum of 8 GeV. The p 's are then 
focussed in a lithium lens and fed into the Debuncher, which is a ring with a circumference of 
500 meters. The p bunches are debunched, by turning the narrow time spread of the bunch 
and large momentum spread into a large time spread and a narrow momentum spread. After 
this bunch rotation, the p's are stochastically cooled (van der Meer, 1972), further reducing 
the phase space of the beam. The beam is then transferred to the Accumulator, added to the 
anti-protons already stored there using a process called stochastic stacking. When enough 
p's have been accumulated to make an intense beam, the beam is extracted and injected 
back into the Main Ring, where it is then accelerated and injected into the Tevatron. 

During the period when the top data was collected, the accelerator operated with 6 
bunches of protons and 6 bunches of anti-protons circulating in the machine. It takes 
approximately 2.5 hours of shining protons on the p source target every 2.4 seconds in order 
to accumulate enough p's to ensure a reasonable luminosity for pp collision. To optimize 
collider livetime, this process of creating and stacking the anti-protons took place during 
collider operation. This meant that there was often beam in the Main Ring while the collider 
experiments (CDF and DO) were taking data. The halo from this beam interacted with the 
walls of the accelerator and sprayed particles into the two collider detectors, depositing large 
amounts of energy in the calorimeters. In the case of CDF, the Main Ring was diverted via 
a dogleg up above the detector and shielded by a steel structure approximately one meter 
thick. Nevertheless, occasionally extra energy was deposited in the calorimeters. These 
events were rejected offline. The DO experiment had no such dogleg or beam diversion and 
consequently the Main Ring ran straight through the DO detector, approximately 2.5 meters 
above the Tevatron beampipe. The DO collaboration was forced to turn off data acquisition 
while the Main Ring beam was passing through their detector, resulting in a 15% livetime 
loss. 

C. Collider 

The p and p bunches are approximately 50 cm in length due to the RF frequency of 
the accelerator. This bunch length determines the long luminous region at the interaction 
points, which is roughly Gaussian and has a sigma of approximately 30 cm. The length of the 
luminous region has many limiting features in terms of triggering and solid angle coverage 
by the detectors. For instance, as will be discussed later in this review, the silicon vertex 
detector in the CDF experiment was 50 cm long and thus, due to the size of the luminous 
region, had a limited acceptance of roughly 60% for bottom quarks from top quark decays. 
Furthermore, the measurements of transverse energies at the trigger level was smeared by 
the lack of knowledge of the event-by-event interaction position. 
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Accelerator radius 


1000 m 


Maximum beam energy 


900 GeV 


Injection energy 


150 GeV 


Peak luminosity 


2 x 10 31 cm" 2 s" 1 


Number of bunches 


6p, 6p 


Intensity per bunch 


ps 10 n p, 5 x 10 10 p 


Crossing angle 


0° 


Bunch length (1 a) 


50 cm 


Transverse beam radius (1 <r) 


~ 25 /im 


Energy spread 


0.15 x 10" 3 GeV 


RF frequency 


53 MHz 


p stacking rate 


«3.5 x 10 10 /hour 


Beam crossing frequency 


290 kHz 


Period between crossings 


3.5 /is 



TABLE VII. Parameters of the Fermilab Tevatron Collider. 



The Tevatron Collider complex was first commissioned with a short run in 1985. The first 
high-luminosity run took place in 1988-89, at which time only the CDF detector was ready 
for data-taking. The peak luminosity for that run was approximately 2 x 10 30 cm~ 2 sec -1 , a 
factor of two higher than the initial design. The second set of high luminosity collider runs, 
Run la and Run lb, begun in 1992 with data taken by both the DO and CDF detectors. 
During these runs, the peak luminosity was 2 x 10 31 cm~ 2 sec~ 1 , while the average luminosity 
was about 1.4 x 10 31 cm _2 sec _1 . The lifetime of the beams in the Collider and the p stacking 
rate was such that new protons and anti-protons were injected into the Collider once a day. 
Table VII lists the relevant parameters of the Collider. 



VIII. DISCOVERY OF THE TOP QUARK 

In this Section we review the recent results that finally led to the establishment of the 
existence of the top quark. The discovery of the top quark was made possible by the 



remarkable success of the Tevatron Collider project, see Section [VII] . These results come 
from the two collider experiments (CDF and DO) at Fermilab's Tevatron. They are based 
on data from the 1992-93 (Run la) and 1994-1995 (Run lb) runs of the collider. Analysis of 
data from the remainder of Run lb was in progress at the time that this review was being 
written. Preliminary results from both CDF and DO are fully consistent with those from 
the earlier data sets. 

The data sets that were used to discover the top quark were collected during Run la 
and the first half of Run lb. The total integrated luminosities were 67 pb _1 for CDF, and 
between 44 and 56 pb _1 , depending on top decay mode, for DO. The difference in integrated 
luminosities between the two experiments is due mostly to the fact that the Main Ring 
accelerator at Fermilab which operates asynchronously from the Tevatron, runs through 



74 



the DO calorimeter, see Section |VII| . Data taking in the DO experiment must be disabled 
whenever a Main Ring proton bunch crosses the detector. 

Given the ti production cross-section at the Tevatron (see Section [TVD , the number of ti 
pairs produced in an exposure of 67 pb -1 is expected to be between approximately 6800 (for 
M top = 100 GeV/c 2 ) and 150 (for M top = 200 GeV/c 2 ). This is a small number of events, 
especially in light of the fact that not all the the possible ti decay channels can be exploited, 
see the discussion in Section M. Thus, the high luminosity delivered by the Tevatron was 
crucial in enabling the experimenters to isolate the very rare top signal. 

Both the CDF and DO top searches assumed Standard Model decay (t — > Wb) of the top 
quark, and were based on the ti dilepton and lepton + jets signatures discussed in Section [V]. 
In all cases, the number of observed events, after a selection procedure designed to maximize 
the acceptance to top quarks, was compared with the number of expected events from non-tt 
sources. An excess of events over the background prediction then constitutes evidence for 
the top quark. 

At the beginning of Run la, the best lower limit on the mass of the top quark was 91 
GeV/c 2 , from the 1989 CDF top search, see Section [VT|. The initial top selection criteria 
for both CDF and DO were optimized for detection of a low mass (M top ~ 120 GeV/c 2 ) top 
quark. A lower limit of M top > 131 GeV/c 2 at the 95 % C.L. was established by DO from 
the Run la data with an integrated luminosity of 13.5 pb -1 (Abachi et ai, 1994). 

Initial evidence for the top quark was reported by the CDF collaboration (F. Abe et al, 
1994a) based on analysis of the Run la data set only, which had an integrated luminosity 
of 19.3 pb -1 . The excess in the number of top candidate events over the background pre- 
diction was 2.8 standard deviations. Additional kinematic features of the data, such as a 
reconstructed mass peak, also supported the ti hypothesis. However, these results were not 
deemed sufficient to unambiguously establish the existence of the top quark. Following the 
publication of a lower limit on M iop , the DO requirements were reoptimized, i.e. tightened, 
for higher top masses. With the optimized requirements, a statistically not very significant 
excess of events (1.9 standard deviations) was also found in the Run la DO data (Abachi 
et ai, 1995a and 1995d). With the addition of the data from the first half of Run lb, a 
statistically convincing excess of events emerged from the analyses of the data sets from 
both collaborations (F. Abe et ai, 1995a; Abachi et ai, 1995b). 

The CDF and DO searches for ti in the dilepton channel were similar and are discussed 
in Section |VIII.A| . On the other hand, the search strategies in the lepton + jets channel 
were quite different. At the beginning of Run la, a silicon vertex detector was installed at 
CDF (Amidei et ai, 1994). Given the excellent 6-tagging capabilities of this device (see 
Section V.E.4), 6-tagging was used to separate the top signal in the lepton + jets channel 
from the W+ jets background. Both vertex tagging and lepton tagging (b —>■ e as well as 
b —> fi) were used in CDF. In contrast, a large fraction of the top sensitivity in this channel 
for the DO experiment came from kinematic separation of ti and W+ jets, although lepton 
tagging (b — > fi only) was also employed. Analyses of the kinematic properties of lepton + 
jets events were also performed on the CDF data (F. Abe et al, 1994a, 1995b, and 1995c). 
Excesses of top-like events were also seen in these CDF studies, and were used to confirm 
the results of the 6-tagging observations. 

In the remainder of this Section, we will describe in some detail the results from the 
CDF and DO top searches. Dilepton searches will be reviewed in Section |V111.A], followed 
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by a summary of the searches in the lepton + jets + 6-tag channel in Section [VIII. B| , and 
of kinematic separation of ti and W+ jets in Section VIII. Q . The results from the two 
experiments will then be summarized in Section VIII. D, and their measurements of the 
pp — > ti cross-section will be presented in Section |VIII.E . 

Before beginning the discussion of the CDF and DO results, we wish to point out a 
difference in when the measured jet energies are corrected back to the parent parton momenta 
in CDF and DO. In all DO analyses, correction to the measured jet energies are made 
prior to the selection of the sample. (The jet energy correction procedure will be discussed 
extensively in Section [LX..A.3| ). For most CDF analyses, the measured jet energies are not 
corrected at the event sample selection stage, or in the following analysis. Corrections are 
only applied when the measurement of the quark or gluon energy is needed, for example in 
the measurement of the top mass, see Section |D£|. Unfortunately, this makes it somewhat 
difficult to compare results from the two experiments. For the jet energies relevant to the 
top search (15 GeV < E T < 150 GeV) the multiplicative correction factor for a CDF jet 
energy cluster is of order 1.5. 



A. Dilepton analysis results from CDF and DO 

The searches for top quarks performed in the dilepton channel by DO and CDF are 
similar. We have discussed the dilepton analysis approach at a general level in Section [V.D . 
What follows are the results of these two analyses. The main difference in analysis strategy 
between the two experiments is due to the fact that CDF has momentum determination from 
charged particle tracking in the central region whereas DO does not, and that DO imposes 
more stringent requirements on the transverse energies of the jets. 





Leptons 


Jets 






Channel 


E T (e) 


PtGu) 










en + jets 


> 15 GeV 


> 12 GeV/c 


> 2 


> 15 GeV 


> 20 GeV 


> 120 GeV 


ee + jets 


> 20 GeV 




> 2 


> 15 GeV 


> 25 GeV 


> 120 GeV 


HfM + jets 




> 15 GeV/c 


> 2 


> 15 GeV 




> 100 GeV 



TABLE VIII. A partial summary of selection criteria for ti — > dileptons in DO. Jets must have 
\rj\ < 2.5 and their energies are corrected. The lepton pseudorapidity coverages are \rj\ < 2.5 for 
electrons, \rj\ < 1.7 for muons in Run la, and |ry| < 1.0 for muons in Run lb. The loss of muon 
coverage is due to ageing of the forward muon chambers. The requirement is discussed in the 
text. 



Both analyses require two high leptons, at least two jets and in most cases, high FJ T . 
The exceptional case is the DO search for decays in the channel, where the FJ T requirement 
is removed. The resolution of the DO muon momentum measurement is limited to of order 
20% by multiple scattering in the toroids (see Table [TV]) , and therefore in the \x\x channel the 
measurement of the FJ T is poor and is not used in the DO analysis. The major backgrounds 
to the top quark signature in the dilepton channel are due to events in which jets fake 
leptons, Z — > ee/ hh/tt, Drell-Yan lepton pair production, WW and bb production. These 
backgrounds are reduced dramatically by the requirements listed in Table VIII and Table [TX. 
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Channel 


Leptons 


Jets 


E T 


A0(# T , lepton or jet)> 20° 


E T (e) 


Pt(m) 


Njet 




en + jets 


> 20 GeV 


> 20 GeV/c 


> 2 


> 10 GeV 


> 25 GeV 


-> t T > 50 


ee + jets 


> 20 GeV 




> 2 


> 10 GeV 


> 25 GeV 


-> # T > 50 


\x\i + jets 




> 20 GeV/c 


> 2 


> 10 GeV 


> 25 GeV 


-> # T > 50 



TABLE IX. A partial summary of selection criteria for tt — > dileptons in CDF. Jets must have 
|^| < 2.5 and their energies are not corrected. The lepton pseudorapidity coverage is \rj\ < 1.0. The 
Er requirement is tightened when the missing transverse energy vector is collinear with either a 
lepton or a jet. See the discussion in the text. 



In particular, as discussed in Section |V.D| , asking for the presence of two jets in the event is 
a very powerful discriminator between signal and all sources of background. 

There are two other differences in the kinematic requirements applied by CDF and DO. 
In the DO analysis, there is a requirement on H^, defined as = T,j ets (Ex) + which is 
the scalar sum of all jet transverse energies plus the highest transverse energy electron in the 
event, if there is one. This is a kinematic variable that attempts to discriminate between the 
top signal and background, by exploiting the difference between the total transverse energy 
of the jets in top and background events. In Fig. |62| we show Monte Carlo predictions of HV 
for the principal backgrounds and for ti events with M top = 200 GeV/c 2 . 

For the CDF analysis, in events where the % T is nearly collinear with the Pt of the 
leptons or jets, the JZi T requirement is tightened. As we will show below, this procedure is 
somewhat effective at suppressing a number of background sources. We will now discuss the 
dilepton CDF and DO results, starting with the background sources, and the specific CDF 
and DO selection requirements designed to discriminate against them. 



1. Z — * ee/^Lfi background 

Z decays constitute an important background to the ti signal in the dilepton channel. 
Following the E^(e), Pt(aO and N jet requirements listed in Table VIII and Table |pX| there is 
still a substantial background from Z + > 2 jets, in which the Z decays to e + e~ or fi + fi~. 
Since Z decays to charged leptons do not produce any neutrinos, no JZ! T is expected in these 
events. Therefore, requiring large $ T eliminates these backgrounds completely, except for 
the effects of mismeasurements of the missing transverse energy in the detector. 

To eliminate the Z 11 background, e + e~ and events with dilepton invariant mass 

between 75 GeV/c 2 and 105 GeV/c 2 are rejected in the CDF top search. This requirement is 
expected to be approximately 75% efficient, roughly independent of top mass, for dileptons 
(e + e~ and fi + ^~) from tt decays. In the DO analysis, the requirements used to discriminate 
against the Z background differ in the e + e~ and n + fi~ channels because the resolution on the 
lepton energy/momentum measurement is so much better for electrons than it is for muons. 
In the e + e~ channel, if the invariant mass (M ee ) of the pair is consistent with the Z mass (79 
GeV/c 2 < M ee < 103 GeV/c 2 ), the missing transverse energy requirement is tightened to 
% T > 40 GeV. Note that this is in contrast to the CDF approach of rejecting any event that 
may be consistent with Z decays, regardless of the size of the missing transverse energy in 
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10 




H T (GeV) 

FIG. 62. DO Monte Carlo predictions of the variable Hr for dilepton decays of a 200 GeV/c 2 top 
quark (solid) and principal dilepton backgrounds (dashed). From Abachi et al, 1995b. 

the event. As was discussed earlier, no Hl T information is used by DO in the /i + /U~ channel. 
Z decays into muons are then rejected on the basis of an overall kinematic likelihood fit to 
the Z — > hypothesis. 

Since the Monte Carlo simulation is unlikely to correctly model the tails of the Hl T reso- 
lution, the DO collaboration relies on multijet data to estimate the remaining Z — > ee+ jets 
background after the JZl T requirement. As a result of JZl T studies in the multijet sample, the 
fraction of Z — ► ee events with % T > 40 GeV is estimated to be less than 2 x 10~ 4 (Abachi et 
ai, 1995d). The CDF requirement is more conservative and loses acceptance in return for 
certainty that there will be no Z decays in the top candidate sample. In Fig. |53]we show the 
distribution of the invariant mass of the electron-positron pair versus the J^ T for DO events. 
Expectations for tt are displayed in Fig. 

It is difficult to compare directly the Et resolutions of CDF and DO. The parametrizations 
of the Et resolutions in the two experiments do not have the same functional form. The CDF 
resolution grows linearly with ^Y^E? whereas the DO resolution appears to grow linearly 
with YIEt, see Table |V[ The tails of the distribution in CDF are dominated by cracks in the 
calorimeter while the DO detector design minimized cracks. The radioactive noise from the 
uranium plates in the DO calorimeter adds of order 25 GeV per event to the total energy 
making difficult a direct comparison of the resolution at a particular T,Ex- 

2. Z — ► tt background 

Z — > tt, followed by leptonic decays of both taus, is a very important background 
source. These events cannot be eliminated by an invariant mass cut, due to the presence of 
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FIG. 63. Dielectron mass versus missing transverse energy in ee+ 2 jet events from DO (Grannis, 
1995). The integrated luminosity is 55.7 pb _1 . Also shown is the DO $ T requirement in the 
dielectron channel. 



unmeasured neutrinos from r — > Ivv. Furthermore, these neutrinos can give rise tO Er 111 
the event. The requirements on jet multiplicity, high or leptons, and high Hl T are all 
effective at reducing the background, see the discussion in Section [V.D| , and also Fig. |38. 



In Fig. ^5| we show the J?! T versus dilepton mass for leptons from Z — > rr decays from the 
ISAJET Monte Carlo. This clearly indicates that the background is not eradicated by the 
Ex requirement alone. 

In the decay of a high momentum r, the decay products are highly collimated, due to 
the small mass of the r lepton. As a result, in these events the Ex will often point in the 
direction of one of the two leptons in the transverse plane. Therefore, in the CDF analysis 
the Ex requirement is tightened to Er > 50 GeV for events in which the azimuthal angle 
) of the Er vector is within 20° of the azimuthal angle (4>i ep ton) of either of the two 
leptons, (see Fig. |66|) . 

There is no such A</> requirement for DO, see Table VIII. However the requirement made 
on H^, see Table VIII, has a rejection power of approximately 2.5 for Z — > rr+ two or more 
jets. 
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FIG. 64. Dielectron mass versus missing transverse energy in tt — > ee+ 2 jets events from the DO 
Monte Carlo (Grannis, 1995). Also shown is the DO $ T requirement in the dielectron channel. 
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FIG. 65. Missing transverse energy versus dilepton mass for Z — > tt + 2 jets Monte Carlo events. 
From the DO collaboration, Grannis, 1995. Also shown is the DO $ T requirement in the dielectron 
channel. 



80 



180 






^ 160 


_ a 


• 


2 140 
























X 100 






1 80 








• 

- • • 


• 

• 


W 60 


— • 

— • 




40 


• 


• 

• • 

• 


< 20 


• 

• • 


• 

• • 
• • 


• 

• 

• 

• 





" .1 . 1 


• • 

■ \ • 


1 

•j 1 1 1 1 1 1 1 1 







50 100 

E T (GeV) 



150 



FIG. 66. CDF Monte Carlo distribution of the azimuthal angle between the $ T and the closest 
lepton as a function of $ T in Z — > tt events with two jets. Also shown is the boundary of the 
CDF $ T requirement. From the CDF collaboration, F. Abe et al, 1994a. 
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3. Drell-Yan Background 



Continuum Drell-Yan production of lepton pairs accompanied by two jets is a background 
not addressed by the Z mass requirement. In the CDF analysis, the A0 requirement between 
the F>r and the closest jet (see Table |IX|) is designed to reduce this background. There is no 
intrinsic FJ T in these events, however jet energy can be lost in cracks, and the measurement of 
the jet energy can fluctuate, faking a FJ T signal. In these cases, the FJ T tends to point towards 
one of the two jets in the event. In order to reduce this background, the % T requirement in 
CDF is raised to 50 GeV when the azimuthal angle between the FJ T direction and a jet is 
less than 20°, see Fig. |67. 

DO's approach is different. First of all, the Ht requirement is more effective than a 
minimum number of jets requirement because the jets in Drell-Yan + jets events are in 
general softer than the jets in ti events, for sufficiently high top mass. Furthermore, to 
reduce the Drell-Yan background, the % T requirement is raised to 25 GeV for dielectron 
events. Recall that in the fifi channel, no cut is made on FJ T because the muon energy 
measurement is poor. 




E T (GeV) 

FIG. 67. Distribution of the azimuthal angle between the JZi T and the closest jet as a function 
F^t for Z — > 11 events with two jets. The $ T characteristics of these events should closely resemble 
those of Drell-Yan events. This is CDF data with an integrated luminosity of 19 pb _1 . Also shown 
is the boundary of the CDF $ T requirement. From the CDF collaboration, F. Abe et al, 1994a. 



4. WW background 

The expected background from diboson production, WW, WZ, and ZZ, is small but not 
at all negligible. The WW process is the dominant diboson background to the ti dilepton 
signature. The WZ and ZZ production cross-sections are expected to be lower and are 
further suppressed by the small Z —> 11 branching ratio, see the discussion in Section |V.D| . 
Neither collaboration has published evidence for the process pp —>■ W + W~, however the 
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cross-section expected from theory is slightly higher, of order 10 pb, than the equivalent ti 
cross-section, which is expected to be a{ti) « 5 pb for M top = 175 GeV/c 2 , (see Fig. [HI). 

WW events include large transverse momentum leptons and % T . In Fig. |68| we show the 
expected distributions of FJ T and A0 for ti and WW events. The kinematic properties of 
the two processes are very similar. The most effective way to reduce these backgrounds is to 
require the presence of high transverse momentum jets. The requirement imposed by DO 
is more effective in discriminating against WW than the simple minimum-number-of-jets 
requirement employed by CDF. The remaining background is estimated entirely from Monte 
Carlo, using the theoretical expectations for the WW cross-section. 




FIG. 68. CDF Monte Carlo expectations for A(f> versus F T for (a) tt with M top = 180 Gev/c 2 , and 
(b) WW + > + 2 jet production. The vertical axis is the smallest of the azimuthal angles between 
the Ft vector and the leptons or jets. The symbol MET in the vertical-axis label refers to the F T . 



5. bb and fake lepton backgrounds 

Doubly semileptonic decays in bb events and events with fake leptons constitute another 
important background to the dilepton signature. These backgrounds are considerably re- 
duced by requiring two isolated high transverse momentum leptons, and high FJ T . The bb 
background is calculated from Monte Carlo, normalized to the rate of lower momentum 
bb — > dileptons events. The probability of finding two high momentum isolated fake leptons 
in an event is exceedingly small. The main fake dilepton background is due to W+ jets 
events, where one of the leptons is from the W — > lu decay, and one of the jets is misidenti- 
fied as a lepton. The background is then estimated by multiplying the number of observed 
W+ jets events by the probability for a jet to fake the isolated electron or muon signature 
as determined from samples of jet events. This probability is typically of order 10 -4 . 



S3 





3 


c 






2.5 






CD 


2 






o 




c 


1 .5 


o 




-) — ' 




Q- 


1 


0) 


1— * 
1 i 1 


0.5 


(Ml 





9- 




< 











ENTRIES 233 








All jet multiplicities 


- • 

— • 
_ • 




• 


ij 


lit* * 

-' 


. /2/^7+1jet 
• 2j 




i i I • r 


•0j 


i i i i i i i i i i 



20 



40 



60 80 100 120 140 
Missing E, (GeV) 



o 



CD 

s_ 
O 

c 
q 

CD 



3 
2.5 

2 
1 .5 

1 

0.5 



N jets * 2 



Signal Region 



/^ity + 1 jet 

jjLjj. + 2] candidate 





9- 



20 40 60 80 100 120 140 

Missing E, (GeV) 



FIG. 69. CDF dilepton results in the \i\x channel. The integrated luminosity is 67 pb _1 . From 
Roser, 1995. The vertical axis is the smallest of the azimuthal angles between the Ex vector and 
the leptons or jets. The symbol MET in the vertical-axis label refers to the % T . 
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FIG. 70. CDF dilepton results in the e/x channel. The integrated luminosity is 67 pb _1 . From 
Roser, 1995. The vertical axis is the smallest of the azimuthal angles between the $ T vector and 
the leptons or jets. The symbol MET in the vertical-axis label refers to the Ft- 
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(TDF 




Background Source 






Z — > TT 


0.38 ± 0.07 


0.24 ± 0.07 


WW 


0.21 ± 0.07 


0.06 ± 0.03 


Fakes, bb 


0.26 ± 0.16 


0.12 ± 0.04 


Z — ► ee or /x/x 




0.24 ± 0.03 


Drell Yan 


0.44 ± 0.28 




Total background 


1.3 ± 0.3 


0.65 ± 0.15 


it expectation, M top = 180 GeV/c 2 


2.4 


1.2 


Data 


5 e/x 1/x/x 


2 e/x 1 /t/x 



TABLE X. Background sources, ti expectations, and event yields in the CDF (Roser, 1995) and 
DO (Narain, 1995) dilepton analyses. One additional CDF candidate event, consistent with the 
Z — ► fjbfxj hypothesis, is not included in this table, see the discussion in the text. The integrated 
luminosity for the CDF data is 67 pb _1 . For the DO data, the integrated luminosity is 55.7 pb _1 for 
ee, 44.2 pb _1 for /x/x, and 47.9 pb _1 for e/i. The ti expectations are obtained using the calculation 
of a(tt) from Laenen, Smith, and van Neerven, 1994. 

6. Results 

The acceptance for the CDF analysis for M top = 180 GeV is e = (0.87 ± 0.10)%, where 
this acceptance includes the branching ratio for the dilepton mode (4/81, see Table PTT] ). 
The equivalent acceptance for DO is (0.55 ± 0.04)%. With these acceptances, and given 
the integrated luminosities and the theoretical expectations for a(tt), both experiments are 
sensitive to ti production for top masses up to M top ps 200 GeV/c 2 . The results of the CDF 
and DO ti searches in the dilepton channel are summarized in Table [X|. 

CDF finds 7 dilepton candidates in 67 pb _1 with 1.3 ± 0.3 expected background events. 
Of these 7 candidate events, five are e/x+ jets events and two are /x/x+ jets events. One of 
the two /x/x events looks very much like a radiative Z decay, Z — > /i/17, since the invariant 
mass of the /x/17 is 86 ± 7 GeV/c 2 . Although the background from radiative Z decays was 
estimated to be less than .04 events, CDF conservatively a-posteriori removes this event 
from the top candidate sample and is left with 6 events (5 e/x and one /x/x). Acceptances for 
top are such that 60% of ti are expected in e/x and 40% in ee//x/x. Figs. and |70| show 
the CDF data in the A0 — % T plane. As can be seen from these figures, the e/x channel is 
much cleaner due to the absence of the Drell- Yan contribution. 

DO finds 3 candidates in approximately 50 pb _1 with .65 ± 0.15 expected background 
events. Of these three events, two are e/x+ jets and one is /x/x+jets. This is consistent with 
the ratio of expected number of events, e/x : ee : /x/x = .34 : .25 : .11 for a 200 GeV top 
quark. 

Thus, both experiments see an excess of dilepton events over the total background pre- 
diction. The size of this excess is loosely consistent with expectations for ti. It is interesting 
to notice that there are some significant differences in the CDF and DO background sources, 
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FIG. 71. Distributions of for dilepton data for DO (Grannis, 1995) The integrated luminosity 
is ~ 50 pb _1 . The solid line is the expected Hr distribution for top events in the dilepton channel. 
The dashed line is the expected backgrounds and the solid histogram is the data. See text for 
details. 



see Table [X[ The highest background source in the CDF analysis is due to Drell-Yan events. 
This background is negligible in DO, probably due to the DO requirements, the stricter 
Et requirement in the ee channel, and the hermeticity of the DO detector. Also, Z — > ee 
and fifi are totally eliminated by the CDF invariant mass requirement, while they make up 
a substantial fraction of the total background in DO. 

As partial evidence that the excess of events is not due to a background fluctuation or to 
underestimation of the backgrounds, we present in Fig. [71] the fir distribution of dilepton 
events from DO. As was discussed earlier, in DO IdV is defined as the scalar sum of the 
transverse energy of all of the jets in the event + the Ey of the highest Et electron (for the 
ee and e/i channels). As we argued throughout this Section, and as shown in Fig. |62|, 
is a powerful discriminant between signal and background. The DO candidate events are 
distributed in HV m a manner more consistent with ti than background. 

In addition, in the six CDF dilepton candidate events, there are five jets (in three events) 
that are tagged as fe-jets by the algorithms which will be described in Section |V111.B[ This 
is to be compared with the expectation of approximately 0.5 6-tagged jets if all six events 
were background, and approximately 3.6 6-tagged jets if all six events were ti. This suggests 
that the excess of dilepton + jets events over the background prediction is correlated with 
the presence of 6-quarks in the event, as expected if the excess were due to ti. 
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B. Lepton + jets + 6-tag 



As was discussed in Section |V.E.4| , one of the most effective ways of isolating a top signal 
in the lepton + jets mode (ti — > qq Iv bb) is to tag the 6-quarks in ti events. This can be 
achieved by searching for leptons from the semileptonic decay of 6-quarks (lepton tagging), 
or by exploiting the long lifetime of 6-hadrons (vertex tagging). Lepton tagging has been 
used by both CDF and DO to extract a top signal; vertex tagging, which requires momentum 
analysis, as well as precise vertex detection capabilities, has been performed by CDF only. 
In Section |VIII.B.l| we will review the CDF and DO selections of lepton + jets data, before 
the 6-tag requirement. The results of the 6-tag analyses on these data are then discussed in 
Sections |VIII.B.2| and [VIII. B. 3] for lepton and vertex tagging respectively. A brief summary 
of the tagging searches, as well as results of cross-checks performed on samples of Z+ jets, 
are presented in Section |VIII.B.4. 



1. Selection of lepton + jets data before 6-tagging 

The CDF and DO selections of the isolated e or fi + jets data for the 6-tagging analyses 
are summarized in Table [XT. The event selection employed by the two collaborations are 



similar. In principle, one would expect to detect four jets in a tt — > lepton + jets event; 



however, as discussed in Section |V.B.2| , the detected number of jets can be smaller (see for 
example Fig. |35|). This is because jets from top decay sometime have Et below threshold, 
or two or more of these jets can be close enough to each other that they are reconstructed 
as a single jet. Therefore, in order to maintain high efficiency, the minimum number of 
jets required by both CDF and DO is three rather than four. The one significant difference 
in the selection requirements between the two experiments is that DO imposes a further 
requirement on the minimum scalar sum of the transverse energy of all jets in the event 
(H<r, see Table Kip. This requirement is expected to improve the rejection against the 



dominant W + jets background, especially for high top quark masses (see Fig. [16]). No such 
requirement is imposed on the CDF data in order to maintain high ti detection efficiency for 
M top as low as 100 GeV/c 2 , and because of the superior background rejection capabilities of 
the CDF vertex tagging analysis. 

The selections summarized in Table |XJ yield data samples consisting mostly of W+ jets 
events, as well as a contamination from QCD events, and hopefully also a ti component. 
The QCD contamination is due to events with fake leptons, as well as semileptonic decays 
of 6-quarks in pp — > bb events. In almost all of the channels, these events are estimated to 
contribute approximately 10% to the event sample; the exception is the DO fi+ jets channel, 
where this background is a factor of 2-2.5 larger, due to the poor resolution of the DO muon 
momentum measurement. These data samples also include small contributions from Z and 
diboson events. 

The numbers of expected ti events in the data samples (see Table |X| ) depend on M top , 
since both the ti production cross-section and acceptance depend on M top . Given the theo- 
retical expectations for the ti production cross-section, the expected signal-to-background in 
these data samples before applying 6-tagging (pre-tag samples) varies between approximately 
1/16 and 1/3, depending on top mass, and selection details. 



CDF DO 

Lepton E T or P r > 20 GeV > 20 Gev (e); > 15 GeV (//) 

Lepton rapidity coverage \r]\ < l (a) \v\ < 2 (e); \r]\ < 1.7 (n)^ 

t T > 20 GeV > 20 GeV 

Number of jets > 3 > 3 

Jet Et > 15 GeV (uncorrected) > 20 GeV (corrected) 

Jet rapidity coverage \r]\ < 2 \rj\ < 2 

Jet cone clustering radius 0.4 0.5 

H T - > 140 GeV 

Integrated luminosity 67 pb _1 48 pb _1 (e); 44 pb _1 (fi) 

No. of events in data sample 203 66 

No. of expected ti events, M top = 140 GeV/c 2 ps 95 « 22 

No. of expected ti events, M top = 200 GeV/c 2 w 13 « 6 



TABLE XI. CDF and DO lepton + jets requirements for the 6-tag analysis. The variable H-j- in DO 
is defined as the scalar sum of the Et of all the jets. The expected number of top events are derived 
from the theoretical estimate of the ti cross-section (Laenen, Smith, and van Neerven, 1994) and 
the calculated acceptances. Notes : ( a ) the CDF muon chambers only cover ~ 2/3 of the solid 
angle for 0.6 < \rf\ < 1; ^ Muons in DO are restricted to \i]\ < 1 for the last ~ 70% of data, due to 
ageing of the forward muon chambers. 



The dominant physics background in the lepton + jet + 6-tag channel is due to the 
associated production of a If boson and a pair of heavy quarks (WQQ, Q = b or c, see 
Section |V.E.4| , Fig. ^2]). As discussed in Section |V.E.4| , the fraction of W+ > 3 jets events 
containing a QQ pair is expected to be of order 3% for Wbb and 5% for Wcc. Therefore, 
given the signal-to-background levels in the lepton + > 3 jets samples before demanding a 
6-tag, a ti signal is expected to stand out after application of a 6-tag requirement, provided 
that instrumental backgrounds can be kept under control. 



2. Lepton tagging in CDF and DO 

Leptons from 6-quarks in ti arise from direct (6 — > civ) or cascade (6 — > c — > slv) decays. 
Given the 6- and c-quark semileptonic branching ratios, and including the contribution from 
semileptonic decays of c-quarks from W — > cs, there is on average approximately one lepton 
(e or fi) from 6 or c decay in each ti — > lepton + jets event in addition to the lepton from W 
decay. These additional leptons tend to have low transverse momentum, see Fig. and to 
be non-isolated because of the presence of nearby hadrons from the 6-quark fragmentation 
and the 6-hadron decay, (see Fig. |72|). Detection of these leptons is therefore more difficult 
than detection of the isolated high leptons from W decays. 

There are two main differences in the lepton-tagging algorithms developed by the two 
collaborations. The first difference is that the DO analysis only tags muons, whereas CDF 
tags both muons and electrons. This is because detection of non-isolated electrons in DO is 
not as effective as in CDF, mostly due to the absence of a magnet for momentum measure- 
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O 20 40 60 80 1 OO 1 20 140 
Scalar sum of P T in a cone of 0.4 (GeV/c) 

FIG. 72. The expected scalar sum of the transverse momenta of all particles in the ti final state 
within a cone of radius 0.4 of the lepton from b — > civ or b — > c — > s/za From the ISAJET Monte 
Carlo, for Mi 0p = 175 GeV/c 2 . The momenta of the lepton and the neutrino are not included in 
the sum. 



ment. In CDF electrons from 6-decays can be identified with sufficient background rejection. 
However, the efficiency for detecting these electrons is lower than the efficiency for detecting 
muons, and only of order 30% of the lepton tagging efficiency in CDF is due to electron tags. 
The second difference between the two experiments resides in the minimum Pt requirement 
for these tagging leptons. In DO, this requirement is set at 4 GeV/c, which corresponds to 
the minimum Pt for a muon to penetrate the DO muon detectors; in CDF this minimum 
Pt requirement is set at 2 GeV/c. The lower Pt requirement improves the efficiency for 
detecting muons from cascade decays, especially for low top quark masses, (see Fig. [49|) . 
however it also results in higher background levels. 

Additional requirements are imposed in the DO analysis to reject events with tagging 
muons collinear or back-to-back with the direction of the E! T vector. These requirements 
are designed to reject QCD events. The rapidity coverages for tagging leptons in the two 
experiments are the same as those for the high Pt leptons from W decays, see Table XI . 
For M top > 140 GeV/c 2 , the ti lepton-tagging efficiency is 20% for both experiments. The 
greater muon coverage of the DO detector makes up for the higher Pt threshold and the 
absence of electron-tags. 

Backgrounds to 6-tagging were discussed extensively in Section [V.E.4| . The WQQ back- 
grounds in both CDF and DO are calculated by assuming that the heavy flavor content of 
jets in W + jets events is the same as that of generic jets. This corresponds to the Method 
I background estimate of Section |V.E.4| , and is expected to yield an over-estimate of the 
background. In DO, a tagging rate per jet as a function of jet Et is defined as the ratio 
of the number of tagged jets divided by the total number of jets as a function of Et, (see 
Fig. |73|). The background to the lepton tagged signal is then calculated by convoluting 
this tagging rate with the Et spectrum of jets in the W+ > 3 jets sample. As discussed 
in Section [V.E.4| , this procedure yields simultaneously an estimate of the instrumental and 
WQQ backgrounds. 

A similar procedure is used in CDF, with the difference that track-tag rates instead of 
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jet-tag rates are used. The track-tag rate is defined as the ratio of the number of tracks in 
jet events that are tagged as leptons divided by the total number of tracks. To calculate 
the background, electron and muon track-tag rates as a function of P-r are convoluted with 
the Pt spectrum of tracks in the W+ > 3 jets sample. The CDF muon track-tag rate is 
displayed in Fig. [73]; the analogous electron track-tag rate is considerably smaller, see F. 
Abe et a/., 1994a. 
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FIG. 73. Tagging rate for jets as a function of in the DO detector as measured from a sample 
of fake electron + jets events. From Snyder, 1995a. 

An additional physical background to the tagged W+ jets search is due to production of 
a W and a single c-quark, W+ charm, (see Fig. |5T|) . In the DO analysis it is assumed that 
this background is automatically included when the jet-tag rate is convoluted with the jet 
Et spectrum. On the other hand, this background is calculated explicitly by CDF and is 
added in separately. 

A background contribution from Z — > fifi events is also present in DO. The two muons 
can result in misinterpreting such an event as W —>■ \w event with a fi tag. The poor muon 
momentum resolution does not allow for a clean removal of these events via an invariant 
mass cut, as is done in CDF. Just as in the DO dilepton analysis (see Section |VIII.A| ), these 
events are identified using a global event x 2_ test for consistency with the Z hypothesis. 
Since the Z removal procedure is not 100% efficient, the DO background estimate in the 
lepton + jets + 6-tag search also includes a contribution from residual Z — > /i/i events in 
the sample. Additional small backgrounds are due to all sources of dileptons discussed in 
Section [VIII. A| (e.g. Z — > tt, dibosons, bb). These backgrounds are expected to be small. 
They are included in the CDF background estimate, and taken as negligible by DO. 

Both CDF and DO see an excess of tags in the W+ > 3 jets samples over their respective 
background estimates. The DO collaboration finds 6 events on an expected background of 
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FIG. 74. Muon track tagging rate as a function of in the CDF detector as measured from a 
sample of jet events. From F. Abe et al, 1994a. 



1.2 ± 0.2; the CDF collaboration finds 23 tags in 22 events, with an expected background 
of 15.4 ± 2.0 tags. The most powerful check of the background calculation procedure is to 
repeat the exercise for W events with only one or two jets, where the ti content of the data 
sample is expected to be very small. This is summarized in Fig. |7S| for DO and Fig. 7B for 
CDF. The background calculations reproduce the expected tagging rates in W events with 
low jet multiplicity. 

In this Section and in Section |V.E.4j , we have argued that the WQQ background is 
overestimated in both the CDF and DO analyses, yet the background predictions in the W+ 
1 and W+ 2 jets samples are in good agreement with the data. The reason for this is that 
the dominant background in these analyses is due to fake leptons, not WQQ. For example, 
if the CDF backgrounds were to be calculated using the best theoretical input for WQQ, 
the background estimate would have been reduced by only about 10%, independently of 
jet multiplicity (F. Abe et al, 1994a). This is well within the 20% systematic uncertainty 
assigned to the background calculation. 

The sizes of the excesses of events seen by the two collaborations are loosely consistent 
with expectations for ti. The DO collaboration sees an excess of 4.8 events; the number of 
expected tagged ti events varies between 4.4 and 1.2 for M top between 140 and 200 GeV/c 2 . 
The CDF excess is 7.6 tags, to be compared with a ti expectations of between 19 tags (M top 
= 140 GeV/c 2 ) and 2.6 tags (M top = 200 GeV/c 2 ). 

Despite the fact that the lepton tagging efficiencies are comparable in the two exper- 
iments, the signal-to-background is better for the DO analysis. This is mostly because 
backgrounds to the detection of low muons are lower in DO than in CDF. This fact is 
best illustrated by comparing Fig. [73] with Fig. [74. The muon tagging rate in DO is a fraction 
of 1% per jet; in CDF it is a fraction of 1% per track. There are three main reasons for 
the lower muon background in DO : (i) muons in DO have to traverse more steel than in 
CDF, so that background from hadronic punch-through is lower in DO than CDF; (ii) the 
DO detector is more compact, resulting in a lower probability for decays in flight of pions 
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FIG. 75. Comparison between the observed number of muon tags and the background expecta- 
tion in W+ jets event as a function of jet multiplicity for DO data. The requirement (see 



Table XI) has been removed for these data. Note that the horizontal axis is in terms of minimum 
jet multiplicity, e.g. the 2-jet bin includes all events with 2 or more jets. From Snyder, 1995a. 



and kaons; (iii) the momentum of the muon is measured after the decay-in-flight in DO, 
whereas in CDF, some average of the momentum of the parent pion and daughter muon is 
often measured in the drift chamber. The background rate per event, dominated by decays 
in flight and punch-through, is 1.8% in DO and 7.6% in CDF. 

In addition to these instrumental effects, the different choice of requirements between 
CDF and DO contribute to differences in signal-to-background. The DO pre-tag event sample 
with the requirement is expected to contain a higher fraction of ti events, for sufficiently 
high top mass. On the other hand, there is no requirement on in CDF, in order to 
maintain good efficiency for low top masses {M top < 120 GeV/c 2 ). For example, for M top 
= 200 GeV/c 2 , the ti contents of the DO and CDF pre-tag samples are expected to be 
approximately 9% and 6% respectively, see Table [XI|. 

Furthermore, the minimum requirement for lepton tags is lower in the CDF analysis 
than in the DO analysis (2 GeV/c vs. 4 GeV/c). Again, the CDF requirement was chosen 
to maintain efficiency for low top mass, where leptons from cascade decays have very low 
transverse momenta, (see Fig. [49|) . If the CDF requirement were to be raised to 4 GeV/c, 
the background would be reduced by a factor of ~ 1.8, while the tagging efficiency for 
M top > 130 GeV/c 2 would be lowered by approximately 20% only (F. Abe et al, 1994a). 
An excess of events in the CDF lepton tag analysis is also present when the minimum Py 
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FIG. 76. Comparison between the observed number of lepton tags and the background expectation 
in W+ jets event as a function of jet multiplicity for CDF data (Kestenbaum, 1996). 

requirement is raised to 4 GeV/c. In this case, 15 tags in 14 events are observed, with a 
background of 8.7 ± 1.8 tags (Kestenbaum, 1996). 
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3. Displaced-vertex tagging in CDF 

The most powerful method employed by the CDF collaboration to extract a top signal 
in the lepton + jet channel is to search for secondary vertices from 6-quark decay (vertex 
tagging). This is made possible with the precise tracking information obtained from a silicon 
vertex detector. 

CDF is the first, and so far the only, experiment to operate such a detector in a hadron 
collider. The first silicon vertex detector (SVX) at CDF was installed in 1992 prior to the 
beginning of Run la. The SVX (Amidei et ai, 1994) is a four-layer cylindrical detector 
which is 51 cm long. The four layers are at distances of 3.0, 4.2, 5.7, and 7.9 cm from the 
beamline. Axial microstrips with 60 /im pitch on the three innermost layers and 55 /im 
pitch on the outermost layer provide precision track reconstruction in the plane transverse 
to the beam. The single hit resolution is 13 /im, and the impact parameter resolution, for 
high momentum tracks is 17 /im. The SVX detector suffered significant radiation damage, 
and was replaced in 1993, during the accelerator shutdown period between Runs la and lb, 
by a very similar detector (SVX', Azzi, 1994) equipped with radiation-hard electronics. 

The vertex tagging algorithm used by CDF is based on reconstruction of displaced ver- 
tices using information from both the central tracking chamber and the SVX (or SVX'). At 
the Tevatron, pp interactions are spread along the beamline with standard deviation a ~ 30 
cm, so that the geometrical acceptance of the vertex detector is about 60% for pp interac- 
tions. This geometrical effect turns out to be the largest source of efficiency loss in tagging 
ti events. 

The position of the primary vertex is needed before searching for a possible secondary 
vertex. The transverse spreads of the colliding p and p beams result in a luminous region in 
the transverse plane which is Gaussian in shape with a ~ 36 /xm. The position and size of 
this region varies somewhat from store to store, and is monitored with an accuracy of order 
10 /zm. On an event-by-event basis, knowledge about the position of the primary vertex is 
improved by performing a fit using information from tracks consistent with originating from 
the primary vertex. The accuracy of the event-by-event determination of the position of the 
primary vertex in the transverse plane depends on the number of tracks available to the fit, 
and varies between 6 and 36 jum. 

The vertex tagging algorithm operates on combinations of at least two tracks with impact 
parameter at least three standard deviations different from zero. Constrained vertex fits are 
performed on these track combinations in an attempt to find one or more sets of tracks 
which are consistent with originating from a secondary vertex. Selection criteria are applied 
to reject tracks from decays of A and K°. Results of the vertex fit include the distance in the 
transverse plane of the secondary vertex from the primary vertex (L^) and its uncertainty 
((Ti X y). For a good secondary vertex which results in a 6-tag, L xy /ai xy is required to be > 3. 
The typical accuracy on the determination of the position of a secondary vertex is o~\ xy ~ 130 
/xm. This is much smaller than the distance travelled by a 6-hadron in a top event (typically 
a few mm, see Fig. [5U]), thereby allowing for efficient identification of secondary vertices 
from 6-hadron decay. There are two possible kind of vertex tags, positive tags and negative 
tags, (see Fig. [77|). Only positive tags are consistent with originating from the decay of a 
long-lived particle produced at the primary vertex. Negative tags, however, provide useful 
information on the performance of the 6-tagging algorithm, as we will discuss below. 
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(b) 



FIG. 77. An idealized sketch of vertex tagging in the transverse plane. The tracks from the primary 
vertex (a) are used to improve on the accuracy of the determination of the position of the primary 
vertex. Vertex (b) is a secondary vertex which is consistent with the hypothesis that it originates 
from the decay of a long-lived particle produced at the primary. This is a positive tag. On the 
other hand, vertex (c), although detached from the primary, is not consistent with the decay of 
a long-lived particle from the primary. Vertex (c) is an example of a negative tag, which is due 
to track mismeasurements. In a tt — > lepton + jets event, the tracks from vertex (a) in general 
originate from the underlying event or from the hadronization of the q or q from the decay of the 
W in t -» Wb, W -> qq. 

The capabilities of the CDF silicon vertex detector for detection of secondary 6-vertices 
are best illustrated in Fig. [78| which summarizes the CDF measurement of the lifetime of 
6-hadrons from the decay B — > J/ty + X; J j^f — > mi. The proper decay length (A) of the 
6-hadron is reconstructed from the position of the vertex as 



where Mj/^, is the J/^f mass, P T is the J/\l/ transverse momentum, and F is a Monte-Carlo 
determined correction factor which accounts for the undetected particles in the B — > J/^+X 
decay. The result of the lifetime fit, r# = 1.46 ± 0.06 ± 0.06 ps, is one of the world's most 
precise measurements of this quantity, and is consistent with the results from LEP. 

The CDF vertex tag efficiency in top events is e tag = (42 ± 5)%. This efficiency is denned 
as the probability of finding at least one (positive) displaced vertex within one of the jets in 
a tt event with > 3 jets. The vertex tag efficiency is a factor of two larger than the analogous 
efficiency of the CDF and DO lepton tag algorithms discussed in the previous Section. 

The value of et ag is determined from a Monte Carlo simulation of the response of the 
detector to tt events. In order to verify the reliability of the detector simulation, a number 
of studies of 6-tagging are performed in samples of b — *■ I events. The results of these studies 
are compared with expectations from a Monte Carlo simulation of this process, and the 



A = L 



xy 
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FIG. 78. CDF measurement of the 6-hadron lifetime from B — > J + X decays, (a) The distri- 
bution in A, the proper decay length, of data in the J/\E' sideband regions. The solid line shows 
the result of the fit. (b) The distribution in A, of data in the J/^ region. The curves are the result 
of the fit. The lightly shaded area is for B — > J/*S> + X events. The darkly shaded area is for the 
non-J/^ background. The large J/^f component at A ~ is due to prompt J/^> production and 
to x -> J/® decays. From F. Abe et al, 1993d. 



level of agreement found between the data and Monte Carlo is used to set the systematic 
uncertainty on e tag - The largest source of inefficiency is due to the long luminous region of 
the Tevatron (a = 30 cm). Since the vertex detector is only 51 cm long, approximately 40% 
of the interactions occur outside its geometrical coverage. 

Note that the tagging efficiency reported in the first CDF publications on vertex tagging 
in top events (F. Abe et al, 1994a) was considerably smaller, e tag = (22 ± 6)%. This 
value of the tagging efficiency was underestimated by 15%, due to a mistake in the Monte 
Carlo simulation. In addition, the tagging algorithm used in the earlier analysis has been 
substantially improved, and the performance of the SVX' in terms of efficiency and signal- 
to-background is somewhat better than that of SVX. 

In the lepton tag analyses described in the previous Section, the instrumental and WQQ 
backgrounds were estimated simultaneously under the assumption that the heavy flavor 
content of generic jets is the same as the heavy flavor content of jets in W events (Method 
I, see Section |V.E.4| ). As has been argued before, this method overestimated the size of the 
WQQ background. On the other hand, the Method II background calculation, which was 
also discussed at length in Section [V.E.4j , relies on a theoretical estimation of WQQ. To the 
extent that the instrumental background in the lepton tag analysis is larger than the WQQ 
background, numerically the Method I and Method II backgrounds are not very different. 
For the vertex tag analysis, however, this is not the case. Therefore, the CDF collaboration 
has chosen to calculate the background in this channel using Method II. 
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The instrumental background, i.e. the background due to false tags in light quarks or 
gluon jets, is estimated from the negative tagging rate as measured in a sample of generic 
jets. Since negative tags are almost exclusively due to tracking mismeasurements, and since 
these mismeasurements are equally likely to produce a positive or a negative tag, the negative 
tagging rate in generic jet is a measure of the mistag probability. This mistag probability is 
parametrized as a function of jet Ey, jet pseudo-rapidity, and track multiplicity, (see Fig [f9|). 
Based on this probability, the number of expected mistagged jets in the W+ jets samples is 
calculated by summing the mistag probabilities for all the jets in the sample. 
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FIG. 79. The tagging rate, defined to be the number of tagged jets divided by the number of jets 
having two or more tracks with > 2 GeV/c reconstructed in the SVX, as a function of (a) the 
jet Ey, and (b) the number of tracks associated with the jet. Shown are the rates for positive tags 
(circles) and negative tags (triangles) as measured in a sample of pp — > jets events. The positive 
tag rate is higher due to the fact that these jets contain a small fraction of heavy quarks (b and 
c). The positive tag rate is used to calculate the sum of backgrounds due to mistags and WQQ 
in the Method I background calculation. The negative track rate is used to calculate the mistag 
background separately in the Method II version of the background calculation. From the CDF 
collaboration, F. Abe et al, 1994a. Note that these tag rates are not for the final version of the 
vertex tagging algorithm used by CDF. 



The Method II WQQ backgrounds is calculated using theoretical expectations for the 
rates of WQQ. The fraction of W events containing a QQ pair is taken from theory, and 
multiplied by the number of observed W events to estimate the number of WQQ in the pre- 
tag samples. The product of this number with the expected WQQ tagging efficiency yields 
the expected tagged WQQ background. Note that the theoretical input is the fraction of W 
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Source 


W+ 1 jet 


W+ 2 jets 


W+ 3 jets 


W+ > 4 jets 


WQQ 


13.8 ± 11.1 


7.8 ± 6.2 


2.0 ± 1.6 


0.5 ± 0.4 


Mistags 


14.8 ± 3.0 


5.3 ± 1.1 


1.4 ± 0.3 


0.5 ± 0.1 


W + c 


15.3 ± 4.6 


4.2 ± 1.3 


0.9 ± 0.4 


0.2 ± 0.1 


Z — ► tt, Dibosons 


0.8 ± 0.3 


0.8 ± 0.3 


0.2 ± 0.1 


0.04 ± 0.2 


bb 


5.7 ± 1.4 


3.0 ± 0.8 


0.8 ± 0.2 


0.18 ± 0.04 


Total 


50 ± 12 


21.1 ± 6.5 


5.2 ± 1.7 


1.45 ± 0.43 


Tagged jets 


40 


34 


17 


10 



TABLE XII. Summary of Method II backgrounds and tags in the CDF vertex tag analysis 
(Carithers, 1995). The 27 tags in the > 3 jets sample occurr in 21 events. Assuming Stan- 
dard Model top production, one would expect between approximately 40 (M top = 140 GeV/c 2 ) 
and 5 (M top = 200 GeV/c 2 ) tagged tt events in the > 3 jets sample. 



events that include a QQ pair rather than the much more uncertain absolute rate prediction, 
see the discussion in Section |V.E.4j . The W + c background (see Fig. |5T]) is calculated in 
a similar manner, i.e. by multiplying the number of W+ jets events by the fraction of W 
events which are expected to contain a single c-quark, and by the tagging efficiency for these 
events. Other backgrounds, such as Z — > tt, dibosons, and bb are computed mostly from 
Monte Carlo. The background expectations and event yields are summarized in Table XII , 
and displayed in Fig. |8(] 

In the W+ > 3 jets sample, CDF finds 27 tagged jets in 21 events, with a background 
expectation of 6.7 ±2.1 tags. The power of the vertex tag algorithm is such that only 
1.9 ± 0.4 of these tags can be attributed to mistags, see Table [XII] . The six events with 
two tagged jets can be compared with four expected for the top + background hypothesis, 
and < 1 for background alone. Furthermore, six of the vertex tagged events also include 
a lepton tag. This is in much better agreement with expectations for top + background 
(about four events) than with background alone (about one event). As further evidence of 
a 6-contribution to this sample, the proper time distribution for tagged jets is also found 
to be consistent with expectations from 6-jets, see the inset of Fig. ^Q. As we will show in 
Section |VIII.E| , the size of the excess is consistent with expectations from tt for a top mass 
in the neighborhood of 160 GeV/c 2 . 

The Method II background calculation in the W + 1 jet sample is in good agreement with 
the data, providing a very important check of the reliability of the background estimation. 
We note that the Method I background estimate in the W + 1 jet sample yields 80 ± 10 
tags, in clear disagreement with the observed 40 tags. In the W+ 2 jets sample, the Method 
II background calculation is lower than the number of tags observed in the data. However, 
approximately 5 tags from tt are expected in this sample based on the excess of tags seen 
in the higher jet multiplicity samples. After accounting for these events, the background 
calculation in the W+ 2 jet sample is then found to be in satisfactory agreement with the 
data. 
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FIG. 80. Number of events before vertex tagging (circles), number of tags observed (triangles), and 
expected number of background tags (hatched) as a function of jet multiplicity. The inset shows the 
secondary vertex proper time distribution for the 27 tagged jets in the W+ > 3 jets data (triangles) 
compared to the expectation for 6-quark jets from tt decays. From the CDF collaboration, F. Abe 
et al., 1995a. 



4. Summary and cross checks of the tagging background calculation on Z+ jets 



Samples of Z+ jets in principle provide an ideal testing ground for the calculation of 
tagging backgrounds in the W+ jets sample. The properties of jets in W and Z events are 
very similar, due to the similarities in the W and Z production mechanism. One exception 
is that whereas in W events QQ pairs can be produced only through gluon splitting, see 
Fig. [52], in Z events additional mechanisms are also expected to contribute (e.g. gg — ► Zbb, 



see Fig. ^T|). A study of expected tagging rates in W or Z + 4 jets, including a model that 
approximates experimental efficiencies and backgrounds, indicates that the probabilities of 
tagging a W or a Z event are expected to be the same within 10-15% (Barger et al, 1994). 
Unfortunately, the combined cross-section times branching ratio for pp — > Z — > 11 is 
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FIG. 81. Feynman diagram for gg — > Zbb. 
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Z+ 1 jet 


Z+ 2 jets 


Z+ > 3 jets 


CDF lepton and vertex tag background expectation 


17.5 


4.2 


1.5 


CDF lepton and vertex tag data 


15 


3 


2 


DO muon tag background expectation 


0.97 ± 0.08 


0.35 ± 0.05 


0.09 ± 0.03 


DO muon tag data 












TABLE XIII. Comparison between the number of observed tags in Z+ jets data and the back- 
ground expectations. The CDF results are from Gerdes, 1995, and correspond to an integrated 
luminosity of 67 pb _1 . The DO results are from Abachi et al., 1995d, and are based on the Run la 
data set only, with an integrated luminosity of 13.5 pb _1 . 



an order of magnitude lower than the cross-section for pp — > W — > lv, giving only limited 
statistics in the Z channel. Nevertheless, it is still instructive to compare the yield of tags 
in the Z sample with the results of the background calculations. This is especially true in 
light of the fact that two out of the three Z+ > 3 jets events collected by CDF during run 
la contained one vertex-tagged jet (F. Abe et al, 1994a). 

The results of the tagging algorithms on Z+ jets events are displayed in Table [XIII| . With 
the higher luminosity data sample, the number of tags in the Z+ jets data is fully consistent 
with the background expectations, within the limited statistics. There is no evidence for an 
anomalously high tagging rate in Z events. 

To summarize, the CDF and DO background calculations have been proved to be reliable 
and, within the statistical and systematic uncertainties, they have been shown to be able to 
account for the rate of tagged jets in both Z+ jets and W+ 1 and 2 jets. There is significant 
evidence for the presence of an excess of 6-jets in the W+ > 3 jets sample. The largest 
excess is seen in the CDF vertex tag analysis, which has the highest efficiency and the best 
signal-to-background. The statistical significances of these excesses, which are consistent in 



size with the expected tt contribution, will be discussed in Section [VIII. D 



It is very natural to attribute these excesses to a ti component in the data, although 
they could also be due to some source of W+ heavy flavor production beyond standard 
QCD processes. As we will show in Section |V111.C| , additional evidence for the existence of 
the top quark can be obtained by kinematic studies of lepton + jets events. Furthermore, 
studies of invariant masses in lepton + jets events show evidence for both t — ► Wb; W Iv 
and t — > Wb; W — > qq, see Section |IX|. 
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Leptons 


Jets 








Channel 


Er(e) 


Pt(m) 








Hx 


A 


e + jets 


> 20 




> 4 


> 15 


> 25 


> 200 


> 0.05 


fjL + jets 




> 15 


> 4 


> 15 


> 20 


> 200 


> 0.05 



TABLE XIV. DO kinematic requirements for the standard event selection (energies in GeV, mo- 
menta in GeV/c). From Abachi et al., 1995b. 



C. Lepton + jets 

Both CDF and DO have performed analyses based on event shapes or kinematic variables 
rather than 6-quark identification in order to increase their tt acceptance. They have both 
used some form of a powerful discriminator between W+ jets background and top signal, 
which is the scalar sum of the transverse energies of the jets and in some cases leptons. DO 
has performed an analysis on a data set complementary to the lepton + jet + 6-tag data 
set. That is, it starts with the sample of lepton + jet events in which a 6-tag is not found. 
The results from this counting experiment are used in calculating the significance of the top 
signal. CDF on the other hand, performs two separate kinematic analyses on the complete 
lepton+jet data sample which includes 6-tagged events. These analyses are however not used 
in calculating the significance of the top signal both because their results are correlated with 
the b-t&g result, and because kinematic analyses depend to a greater extent on Monte Carlo 
generation details and theoretical assumptions. 

In comparing the data with the theoretical expectations for W+ jets, the CDF and 
DO analyses use the VECBOS Monte Carlo. As discussed in Section |\TEUl VECBOS 



is a leading order QCD parton Monte Carlo. Models of the underlying event and of jet 
fragmentations have been added to VECBOS by both collaborations to allow for comparisons 
with experimental data. 



1. DO lepton + jet kinematic analysis 

The requirements for this analysis are given in Table |X1V|. This analysis requires that 



there be at least four jets in the event, in contrast to the 6-tag analysis which includes 
events with three jets. The two main kinematic requirements are that the scalar sum of the 
transverse energies of the jets, H^, be > 200 GeV and that the event be aplanar such that, A 
> 0.05. Aplanarity, A, was defined in Section [V.E.3| as A = 3/2Ai, where Ai is the smallest 



eigenvalue of the matrix = Y^PaPb/HP 2 -, where Pi are the cartesian components of 
momentum of the parton, P is the magnitude of the three-momentum, and the sum is over 
all partons. In the DO analysis, A is calculated from the jets in the event. Note that the 
definition of Ht in the DO lepton + jets analysis is different than the definition of in the 
dilepton analysis, where the E^ of the highest transverse energy electron was also included, 
see Section [VIII. A| . Fig. |2| shows the distribution H T for top quarks produced by Monte 



Carlo and expectations for principal backgrounds. There is expected to be a clear separation 
between signal and background. 
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FIG. 82. Expected distributions in DO for top Monte Carlo events with M top = 200 GeV/c in 
the lepton + jets + no 6-tag sample (solid) and the principal backgrounds to this channel (dashed). 
From Abachi et al., 1995b. 



This analysis must rely on the correctness of both the Monte Carlo kinematics and the 
energy scale of the calorimeters. To check the combined effect, comparisons are performed 
between data and Monte Carlo distributions of H-r for two samples dominated by back- 
ground, W+ > 2 jets and W+ > 3 jets. Fig. |83| shows the data and VECBOS Monte 
Carlo prediction for these two samples. These samples may contain some top events, but 
are dominated by W+ jets background if the theoretical expectations for the ti production 
cross-section are correct. The agreement between the data and the Monte Carlo prediction 
is quite good. 

The acceptance, including branching ratio (24/81, see Table |T|), for the DO ti kinematic 
search ranges from 0.50% for M top = 140 GeV/c 2 to 1.70% for M top = 200 GeV/c 2 . The 
expected background, the expected signal and the data are shown in Table |XV| , and Fig. [84 . 

The two main sources of background in this analysis are QCD events in which one jet 
fakes the lepton signature and, especially, W+ > 4 jet events. The QCD background is 
calculated by studying a sample of events containing electromagnetic clusters that fail the 
lepton selection requirements. The estimation of the W+ jets background is clearly a crucial 
issue. This background is estimated in two ways (Abachi et al, 1995d; Grannis, 1995). In 
the first method, the number of W+ > 4 jets events in the data is estimated by extrapolating 
from the number of W+ > 2 and 3 jets events assuming a scaling law, i.e. assuming that the 
ratio N n /N n _i is independent of n, where N n is the number of QCD W+ > n jets events, 
(see Fig. |5| and Fig. |86]). Then, the number of W+ > 4 jets events expected to satisfy the 
aplanarity and Hr requirement is obtained from N4 and the shape of the VECBOS Monte 
Carlo distribution for W+ > 4 jets in the A vs. Hr plane, (see Fig. |34j). An alternative 
estimate is obtained by fitting the number of events in the four regions of the A vs. Ht 
plane indicated in Fig. |84] to contributions from ti, W+ > 4 jets and QCD events, where 
the shapes of the ti and W+ > 4 jets components are taken from Monte Carlo. The two 
methods are in reasonable agreement. The total expected background is 1.9 ± 0.5 events 
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FIG. 83. DO distributions for (a) e + > 2 jets, and (b) e + > 3 jets The curves are expectations 
from the VECBOS Monte Carlo normalized to the data. From Snyder, 1995a, for an integrated 
luminosity of 48 pb _1 . 



for the scaling method, and 2.6toJ ^ 0-5 events for the fit in the A vs. H T plane (Grannis, 
1995). 
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Events 


e + jets 


Expected background 
Expected for M top = 140 GeV/c 2 
Expected for M top = 200 GeV/c 2 
DO data 


1.22 ±0.42 
4.05 ± 0.94 
1.8 ±0.31 
5 


fi + jets 


Expected backgrounds 
Expected for M top = 140 GeV/c 2 
Expected for M top = 200 GeV/c 2 
DO data 


0.71 ±0.28 
2.47 ± 0.68 
0.95 ±0.24 
3 



TABLE XV. DO Lepton + jets results. From Abachi et al, 1995b. The integrated luminosities 
are 48 pb~ x for e + jets and 44 pb^ 1 for n+ jets. The ti expectations are normalized to the a(ti) 
calculation of Laenen, Smith, and van Neerven, 1994. 
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FIG. 84. Aplanarity (^4) versus H<r for single-lepton events for data, tt Monte Carlo, multijet 
background from data (with an effective luminosity = 60 x data luminosity), and background from 
W+ jets VECBOS Monte Carlo. From the DO collaboration, Abachi et al, 1995f. 
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W events : Run 1a + 1b data and Monte Carlo(VECBOS) 
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FIG. 85. Inclusive jet multiplicity spectrum for W — > eu + jets events for several jet energy 
thresholds. Data are shown by the solid symbols. Monte Carlo predictions are shown by the open 
symbols. From the DO collaboration, Abachi et al, 1995f. The integrated luminosity is 48 pb _1 . 
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FIG. 86. Inclusive jet multiplicity spectrum for Z — ► ee+ jets events. From the DO collaboration, 
Snyder, 1995a. The integrated luminosity is 56 pb _1 . 
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Sample 


Threshold 




E T (GeV) 


\Vjet\ 


Events 


(!mTtvrnl 


low 


=3 
Veto jet 4 


> 8 

> 8 


< 2.4 

< 2.4 


814 


Control 


high 


=3 
Veto jet 4 


> 15 

> 8 


< 2.0 

< 2.4 


104 


Signal 


low 


> 4 


> 8 


< 2.4 


267 


Signal 


high 


> 3 

> 1 


> 15 

> 8 


< 2.0 

< 2.4 
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TABLE XVI. CDF definition of the two control and two signal samples in the H analysis. The 
third, fourth and fifth columns list the criteria placed on the jets in each event. The integrated 
luminosity is 67 pb _1 . 



2. CDF lepton + jets kinematic analysis 

Two kinematic analyses of lepton + jets data have been published by the CDF collab- 
oration. The first analysis (H analysis, F. Abe et al, 1995c) is similar to the DO analysis. 
The second analysis (cos8* analysis, F. Abe et al, 1995b) relies on the fact that jets in tt 
events are more central than jets in W+ jets. 

The variable H in CDF is denned as the scalar sum of the lepton transverse momentum, 
the JZl T (i.e. the of the v) and the E^'s of all jets with > 8 Gev and pseudo-rapidity 
\v\ < 2.4. Note that while the event selection uses uncorrected jet transverse energies, H is 
calculated using corrected energies. The difference between the variables H used by CDF 
and Ht used by DO is the inclusion of the measurements of % T and lepton momentum in H. 

The H analysis is performed by comparing the observed H distribution in the data 
with prediction from the W+ jets and tt Monte Carlos. A deviation from the W+ jets 
expectation would then signal the presence of tt in the data. Note that this approach is 
totally independent of theoretical expectations for the W+ jets cross-section since it relies 
only on the predicted shape of the H distribution. In order to check the reliability of the 
theoretical prediction, the data H distribution is compared with theoretical predictions in 
four different samples, two control samples and two signal samples, see Table [XVI. 

The two control samples reject events with a fourth jet. These samples are designed to 
be dominated by background, and are a good testing ground for the VECBOS W+ jets 
calculation. There are less than 10 ti expected in each sample. The first control sample has 
a low requirement on the transverse energy of the jets, to further increase the proportion of 
W+ jet background relative to top events. The two signal samples require the presence of a 
fourth jet. The first signal sample is obtained with a low transverse energy threshold on the 
jets (8 GeV); in the second signal samples, which is expected to be enriched in tt events, at 
least three of the four jets must have transverse energies greater than 15 GeV. 

Data-Monte Carlo comparisons for the control samples and the signal samples are dis- 
played in Fig. |87] and Fig. The shape of the H distribution in the control samples agrees 
with the W+ jets VECBOS Monte Carlo, while the two signal samples are not consistent 
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FIG. 87. Comparison of the H distribution for the control data samples and the VECBOS Monte 
Carlo prediction, (a) W + 3 jets events passing the low Et threshold requirements; (b) W + 3 jets 
events passing the high Ey threshold requirements. The VECBOS prediction, including a 1% ti 
contribution for (a) and 10% ti contribution for (b), has been normalized to the data. From the 
CDF collaboration, F. Abe et al, 1995c. The integrated luminosity is 67 pb _1 . 
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FIG. 88. Comparison of the H distribution for the signal data samples and the VECBOS Monte 
Carlo prediction, (a) W+ > 4 jets events passing the low E^ threshold requirements; (b) W+ > 4 
jets events passing the high Ej- threshold requirements. The VECBOS prediction is normalized to 
a fit to the sum of ti and W + 4 jets Monte Carlo predictions. From the CDF collaboration, F. 
Abe et al., 1995c. The integrated luminosity is 67 pb _1 . 
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FIG. 89. The binned likelihood fit of the high threshold signal sample (solid line) to a linear 
combination of the VECBOS W + 4 jets and HERWIG ti Monte Carlo prediction. The H distri- 
bution of 6-tagged events is shown in the shaded histogram. From the CDF collaboration, F. Abe 
et al, 1995c. The integrated luminosity is 67 pb _1 . 

with W+ jets production alone. The differences between the high threshold signal sample 
and the VECBOS Monte Carlo is too large to be explained entirely by a systematic error in 
the experimental energy scale. 

The shape of the H distribution for ti events is expected to be a function of M top . 
Fitting the data to a linear combination of W+ 4 jets and tt Monte Carlo prediction yields 
M top = 180 ± 12(stat.)tj;5(syst.) GeV/c 2 . This value for M top is in good agreement with 
both the expected value from the Standard Model and with the directly measured value, 
see Section [IX|. Furthermore, the ti cross-section extracted from this measurement is also 
in good agreement with the cross-section extracted from studies of dilepton and 6-tagged 
events, see Section |V111.E| . 



As a further evidence for tt, the events containing 6-tags are highlighted in Fig. |89. 
These events are concentrated at high values of H, as would be expected for ti. Note that 
in background events there is a small dependence on the b-tag probability as a function of 
jet Et, and hence as a function of H, which biases the tagged events to high H. However, if 
one assumes that there are no ti events in the sample, this bias is not sufficient to account 
for the concentration of tagged events at high values of H. 

The second CDF kinematic analysis of lepton + jets data, (cos^* analysis, F. Abe et 
al, 1995b), is based on the fact that jets in top events are expected to be more central 
than in W+ jets events (Cobal, Grassman, and Leone, 1994). This analysis is also based 
on control and signal samples, see Table [XVI]] . The only difference in the two samples is in 



the requirement on the cosine of the polar angle (9*) of the jets in the center of mass frame 
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TABLE XVII. CDF kinematic requirements for the cos#* kinematic analysis. Cos#* is the cosine 
of the polar angle of the jet in the center of mass frame of the lepton + neutrino + jets system. 
The cos#* requirement applies to the three highest E^ jets only. Unlike all other CDF analyses, in 
this analysis jet energies are corrected at the event selection stage. Energies are in GeV, momenta 
are in GeV/c. From F. Abe et al., 1995b. The integrated luminosity is 67 pb _1 . 



of the two incoming partons. Events in the signal sample have three high central jets, 
events in the control sample have at least one jet that is emitted in the forward or backward 
region. Monte Carlo studies indicate that top events should be approximately equally split 
between signal and control samples, whereas of order 75% of W+ jets events should be in 
the control sample. Thus, the signal sample should be enriched in top quarks with respect 
to the control sample. 

An event-by-event Monte Carlo-based relative likelihood for the tt versus the W+ jets 
hypothesis is constructed for all events based on the E T of the second and third highest 
E^ jets in the event. Distributions of log-likelihood for data and Monte Carlo, for both 
signal and control sample, are displayed in Fig. |90| and Fig. |9Tj. There is a clear excess of 
events at high log-likelihood in the signal sample. The probability that the shape of the 
data log-likelihood distribution can be explained by W+ jets alone is < 0.26%, assuming the 
VECBOS model of W+ jets production. Furthermore, the 6-tagged events are concentrated 
in the high log-likelihood region, as one would expect for tt. 



3. Summary of lepton + jets kinematic analyses 

The CDF and DO collaborations have performed extensive studies of the kinematic prop- 
erties of lepton + jet events. The W+ jets data are inconsistent with the expectations from 
leading order QCD as implemented in the VECBOS Monte Carlo program, for high jet 
multiplicities and high jet transverse energies. On the other hand, there seems to be sat- 
isfactory agreement between data and theory in regions of low jet multiplicity and/or low 
jet transverse energies. The conclusion that can be drawn is that either the available QCD 
calculation is incorrect or incomplete, or that there is a source of events beyond standard 
QCD production of W+ jets. The most natural interpretation of the data is to ascribe this 
discrepancy to a tt component in the data. The correlation between high jet transverse 
energies and fe-tags reported by the CDF collaboration strengthens this conclusion. 



As was mentioned in Section |V.E.3| , samples of Z+ jets events are ideal to test the 
predictive power of the vector-boson + jets QCD calculation. A number of such tests 
have been performed, and they seem to yield results consistent with the QCD calculation. 
However, the statistics in the Z samples are not sufficient to decisively validate the QCD 
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FIG. 90. The distribution in natural log-likelihood for events in the signal sample, (a) Monte Carlo 
expectations for W+ jets (VECBOS) and ti (HERWIG); (b) Data. The lightly shaded histogram 
shows the 6-tagged events. The darkly shaded histogram shows events with two 6-tags. From the 
CDF collaboration, F. Abe et ai, 1995b. The integrated luminosity is 67 pb _1 . 
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FIG. 91. The distribution in natural log-likelihood for events in the control sample, (a) Monte 
Carlo expectations for W+ jets (VECBOS) and it (HERWIG); (b) Data. The lightly shaded 
histogram shows the 6-tagged events. The darkly shaded histogram shows events with two 6-tags. 
From the CDF collaboration, F. Abe et at, 1995b. The integrated luminosity is 67 pb _1 . 
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calculation in the kinematic region most relevant for the top search. 



D. Significance of the top signal 

As discussed in Sections |V111.A|, |V111.B|, and |V111.C|, the CDF and DO top searches in 



a number of different channels find excesses of top-like events over the background predic- 
tions. The significance of the excess for a given channel is defined as the probability that a 
background fluctuation would yield a number of events equal to or larger than the number 
of observed events. This probability is calculated by convoluting the Poisson probability 
for the fluctuation of the mean expected number of background events with its uncertainty, 
which is assumed to be Gaussian. 

When combining results from more than one channel, DO and CDF use different proce- 
dures. In the DO case, events from all channels as well as backgrounds from all channels are 
added, and the combined significance is defined as the probability that the sum of all the 
backgrounds fluctuates to give a total number of events greater than or equal the number 
of observed events. On the other hand, the CDF combined significance is defined as the 
probability that the product of the significances for the different channels be less than or 
equal to the measured value for this product. This is necessary since the CDF vertex-tag 
channel has so much better efficiency and signal to noise than the dilepton and (especially) 
the lepton-tag channels. In fact, if the DO prescription were to be applied to compute the 
total significance of the CDF excesses in the three channels, the combined excess would 
appear less significant than the excess seen in the vertex-tag channel alone, despite the 
fact that excesses are seen in the other two channels as well. When combining different 
channels, correlations are accounted for. For example, both the CDF lepton and vertex tag 
searches are sensitive to the size of the WQQ background, resulting in a correlation between 
significances in the two channels. These effects are taken into account using Monte Carlo 
simulations. 

The significances for the CDF and DO results are displayed in Table [XVII]] . These 
significances, which as discussed above are just probabilities, are also expressed in number 
of standard deviations (No. of a). The number of a corresponds to the point in a Gaussian 
probability function of mean zero and unit standard deviation where the integral of the 
probability function between that point and infinity is equal to the significance. 

The results from both experiments are such that the probability that the data can be 
explained as background fluctuations is exceedingly small. We note that, in the calculation of 
the combined significance of their top observation, the CDF collaboration very conservatively 
does not include the excess of top-like events seen in their purely kinematic analyses of 
lepton + jets data (see Section [VIII. C.2|) . Furthermore, as will be discussed in Section [IX|, 
both CDF and DO find additional evidence for a ti contribution to their lepton + jet data 
samples in the top mass reconstruction analyses. To summarize, the data from the counting 
experiments at the Tevatron provide overwhelmingly convincing evidence for the existence of 
new physics. We shall show in the following section that the invariant mass and kinematics 
of this new physics are consistent with a top quark of M top ~ 175 GeV/c 2 . 
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TABLE XVIII. CDF (F. Abe et al, 1995a) and DO (Abachi et al, 1995b) event yields, background 
expectations, and significances of the observed excesses of top-like events in different channels. In 
the CDF 6-tagging channels, we show both the number of events and, in brackets, the number of 
tags. There are 6 events which are tagged with both a vertex 6-tag and a lepton 6-tag in CDF. 
The two all combined entries for CDF data refer to 43 events and 56 objects, where an object is 
defined as either a dilepton event or a tag in a lepton + jets event. The CDF 6-tag backgrounds are 
expressed in number of tags; the CDF all combined background is expressed in number of objects. 
Tags in CDF dilepton events are not included anywhere in this table. 



E. Measurement of the pp — ► tt cross-section 

The ti production cross-section is calculated from the luminosity, the background- 
subtracted event yields, the acceptances, and the tt decay branching fractions for the various 
channels. The calculated cross-section depends on the top mass, since the acceptances are 
in general also a function of the top mass. 

The measurement is complicated by the fact that in some cases the expected number 
of background events depends on the cross-section itself. For example, in the CDF 6-tag 
analysis the background due to WQQ, which is quoted in computing the significance of the 
tt observation, is calculated under the assumption that the pre-tag W+ > 3 jets sample 
does not contain any top events. Having established the existence of the top quark, this 
background estimate needs to be revised, since the untagged sample does indeed include a 
tt component. This is accomplished by first calculating the cross-section using the overes- 
timated background contribution, then recalculating the background based on the value of 
the cross-section, and iterating until the result becomes stable. 

The CDF and DO tt production cross-sections, are listed in Table |X1X|, and displayed in 



Fig. As discussed at the beginning of this Section, the cross-section calculation depends 
on the assumed value of M top . Here we report a(pp — > ti) computed at the values of M top 
measured by the two collaborations, see Section |X|. The CDF and DO cross-sections reported 
here are obtained from data sets of 110 pb -1 and 100 pb -1 respectively. 

The top mass dependence of the cross-section measurement is 5a/5M t0 p ~ —0.05 
pb/GeV/c 2 (DO) and —0.04 pb/GeV/c 2 (CDF). For both experiments the dominant un- 
certainty on a{pp — > ti) is due to the limited statistics of the data samples. Also, for both 
experiments the separate cross-section measurements in the different channels are nicely 
consistent with each other. The CDF measurement is lower than, although still consistent 
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TABLE XIX. CDF (Tartarelli, 1996) and DO (Narain, 1996; Klima, 1996) measurements of a(pp -> 
tt) at a center-of-mass energy of 1.8 TeV. The integrated luminosities for these measurements 
are 110 pb -1 for CDF and 100 pb -1 for DO. Note that the most recent DO top mass value is 
Mt p = 170 ± 18 GeV/c 2 , see Section IX. A. 4. Cross section values for this value of M top for the 
individual channels have not been released by the DO collaboration. Here we list these cross-sections 
at Mtop = 180 GeV/c 2 to demonstrate the consistency between the measurements in the different 
channels. 

with, the cross-section value of 13.9^1 pb obtained from the analysis of the first 19 pb -1 of 
CDF data (F. Abe et al, 1994a). 

Both the CDF and DO tt cross-section measurements are in agreement with the the- 
oretical prediction. To summarize, both collaborations have not only reported excesses of 
top-like events in several different channels, but the sizes of these excesses are also consistent 
with each other. From these measurements, they have extracted values of a(pp — > tt) which 
are consistent, within the quoted uncertainties, with Standard Model expectations. 



IX. MEASUREMENT OF THE TOP QUARK MASS 

The CDF and DO top quark searches described in Section |V111| yield an excess of top- 



like events over the background expectations. In this Section we will present measurements 



of the top mass performed on a sub-set of these events. As was discussed in Section |T|, 
the mass of the top quark is one of the free parameters of the Standard Model, and its 
value enters in the calculation of radiative corrections to a large number of electroweak 
observables. Comparison of the measurement of M top with indirect determinations from 
electroweak measurements at LEP and SLC, as well as measurements of the W mass at 
the Tevatron and LEP200, and analysis of neutral current neutrino data, allow for rather 
stringent tests of the Standard Model. 

There are a number of possible ways to extract the top mass from the CDF and DO data. 
An indirect method is to simply compare the cross-section for the observed tt signal with its 
theoretical expectation, since the latter depends very strongly on the top mass, (see Fig. |32|) . 
However, it is more informative to measure both the top mass and the ti cross-section, and 
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FIG. 92. Combined CDF and DO top cross-section and mass measurements. (The mass measure- 
ments will be summarized in Section DC). Also shown, as a band, is the theoretical expectation for 
a(pp — ► tt) as a function of M top with its uncertainty (Laenen, Smith, and van Neerven, 1994). 



then compare with the theoretical prediction for a(tt). This approach results in a test of 
the QCD calculation of top production, and is sensitive to possible non-Standard Model 
production mechanisms for top quarks, see the discussion in Section [X|. 

Both CDF (F. Abe et al, 1994a and 1995a) and DO (Abachi et al, 1995b) have reported 
direct measurements of the top mass from lepton + jets data and these will be reviewed in 
Section IX. A . Other less direct methods to extract the value of M top from both the lepton + 
jets data and the dilepton data will be summarized in Section |1X..B| . Finally, in Section |LX.C| , 
we will present results of attempts to reconstruct the W mass from W — > qq in lepton + 
jets events. 



A. Direct measurements of the top quark mass from lepton + jets events 

The method used by both CDF and DO to measure the top quark mass is the standard 
method used in particle physics to measure particle masses, namely one measures the mo- 
menta of all the decay products, assigns masses to them, and then reconstructs the invariant 
mass of the original particle. In the case of top quarks, the decay products include neutrinos 
and hadron jets. The first example of the reconstruction of resonances using jets was UA2's 
observation of a bump consistent with the W and Z in the jet-jet invariant mass distribu- 
tion (Alitti et al., 1991a). In this case there was an an enormous background from QCD 
production of jets not involving a vector boson. The field of jet-spectroscopy is still in its 
infancy, although much progress has been made in the past year. 
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FIG. 93. Sketch of the tt decay chain in the lepton + jets channel. 

1. Constrained fits, combinatorics, and top mass resolution 

In the ti — > lepton + jets mode, one of the two top quarks (tjj decays semi-leptonicaly 
and the other one (t 2 ) decays hadronicaly, (see Fig. |93]) : 

h -> Wih ; W\ — > Iv 

t 2 -> W 2 b 2 ; W 2 —> qq 

The mass measurement is performed on the sample of events with a lepton, JZl x and four jets. 
The four jets are identified with the four quarks in the final state (6 1; b 2 , q and q). Without 
further information there is no way of knowing which jet originates from which quark. 
All possible combinations must be considered. A constrained fit to the ti hypothesis is 
then performed for each jet-quark assignment in a given event, assuming energy-momentum 
conservation at each vertex in the ti decay chain. The fit uses the following constraints : 

• The invariant mass of the jets assigned to the q and q is constrained to the W mass. 

• The $ T gives the transverse momentum of the neutrino. The longitudinal momentum 
of the neutrino, Pl„, is obtained by requiring the mass of the lepton and the neutrino 
to equal the W mass. This condition results in a quadratic equation for P L „, which 
has in general two distinct solutions. 
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• The invariant masses of the decay products of t\ and ti must be equal, i.e. the invariant 
mass of the three jets assigned to b 2 , q, and q (M 2 ) must be equal to the invariant 
mass (Mi) of the lepton, neutrino and fourth jet (61, see Fig. |93|). The top mass for 
the jet-quark assignment under consideration is then M top = Mi = M 2 . 

All of the components of momentum for the final state particles are measured, except Plv 
With one unmeasured quantity and three constraints, the fit is a two-constraint (2C) fit. 
The constrained fit also yields a x 2 for each combination which is a measure of the goodness- 
of-fit to the tt hypothesis. The fitted value of the top quark mass for a given combination 
is given by M top = Mi = M 2 at the point where x 2 is minimized. 

There are 4! = 24 possible ways of assigning the four jets to the four final-state quarks. 
Since there are two solutions for Pl v , this would result in 24 x 2 = 48 configurations. 
However, the interchange of jet assignments between the q and q from the W has no effect 
(see Fig. |93|), so that the number of truly distinct configurations is 48/2 = 24. If one or 
more of the jets is 6-tagged, the number of configurations can be reduced by allowing only 
configurations where 6-tagged jets are associated with the 6-quarks in the event. With one 
6-tagged jet, the number of combinations is twelve; with two 6-tagged jets, this number is 
reduced to four. 
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FIG. 94. Fit x 2 distributions for (a) best, (b) second-best, (c) third-best, and (d) the correct jet 
permutation. For e+jets events, from the ISAJET Monte Carlo, Mfo p =180 GeV/c 2 . The arrow 
shows the cut value. Note that the plot of correct permutations has less than half as many entries 
because only the cases in which a correct permutation could be found unambiguously from the 
Monte Carlo are included. From the DO collaboration, Snyder, 1995a. 
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The large number of possible jet-quark assignments, the poor jet energy resolution (see 
Section [V.B.2 and Fig. |32| ), and the effects of initial and final state gluon radiation, greatly 
complicate the top quark mass measurement. Monte Carlo studies indicate that there is 
often at least one combination with incorrect quark-jet assignments which yields a better 
fit to the ti hypothesis than the combination with the correct assignment, (see Fig. |94| ). 
Gluon radiation presents a problem because it can give rise to additional jets in the event. 
Both CDF and DO consider only the four highest Ey jets in the event, since inclusion of a 
fifth jet would increase the number of possible combinations by a factor of five. However, 
if one of these four jets is from gluon radiation, the constrained fit will be operating on the 
wrong objects. Because of the poor energy resolution, the goodness-of-fit variable is not 
very effective at eliminating this kind of event from the data sample. An additional effect 
of gluon radiation is that the lepton + 4 jets sample also includes ti events from the wrong 
decay mode, for instance events of the type ti — > lub rub can pass a lepton + 4 jets selection 
if both 6-jets are found, the tau decays hadronically and is reconstructed as a jet, and an 
additional gluon jet is present. 

The size of these effects depends somewhat on the details of the event selection. In the 
CDF analysis, Monte Carlo studies indicate that approximately 7% of ti lepton + 4 jets 
events are from the wrong decay mode, and of order 50% of the events have at least one 
of the four highest jets from gluon radiation. In the remaining events, the combination 
with the lowest x 2 corresponds to the correct parton-jet assignments only about one half of 
the times. 

Constrained fits to incorrect parton-jet assignments or to ti events from the wrong decay 
mode, in general yield incorrect values of M top - The CDF and DO groups have chosen to deal 
with the problem of combinatorics in slightly different ways. In the CDF analysis, only the 
lowest x 2 combination in a given event is considered. In the DO analysis, the top quark mass 
for a given event is taken as the x 2 -weighted average of all combinations (up to three) that 
have acceptable values of x 2 ■ The advantage of the DO approach is that the M top value for a 
given event is more stable under small changes in the measurements or the fitting procedure. 
These changes can cause the fit to converge to a different jet-parton configuration, and result 
in a value of M top which can be considerably different from the original. 

Wrong combinations result in a significant broadening of the expected mass resolution, 
see Fig. |95|. In CDF, the mass resolution for the correct jet assignment is expected to be 
12 GeV/c 2 . As a result of gluon radiation and wrong parton-jet assignments, the mass 
resolution is a factor of two worse, with significant non-Gaussian tails. Very similar results 
have also been reported by the DO collaboration. 

Despite the effects of gluon radiation and the high probability to choose a wrong combi- 
nation, the peak in the mass distribution of Monte Carlo events is not shifted significantly, 
(see Fig. |95|). This is partly because of order one-half of the wrong combinations involve 
interchange of one of the quarks from W decay (q or q) with the b from the hadronic top 
decay (bi, see Fig. For this class of events, the reconstructed top mass distribution is 
broader but still peaks at the correct value. 

It may be possible to reduce the effect of wrong combinations by including more informa- 
tion in the event fitting procedure. Examples of additional pieces of information that could 
be included are the expected angular distributions derived from the V-A structure of the 
top quark decay or the rapidity distributions of top quarks predicted by the ti cross-section 
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FIG. 95. Reconstructed top mass distribution for Monte Carlo events generated with the HERWIG 
Monte Carlo program, and simulated with the CDF detector simulation. The input value of the 
top mass in the Monte Carlo is Mt op = 170 GeV/c 2 . The solid line corresponds to the result of the 
constrained fit when requiring that one of the 6-jets is a b in the fit. The dashed histogram refers 
to the fit with the correct assignment for each of the jets. The width of the top quark, which for 
this top mass is T(t) « 1 GeV/c 2 see Fig. is negligible compared to the experimental resolution. 
From the CDF collaboration, F. Abe et al, 1994a. 



calculation (Kondo, 1991; Goldstein, Sliwa, and Dalitz, 1993). 

Another effect of the large number of combinations and the poor jet energy resolution 
is that the x 2 goodness-of-fit variable does not provide significant background rejection. (If 
it did, this variable would have been used to separate the top signal from the background). 
As discussed in Sections |V.E.1| , |Vlll.fj| , and |V111.(J| , the main background to the lepton + 
jets top signal is due to W+ jets production. In a W+ 4 jets event, it is almost always 
possible to find one jet-quark assignment with \ 2 low enough to be consistent with the ti 
hypothesis. Therefore, in order to perform the top mass measurement, both the background 
contamination of the event sample, and the mass distribution of the background events, 
need to be understood. 



2. CDF and DO top mass measurements 

The CDF and DO top mass measurements are based on the constrained fit procedure 
described in the previous Section. What is measured in these detectors is the energy and 
direction of the jets. In order to measure the top mass, the jet energies must be corrected 
to infer the original momentum of the partons. The correction procedure takes into account 
effects of non-linearities in the hadron energy response of the calorimeters, underlying event 
contributions to the energy of the calorimeter cluster, and gluon radiation outside the clus- 
tering cone. Special corrections are applied to lepton-tagged 6-jets, since additional energy 
is carried away by the neutrino emitted in the semileptonic 6-decay (typically a few GeV). 

A value of M top is calculated for each event using the constrained-fit procedure described 
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above. Events are rejected if the fit x 2 is inconsistent with the tt — > lepton + jets hypothesis. 
The top quark mass is extracted by performing a likelihood fit to the M top distribution 
for the remaining events. This M top distribution is fit to the sum of background and tt 
components, with the value of the top quark mass allowed to vary. The M top distribution 
for the dominant background, (W+ 4 jets events), is obtained by performing the same 
constrained fit procedure on a sample of events obtained from the VECBOS Monte Carlo 
and the detector simulation. The shape of the expected M top distribution, as a function 
of the mass of the top quark, is also taken from Monte Carlo event generators (ISAJET, 
HERWIG, or PYTHIA) and the detector simulation. In the fit, the size of the background 
contribution can be constrained, within errors, to its calculated value, or can be left free to 
float. In the latter case, the size of the background contribution returned by the fit serves 
consistency check of the procedure. 




FIG. 96. Mass distribution for DO top candidate events (histogram) compared with the expected 
mass distribution for 199 GeV/c 2 top quark events (dotted curve), background (dashed curve), and 
the sum of top + background (solid curve) for (a) standard and (b) loose event selection. 6-tagging 
information is used when available. From Abachi et al, 1995b. The integrated luminosity is ~ 45 
plr 1 . 

The results of the likelihood fit from the DO collaboration are shown in Fig. |96|(a) 
and |96](b) for two different event selections, tight and loose. In both cases the four jets 
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are required to have > 15 GeV. The tight selection includes requirements on H-r and 
aplanarity. The tight sample includes the eight lepton + four jets events from the analysis 
described in Section |V111.C| , as well as the 6 6-tagged events from Section |V111.B| . Only 11 
out of these 14 events yield an acceptable fit to the ti hypothesis. For the loose selection, the 
Ht requirement is removed and the aplanarity requirement is loosened. The Hr requirement 
selects events with high E-r jets, and introduces a bias that favors events with high recon- 
structed top mass. By removing this requirement, however, the background contribution 
is enhanced. The number of events increases from 14 to 27; 24 of these have at least one 
combination with good x 2 (Snyder, 1995a). As discussed in Section |V .±3| and illustrated in 
Fig. |33], in ti — > lepton + jets events there is a significant probability for two quarks in the 
final state to merge into a single jet. To minimize this effect, the jet-cone clustering radius 
(AR, see Section |V.B| ) in the DO mass reconstruction is changed from the 0.5 used in the 
selection of the top signal to 0.3. Likelihood fits to the two distributions result in top mass 
values of M top = 199±25 GeV/c 2 and M top = 199i2i GeV/c 2 for the tight and loose selections 
respectively (statistical errors only). The M top data distributions are not well described by 
the background hypothesis alone. This provides further kinematic evidence for the existence 
of the top quark. 

The CDF top observation in the lepton + jets channel was based on a sample of 6-tagged 
events with at least three jets of uncorrected E^ > 15 GeV, clustered with a cone of radius 
AR = 0.4. As it is clear from Fig. |93|, the constrained fit can only be applied to events with 
at least four jets. Therefore, the CDF mass measurement is performed on the sub-sample 
of events with a fourth jet. To maintain high efficiency, the E^ threshold on the fourth 
jet is lowered from 15 to 8 GeV (uncorrected). We stress that uncorrected jet energies are 
used at the event-selection stage only; for fitting purposes, all jet energies are corrected, see 
the discussion in Section [1X.A.3| . In Fig. we show the M top distribution of the pre-tag 
CDF lepton + 4 jets sample. Based on the CDF ti cross-section measurement, this sample is 
expected to be a mixture of approximately 30% ti and 70% W+ jets. The probability for the 
shape of the data M top distribution to be consistent with background only is approximately 
2%. Hence, just as in the DO case, this distribution provides additional evidence for the top 
quark. By demanding a &-tag, the bulk of events at low M top is removed, leaving a cluster 
of events between M top = 150 and 210 GeV/c 2 , see Fig. |98|. The result of the likelihood fit 
to these events is M top = 176 ± 8 GeV/c 2 (statistical error only). 

When including systematic effects, the mass values reported by the two collaborations are 
M top = 199i^(stat.)±22(syst.) GeV/c 2 (DO) and M top = 176 ± 8(stat.)±10(syst.) GeV/c 2 
(CDF). The results are consistent with each other. In the following Section we will discuss 
the important issue of systematic uncertainties. 
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FIG. 97. Mass distribution for the 88 CDF lepton + > 4 jets events with a good x 2 fit to the ti 
hypothesis, before the b-tag requirement (solid line). The darkly shaded histogram is for the 19 
events with a 6-tag. The expected W+ 4 jets contribution to the pre-tag sample is shown in the 
lightly shaded histogram. From F. Abe et al, 1995a. The integrated luminosity is 67 pb _1 . 
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FIG. 98. Mass distribution for the 19 CDF 6-tagged lepton + > 4 jets events with a good x 2 
fit to the ti hypothesis (solid). Also shown are the background shape (dotted), and the sum of 
background and ti Monte Carlo (dashed) for M top = 175 GeV/c 2 . The inset shows the likelihood 
fit used to determine the top mass. From F. Abe et al., 1995a. The integrated luminosity is 67 
pb- 1 . 
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3. Jet energy corrections and systematics on the M top measurement 

Although the data samples of top candidates are rather small, the size of the statistical 
error on the top quark mass measurement is already comparable to the size of the systematic 
uncertainties. It is expected that in the next decade much larger ti samples will be available, 
see Section |X|. The systematic uncertainties will then be the limiting factor in the precision 
of the top mass measurement. 

The dominant uncertainty in the top mass measurement is due to the uncertainty in the 
jet energy scale, i.e. to the transfer function (correction factor) that relates the measured 
jet energies to the energies of the original quarks from top decay. This uncertainty has 
two components, (i) an instrumental uncertainty related to the response of the calorimeter 
to hadrons, and (ii) an uncertainty in the understanding of the fragmentation and gluon 
radiation processes. In order to discuss these systematic uncertainties, we begin by describing 
the many steps in the jet energy correction procedure (F. Abe et ai, 1992d and 1993c; Abachi 
et ai, 1995d). 

The jet energy response is in general not uniform across the detector because, for example, 
of cracks at the boundaries between calorimeter modules. The first task of the correction 
procedure is to equalize response across the calorimeter. The size of the effect is measured 
in-situ by di-jet or photon balancing techniques. Di-jet balancing is performed on a sample 
of pp — > 2 jets events, with one jet restricted to a well-understood region of the calorimeter. 
Since the transverse energies of the two jets are expected to be equal, any inbalance as 
a function of the position of the second jet is a measure of the position-dependent non- 
uniformity of the calorimeter response. A similar study can be performed using pp 7+ jet 
events. In these events the accurate measurement of the of the photon can be compared 
to the measurement of the jet Ey. These effects can be measured, and hence corrected for, 
with high precision, because of the very large number of di-jet and photon-jet events that 
can be used for this study. 

The jet cluster will in general include energy deposited from particles unrelated to the 
parent parton, for example particles from the underlying event. The jet energy is therefore 
corrected by subtracting off the average underlying event deposition as measured in minimum 
bias events. For jets clustered with a cone radius AR = 0.4, this correction amounts to 
approximately 600 MeV for the CDF detector. Note that in ti events the amount of extra 
energy can be higher, because of cross-talk between the many jets in the final state. An 
additional, Monte Carlo based, correction for this effect has been developed by the CDF 
collaboration. A further correction needs to be applied in DO due to noise from radioactivity 
of the uranium plates used in the calorimeter. 

The next ingredient in the jet energy correction procedure involves understanding the 
absolute energy response of the calorimeter. The response to individual hadrons is measured 
in test-beams. In CDF the hadron response of the calorimeter as a function of momentum is 
also obtained from samples of isolated particles from pp collisions. The detector simulation is 
adjusted to incorporate information from these measurements. Finally, the jet energy scale 
is measured by simulating the calorimeter response to jets generated using a QCD-based 
model of jet fragmentation. 

The absolute energy scale can also be derived from photon balancing, since the energy 
of the photon is precisely measured in the electromagnetic calorimeter. The electromagnetic 
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FIG. 99. DO energy scale correction for jets as a function of jet transverse energy in the central 
and forward regions. Results are for jets reconstructed using a cone-size AR= 0.5. The dashed 
curves represent the error bands. From Abachi et al, 1995d. 



calorimeter is precisely calibrated by studying the distribution of E/P for electrons (CDF), 
or from reconstruction of the Z — > ee resonance (DO). As we will discuss shortly, the photon 
balancing technique simultaneously tests instrumental as well as gluon radiation effects. 

Finally, in order to go back to the original parton energy, an additional correction has to 
be applied for energy radiated outside the clustering cone, or carried away by low momentum 
hadrons swept away by the magnetic field (CDF). This correction is based on a Monte Carlo 
model of the QCD process. The overall jet energy correction factors are displayed in Fig. ^ 
and Fig. |100| for DO and CDF respectively. The observed jet energies are to be multiplied 
by these correction factors to obtain the corrected energy. 

Photon-jet balancing provides a powerful probe of the behavior of the overall jet energy 
correction function, (see Fig. |101|) . Unfortunately this test is not free from its own systematic 
uncertainties. Photon samples collected at the collider are contaminated at about the 50% 
level by two-jet events, with one of the jets fragmenting to a high momentum tt° or 77, 
which is misidentified as a single photon. In these cases the transverse energy of the photon 
candidate is not expected to equal the transverse energy of the recoiling jet, because of the 
presence of additional fragmentation hadrons in the jet that fakes the photon signature. 
Furthermore the balancing results can be affected by undetected low-energy initial-state 
gluon radiation. It is estimated that these effects introduce an uncertainty of order 5% on 
the determination of the absolute energy scale. Similar studies are performed using Z+ jet 
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FIG. 100. CDF energy scale correction for jets as a function of jet pseudorapidity. Results are for 
jets reconstructed using a cone-size AR = 0.4. (a) observed Et = 15 GeV; (b) observed Et = 30 
GeV; (c) observed Et = 50 GeV; (d) observed Et = 100 GeV. The cracks between calorimeter 
modules (at r) 0., 1.1, 2.4) are apparent. 



events, (see Fig. 102). The background effects which systematically limit the usefulness 



of photon-jet balancing are not present in this case; however, the available statistics are 
considerably lower. 

The correction procedure described above applies to gluon and light-quark jets. Monte 
Carlo studies show no significant difference for 6-jets, except in the case of semileptonic 
decays where an additional correction for the undetected neutrino needs to be applied. 

The uncertainty in the jet energy scale for both the CDF and DO mass analyses is 
estimated to be 10%. This value is somewhat larger than one would infer from photon-jet and 
Z-jet balancing studies, (see Fig. |101| and Fig. |102|) . There are however additional questions 
concerning the applicability of this correction procedure to the hadronic environment in 
top events which contribute to the 10% uncertainty estimate. The 10% jet energy scale 
uncertainty translates to a top mass uncertainty of 8 GeV/c 2 and 21 GeV/c 2 for CDF and 
DO respectively. The difference between the two experiments in the size of the reported M top 
uncertainties is not fully understood at this time. 

There are additional smaller systematic uncertainties on the top mass measurement. A 
related uncertainty comes from the Monte Carlo modeling of gluon radiation in tt events. 
As was discussed in the previous Section, gluon radiation broadens the reconstructed M top 
distribution. Since the top mass is extracted from a fit to the M top data based on the 



125 



300 



Mean =0.027 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 

(E T (photon) -E T (j et) )/E T (photon) 



FIG. 101. A = (Et (photon) - Et (recoiling-jet))/Er (photon). The jet transverse energy is 
corrected. Jets are clustered in a cone AR = 0.4. From the CDF collaboration, F. Abe et al., 
1994a. 

expected it M top distribution, different assumptions on its shape can result in shifts of the 
measured top mass. For example, differences between the ISAJET and HERWIG models 
result in shifts of 4 and 1 GeV/c 2 for DO and CDF respectively. Further uncertainties at the 
1-2 GeV/c 2 level are present in the CDF measurement from uncertainties in the background 
shape, which is taken mostly from the VECBOS Monte Carlo, as well as the details of the 
likelihood fitting technique. 

The CDF and DO M top measurements are the first examples of the application of jet 
spectroscopy techniques to the determination of the mass of an elementary particle. The 
understanding and control of the systematic uncertainties due to the jet energy measure- 
ments are expected to improve in the future. This will be crucial to allow for more precise 
measurements of the top mass. We will discuss a number of possible approaches in Section [X|. 

4. Updated CDF and DO top mass measurements 

As this review article was being completed, both the CDF (Tartarelli, 1996) and DO 
(Narain, 1996) top mass results have been updated. In both experiments, the statistical 
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FIG. 102. Difference between the transverse energy of the electron pair and the jet in {Z — > ee) + 1 
jet events. For Monte Carlo and data, with and without out-of-cone corrections, but including all 
other corrections. The minimum jet is 10 GeV. Jets are clustered in a cone AR =0.3. From 
the DO collaboration, Snyder, 1995a. 



errors have decreased due to the larger data samples and, more importantly, the systematic 
uncertainties have been reduced as a result of more detailed studies of the jet energy scales 
in the two experiments. Since numerically the DO result is somewhat different than the 1995 
result discussed in Section [IX.A.2| , we briefly include these updated results here. 

The new DO measurement is based on an integrated luminosity of 100 pb _1 . Besides 
the improvement in statistics, there are four differences between the old and new DO mass 
measurements: (i) the selection of the sample on which the mass fit is performed has been 
changed, (ii) an error in the out-of-cone jet energy corrections has been fixed, (iii) the 
HERWIG Monte Carlo is used instead of ISAJET, and (iv) the jet energy scale has been 
shifted downwards by 5%. 

The new sample is derived from the loose lepton + four jets sample, removing the 
aplanarity cut. The sample is then split into two pieces, those events with and without a 
low P T muon 6-tag. The untagged events must satisfy two additional requirements: (i) the 
transverse energy of the leptonic W must be greater than 60 GeV, and (ii) the pseudorapidity 
of the W must be in the interval ± 2. The longitudinal momentum of the neutrino from the 
W is chosen from two possibilities as the one having the lowest momentum. Furthermore, 
the events are required to pass a kinematic liklihood test. This was done in order to reduce 
background for the loose sample, while reducing the mass bias due to the previous 
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Source M top uncertainty 
Jet energy scale (4% ± 1 GeV) 7 GeV/c 2 

Different Monte Carlo tt generators 6 GeV/c 2 

Fitting 3 GeV/c 2 

Background uncertainty 2 GeV/c 2 

Total 10 GeV/c 2 



TABLE XX. Systematic uncertainties in the updated DO mass measurement. From Narain, 1996. 



Source 


M top uncertainty 


Jet energy scale (detector effects) 


3.1 GeV/c 2 


Soft gluon effects 


1.9 GeV/c 2 


Hard gluon effects 


3.6 GeV/c 2 


Different Monte Carlo tt generators 


0.9 GeV/c 2 


6-tagging bias 


2.3 GeV/c 2 


Background spectrum 


1.6 GeV/c 2 


Fit configuration 


2.5 GeV/c 2 


Likelihood method 


2.0 GeV/c 2 


Monte Carlo statistics 


2.3 GeV/c 2 


Total 


7.1 GeV/c 2 



TABLE XXI. Systematic uncertainties in the updated CDF mass measurement. From Tartarelli, 
1996. 



requirement. No further requirements are made on the sample with a 6-tag, except for 
additional improvements to the algorithm for removing Z — > ////, and for the application of 
a tighter low Pt muon selection. 

This analysis results in 30 events which satisfy the fit \ 2 i nve °f which have a 6-tagged 
jet. The background in this sample is estimated to be 17.4 ± 2.2 events; the fit is shown in 
Fig. |103| , and gives M top = 170 ± 15(stat.) ±10(syst.) GeV/c 2 . This is to be compared with 
the earlier result M top = 199l 2 i(stat.)±22(syst.) GeV/c 2 . The systematic uncertainties are 



summarized in Table PCX 



The updated CDF measurement, based on an integrated luminosity of 110 pb _1 , is 
performed in the same manner as described in Section |1X.A.2| . The result is M top = 176 ± 
6(stat.)±7(syst.) GeV/c 2 ; the systematic uncertainties are summarized in Table |XXI . 



The combined top mass measurement from the two Tevatron collaborations is M top = 
175 ± 8 GeV/c 2 . Here we have added the errors in quadrature and neglected correlations 
in the systematic uncertainties. These correlations are due, for example, to the modelling 
of tt production, and to the common assumptions made by the two collaborations in the 
determination of the jet energy scale. 



128 




FIG. 103. Mass distribution for the 30 DO lepton + > 4 jets events with a good \ 2 fit to the ti 
hypothesis, from the 100 pb _1 data sample. The solid line is a fit to top plus expected background; 
the dashed line is the expected background from VECBOS (W + jets) and QCD multi-jet back- 
ground; the dotted line is the fitted top contribution. The inset shows the liklihood distribution. 
From Narain, 1996. 
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B. Top mass from dilepton events and kinematic distributions 

The top mass can also be reconstructed, in a less direct way, from dilepton events (Dalitz 
and Goldstein, 1992; Kondo, 1988 and 1991). In a it dilepton event, both top quarks decay 
semi-leptonicaly : 

h -> W x bi ; Wi -> hvi 



t 2 -> W 2 b 2 ; W 2 -> / 2 ^2 

Because of the presence of two neutrinos, a direct event-by-event reconstruction of the top 
mass based only on the measurements of the momenta of the leptons and the jets is not 
possible. The system is underconstrained, as can be seen from a simple accounting of the 
degrees of freedom. To fully describe the event, one needs the momenta of all the quarks 
and leptons in the final state. The momenta of the charged leptons are measured, and 
those of the quarks are inferred from the jet energies. Six parameters are needed to fully 
describe the two neutrinos; only two measurements (from the two components of E! T ), and 
three constraints are available : Mass(Zi^i) = M^; Mass(72^2) = M^; Mass(ii^ifei) 
= Mass(/2^2&2)- This leaves 6—5 = 1 parameter undetermined. In order to measure the 
top mass, additional information need to be included. Possibilities include the Standard 
Model V-A expectations for the angular distributions in the top decay chain, or theoretical 
expectations for the kinematic properties of the produced t and i. 

The DO collaboration (Snyder, 1995b) recently reported a preliminary measurement of 
the top mass from dilepton events. The DO dilepton mass measurement uses five dilepton 
events. It is based on a data sample with a higher integrated luminosity than the DO dilepton 
top search described in Section [VIII. A| which yielded three candidate events. 



The top quark mass reconstruction for a single event goes as follows. By assuming a 
value of M top , and for a given jet-6-quark assignment, the momenta of v\ and v 2 are deter- 
mined up to a possible four-fold ambiguity from the two quadratic constraints M 2 (/iz/i) = 
M^, and M 2 (/ 2 z/ 2 ) = M 2 ^. For each configuration (two jet-6-quark assignments, and possible 
four-fold neutrino ambiguity), the event is fully reconstructed. A probability (Prob; ep ) is 
assigned to each configuration based on the energy of the leptons in the rest frame of the 
top quarks. This probability is calculated for the assumed top mass, and is based on the 
expected structure of the decay. The momentum fractions X\ and x 2 of the incoming quarks 
qi and q 2 in the reaction q\ + q 2 — > ti are also reconstructed from the invariant mass (M#) 
and momentum (Pfj) of the ti system. Modulo the effects of gluon radiation, x\ and x 2 
can be obtained from M# = ^/x\X 2 s and P t t = 0.5y/s(x\ — x 2 ), where -/a = 1.8 TeV is the 
center-of-mass energy of the pp collision. An additional probability (Prob x ) is assigned to 
the configuration, still as a function of the assumed top mass, based on the parametriza- 
tion of the parton distribution functions at momentum fractions x = X\ and x = x 2 , (see 



Fig. [19]). Because the parton distribution functions are decreasing functions of x, Prob x 
includes a correction to remove biases towards low values of M top . The total probability 
for the configuration is given by the product Prob; ep Prob x . The total probability for the 
event is the sum of the probabilities for all possible configurations. By varying the value of 
the assumed top mass, this procedure yields a probability distribution for each event as a 
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function of top mass. To include resolution effects, the procedure is repeated many times 
after smearing the measurements of the lepton and jet momenta according to their expected 
resolutions. The total probability distribution for a single event is then defined as the sum 
of the probability distributions for the smeared events. An event top mass value is defined 
as the value of M top at which the probability distribution is maximum (M pea k). Finally, the 
top mass is extracted from likelihood fits of the data M pea fc distribution to the superposition 
of the expected tt and background contributions, (see Fig |104|). 
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FIG. 104. Dilepton mass measurement from the DO experiment. The arrows show the M pea fc values 
for the five DO dilepton candidate events. The solid curve is the expected signal distribution for 
M top = 145 GeV/c 2 , the dashed curve is the expected background distribution. See text for details. 
From Snyder, 1995b. 

The preliminary DO result for this procedure is M top = 145±25(stat.)±20(syst.) GeV/c 2 . 
The systematic uncertainties include energy scale effects, Monte Carlo modeling, and un- 
certainties in the background contributions to the fit. Within the large errors, this value is 
consistent with both the CDF (M top = 176 ± 9 GeV/c 2 ) and DO (M top = 170 ± 18 GeV/c 2 ) 
measurements in the lepton + jets channel described in Section |IX.A.4 . 

A number of kinematic quantities in both dilepton and lepton + jet events have also 
been shown to be sensitive to the top mass (see for example Baer et al, 1990). It is then 
in principle possible to perform a measurement of the top mass by comparing kinematic 
distributions to expectations for ti. Some possibilities are the event transverse mass in 
dilepton events, the mean of the lepton-6 mass, the transverse momentum of 6-jets, or the 
total transverse energy, (see Fig. |105|) . 

These more indirect measurements are less statistically powerful than the direct mea- 
surement of the top mass described in Section |IX.A| , since the mass peak used in the direct 
measurement provides optimal discrimination between different M top hypotheses. However, 
indirect measurements are sensitive to different systematic effects. For instance, the combi- 
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FIG. 105. Predictions from the IS A JET pp -> tt Monte Carlo at y/s = 1.8 TeV, as a function of 
the top mass, for (a) the mean transverse momentum of 6-quarks, and (b) the mean of the scalar 
sum of the transverse momenta of the four quarks in lepton + jets events. The vertical error bars 
represent the expected RMS for these quantities. 



natorics problem which plagues the direct top mass measurement, would not play a role in 
most of these techniques. A measurement of the top mass based on the 6-quark transverse 
momentum spectrum or the 6-lepton invariant mass, is almost entirely independent of initial 
state radiation in ti events. 

The CDF collaboration has also reported measurements of the top mass based on kine- 
matic distributions in both dilepton events (Tartarelli, 1996) and lepton + jets events (F. 
Abe et ai, 1995c). The dilepton mass determination is based on a fit to the jet E T spectrum 
in these events and yields M top = 159l|2(stat.)±17(syst.) GeV/c 2 . The lepton + jet mea- 
surement is based on a fit to the the total transverse energy distribution (see Section |V111.C| ), 
and gives M top = 180±12(stat.)i]J(syst.) GeV/c 2 . Both these results are in agreement with 
the results derived by direct mass reconstruction techniques. 

The direct top mass measurement in 6-tagged lepton + jets events is, and is likely 
to remain, the method of choice for determining the top mass. Nevertheless, alternative 
methods such as the under-constrained dilepton mass reconstruction or methods based on 
kinematic distributions, serve as useful consistency checks. 
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C. Reconstruction of the W mass from hadronic decays in lepton + jets events 



In a tt — > lepton + jet event, one of the two W bosons from the top decay chain is 
expected to decay hadronically (W — > qq, see Fig. |93|). Both the CDF and DO collaborations 
have reported evidence for this decay mode by reconstructing a peak in the invariant mass 
distribution of two jets in the event. Hadronic W reconstruction is interesting for two 
reasons. First of all, the presence of such a peak in the data is further strong evidence 
that the excess of events over the background prediction is indeed due to tt production. 
Furthermore, the hadronic W peaks provide the most ideal calibration for the jet energy 
scale of the two experiments. Note that it is extremely difficult to get the scale information 
from inclusive W decays to two jets, because of the very high QCD di-jet background which 
would both mask the signal and saturate the data acquisition bandwidth. 



Dijet mass of untagged jets in events 
with a b tag and a second loose tag 

CDF PRELIMINARY 100 pb" 



For the eigh; events in peak 
Mear 81 ,4 ± 3.2 GeV 
Sigma 9.1 ± 6.0 GeV 



errors statistical only 
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FIG. 106. The di-jet mass of untagged jets in data events with a double 6-tag (histogram), com- 
pared to the expected signal from the top Monte Carlo (dotted). The integrated luminosity for 
this data set is 100 pb -1 . From the CDF collaboration, Yao, 1995. 



The hadronic W has been reconstructed by the CDF collaboration (Yao, 1995). The 
method uses W+ > 4 jet events with two 6-tagged jets. One of the jets is tagged using 
either of the standard CDF lepton or vertex tagging algorithms. In order to increase the 
statistics of the sample, the second jet is tagged using a looser vertex tag algorithm. Such a 
double-tag requirement considerably reduces many sources of background, and the tt purity 



of the sample is very high. Just as in the constrained fit procedure described in Section [IX. A| , 
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only the four highest Et jets are considered. With this restriction, in a doubly 6-tagged ti 
event there is no ambiguity in assigning jets to the hadronic W -decay. The invariant mass 
distribution of the two untagged jets in this sample is shown in Fig. |106| . In eight of the 
ten events, the invariant masses of the two untagged jets cluster tightly around Mjj ~ 80 
GeV/c 2 . This provides overwhelming evidence for the presence of hadronic decays of W 
bosons in the lepton + jets + 6-tag sample. The two outlying events are most likely due to 
cases where one of the untagged jets considered to be from the W is really a gluon jet from 
initial or final state radiation. 

The DO hadronic W mass reconstruction procedure (Strovink, 1995) in lepton + jets 
events is much more complicated, since the 6-tag in DO is not as efficient as in CDF. The 
analysis uses lepton + > 4 jets data, and considers only the four jets with the highest 
Et- Just as in the constrained fit procedure (see Section [IX. A| ), the neutrino longitudinal 
momentum (Plu) is obtained by constraining the lepton-neutrino invariant mass to the W- 
mass, and only the solution with the smallest \Plv\ is considered. No constraint is placed on 
the di-jet mass, so that there are four different ways of partitioning the event into (luji) and 
(J2J3J4), unless one of the jets is 6-tagged, in which case the number of combinations is two. 
A two-dimensional scatter-plot of di-jet mass versus top mass is filled for each combination 
with weight proportional to exp(— x 2 /2), where x 2 = ln 2 (mass(Zz/j 1 )/mass(j 2 j3j4))- The 
normalization is chosen in such a way that the weights for all combinations in a given event 
sum to unity. The top mass is defined as the average of mass(lvji) and mass(j2j3j4) for 
electron + jets events. For muon + jets events, the top mass is calculated by taking a 
weighted average of the two, with the values of hadronic and leptonic masses weighted in 
the ratio 60 : 40. The di-jet mass is defined as follows. If one of the jets in (J2J3J4) is tagged, 
the di-jet mass is the invariant mass of the two other jets. Otherwise, the most energetic jet 
in the top rest frame is used as the b. But if no jets in {jzjsj^} are tagged, and their energy 
in the top rest frame are such that (E2 — E3) < (E3 — E4), then the permutation is plotted 
twice, with equal weights, with di-jet mass chosen as mass(j3j4) or mass(j2j4)- Studies 
of ti Monte Carlo events, indicate that this procedure results in a peak in the di-jet-mass 
vs. top-mass scatter plot at di-jet-mass = M w and top-mass = M top . On the other hand, 
background events are expected to peak at lower values. Projections from the di-jet-mass 
vs. top mass scatter plot are shown in Fig. |107| . The top mass and di-jet mass distributions 
are peaked around 180 and 80 GeV/c 2 respectively. The probability for these distributions 
to be consistent with the background-only hypothesis is 1.3%. 

In conclusion, both collaborations have shown evidence for hadronic W decays in their 
respective top event samples. With the aid of vertex tagging, the CDF W — > qq peak is 
very clean and straightforward to understand. Future higher statistics samples of W — > qq 
in top events will provide a very important calibration to the top mass measurement, see 
Section [X]. 
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FIG. 107. Distributions of (a) reconstructed top quark mass, and (b) reconstructed di-jet-mass. 
The reconstructed top-mass is plotted only for di-jet-mass > 58 GeV/c 2 . The reconstructed di- 
jet mass is plotted only for top mass > 150 GeV/c 2 . Distributions are shown for data, sum 
of background and HERWIG tt Monte Carlo, background alone, and background normalized to 
match the area of the data. From the DO collaboration, Strovink, 1995. 
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X. FUTURE PROSPECTS 



One of the goals of particle physics in the next decade is to perform detailed experimental 
studies of the properties of the top quark. Since the mass of the top quark is of the same 
order as the scale for electroweak symmetry breaking, it is possible that new physics effects 
will manifests themselves in the top sector. While we do not know what these new effects 
might be, it is clear that top quark physics represents an opportunity to uncover physics 
beyond the Standard Model. 

In this Section we will concentrate on prospects for top physics at the Tevatron. We 
will begin by briefly summarizing plans for the upgrade of the accelerator and the detectors 
in Section [X.A| . Detailed studies of the possibilities for the Tevatron top physics program 
are well under way (Amidei and Brock, 1996). Prospects for improving the accuracy of the 
top mass measurement will be discussed in Section [X.B| ; the potential for study of the Wtb 
vertex will be addressed in Section [XX]; finally, additional tests of the Standard Model, as 



well as searches for new physics involving the top quark will be discussed in Section K.D 



A. Accelerator and detector upgrades 

The CDF (CDF collaboration, 1995) and DO (Tuts, 1996) detectors at the Tevatron will 
be undergoing major upgrades in the next few years. The most significant improvements for 
top physics will be the installation of a magnet for charged particle momentum determination 
in DO and new 3D silicon vertex detectors in both detectors. Since vertex-tagging of 6-quarks 
in top events has been shown to be such a powerful tool, the capabilities of the upgraded 
DO detector for top physics will be considerably enhanced. The new CDF silicon vertex 
detector will also have a significant impact in the top physics program at CDF. With three 
dimensional information, the 6-tagging efficiency will be improved, and the instrumental 
background level will be reduced; furthermore, the new vertex detector will be long enough 
to cover the whole luminous region of the Tevatron, increasing the acceptance for vertex 
tags (recall that the geometrical coverage of the present CDF vertex detector is only about 
60%, see Section |V111.B.3| ). 



At the same time, the Fermilab accelerator complex will be upgraded with the construc- 
tion of a new high intensity 120 GeV proton accelerator (The Main Injector). The Main 
Injector will replace the Main Ring, whose aperture currently limits the luminosity, as the 
injector to both the Tevatron and the p source. In addition, a new 8 GeV permanent magnet 
ring (Recycler) has been proposed to achieve a more efficient accumulation of antiprotons 
by recycling (hence the name) the p from the previous store, store the anti-protons when 
the accumulator is full, and protect them from power glitches in the accelerator (Foster, 
1995; Jackson, 1995). The number of bunches will also be increased from 6 to 36, allowing 
higher luminosity without increasing the number of interactions per crossing. This increase 
shortens the bunch crossing interval from the current value of 3.5 /isec to 396 ns. This re- 
quires an upgraded trigger, data acquisition and front-end electronics systems, which must 
be pipelined to handle the increased rate. After completion of these upgrades in 1998-99, the 
Tevatron luminosity will be increased by one order of magnitude to 2 x 10 32 cm _2 s~ 1 , with 
further luminosity improvements likely to occur in the following years. The center-of-mass 
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energy of the Tevatron is also expected to increase from 1.8 TeV to 2 TeV, resulting in a 
increase of approximately 30% in the ti cross-section. 

This series of improvements in the Fermilab collider program will allow for much more 
detailed studies of the top quark than those that are possible with the low statistics data 
samples that are available now. The projected sizes of the ti samples in each experiment for 
an integrated luminosity of 1 fb _1 will be of order 100 events in the dilepton channel and 
of order 500 events in the lepton + 4 jets channel with one 6-tag, with half of these events 
having both 6-jets tagged (Amidei and Brock, 1996). 

In addition, towards the middle of the next decade, the Large Hadron Collider (LHC) 
at CERN will become operational. The LHC is a very high luminosity (> 10 34 cm -2 s _1 ) 
pp machine with y/s = 14 TeV. At this energy, a(pp — *> ti) ~ 700 pb, a factor of 100 higher 
than a(pp — > ti) at the Tevatron. With the high luminosity and the high cross-section, the 
LHC can be considered to be a top-factory. In the even more distant future, top physics 
will also be pursued at a very high energy electron collider (e.g. The Next Linear Collider, 
NLC), assuming that such a machine will be built. 



B. Improving the top mass measurement 

As was discussed in Section |H|, the top quark mass is a fundamental parameter of the 
Standard Model. Its value enters in the calculation of radiative corrections to a large number 
of electroweak observables. It is therefore very important to measure the top quark mass as 
accurately as possible to allow for precise tests of the Standard Model. 

The best value of M top as obtained from fits to the LEP and SLC measurements, as 
well as the measurements of the W mass from pp experiments and deep inelastic neutrino 
scattering, is M top = 179 ± 9ll\l GeV/c 2 where the second uncertainty comes from varying 
the Higgs mass between 60 and 1000 GeV/c 2 (see Section HI). This is in good agreement 
with the values reported by the CDF (M top = 176 ± 9 GeV/c 2 ) and DO collaborations 
(Mtop = 170 ± 18 GeV/c 2 ), see Section 0. 

The accuracies of the neutrino measurements and of the LEP and SLC measurements at 
the Z are not expected to dramatically improve in the coming years. On the other hand, the 
accuracy in the determination of the W mass will improve by almost one order of magnitude, 
from the Tevatron experiments as well as from e + e~ — > W + W~ at LEP200. Within the 
Standard Model, radiative corrections to the W propagator (see Fig. 11) result in definite 
predictions for the W mass as a function of the top and Higgs mass (see Fig. O). It is 
therefore very interesting to measure both the W mass and the top mass as precisely as 
possible. 

How accurately can the top mass be measured? Experience from CDF and DO indicates 
that the method of choice for measuring the top mass is to perform constrained fits on the 



6-tagged lepton + jets data sample, see Section IX. A . We can extrapolate the statistical 
accuracy of a future top mass measurement from the present CDF measurement. This mea- 
surement has a statistical uncertainty of 6 GeV/c 2 for an integrated luminosity of 110 pb -1 . 
After the first Main Injector run of the Tevatron, we can expect the integrated luminosity 
to be of order 1 fb _1 per experiment, with a 30% increase in top cross-section from running 
at the higher center of mass energy of 2 TeV, and a 40% increase in geometrical acceptance 
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from the new CDF vertex detector. Since the statistical uncertainty varies inversely as the 
square root of the number of events, we can expect a statistical uncertainty of order 1.5 
GeV/c 2 from CDF. A similar uncertainty can be expected from the upgraded DO detector. 
Systematic effects, which at present are at the level of 7 GeV/c 2 , will be the limiting factor 
in the precision of the top mass determination. 

The systematic uncertainties in the top mass measurements have been described in Sec- 
tion |LX..A.3|. The largest uncertainty is due to the understanding of the jet energy scale, 



as well as the related issue of additional jets from gluon radiation. At this time, it is not 
entirely clear what the ultimate precision will be. The dominant component of the energy 
scale uncertainty is related to the reliability of the extrapolation from jet energies to parton 
energies, and it is the understanding of the QCD process, rather than the understanding 
of the detector, which limits the measurement. Higher statistics tests of the understanding 
of the energy scale will be performed in 7 + jet and Z + jet events (see Section |IX.A.3| , 
Figs. |101| , and |102| ). The Monte Carlo modeling of gluon radiation will be more precisely 
checked and/or tuned by examining the energy flow within a jet. There will however remain 
systematic uncertainties related to the transfer of this calibration from the control samples 
to the hadronic environment in ti events. The ultimate size of these uncertainties is at the 
moment not well understood. 

The additional statistics that will become available will be important to reduce the 
systematic uncertainty. With enough statistics, the number of lepton + jets events with 
two 6-tagged jets will be sizable, and these events will provide a very important tool for 
understanding systematic uncertainties in the top mass measurement. In these events, there 
are no ambiguities in assigning the two jets to the hadronic W decay. The invariant mass of 
these two jets, which ideally should reconstruct to the W mass, provides an in-situ calibration 
of W — > jet-jet invariant mass reconstruction. First results from this kind of study are very 
promising, (see Fig. |106|) . Besides the energy scale issue, studies of doubly-6-tagged events 
will provide useful handles on other systematic effects that limit the precision of the top 
mass measurement. Since mis-tag backgrounds in this data sample are very small, events in 
the tails of the jet-jet invariant mass distribution will be ti events with one jet from gluon 
radiation and WQQ events. Therefore doubly-fe-tagged events will be useful in directly 
measuring these components of the data set. This will improve the understanding of the 
top resolution function, which also depends on the number of jets from gluon radiation in 
the sample. The doubly-6-tagged sample will also allow for a test of the modeling of the 
WQQ background component, which affects the top mass measurement since the top mass is 
extracted from a likelihood fit of the data to the sum of ti and W background. The number 
of jet-parton combinations for doubly-6-tagged events is only four, as opposed to twelve 
combinations that must be considered in the present CDF analysis which only demands 
one 6-tag, see Section |IX.A.1| . Since statistics are not expected to be the limiting factor, it 
is possible that a measurement of the top mass using these events will be ultimately more 
accurate than a measurement based on events with one 6-tag. Alternatively, these events can 
be used to study the effect of wrong combinations. Using these events, it will also become 
possible to measure the probability for the constrained fit to converge to a combination with 
a 6-tagged jet assigned to a light quark jet. 

Further improvements in the understanding of the systematic uncertainties will be pos- 
sible due to the high statistics ti samples that will be available. It will become possible to 
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better quantify the effects of gluon radiation by measuring the relative rates of tt — > dilepton 
+ 2 vs. 3 jets, and tt — ► 1+ 4 vs. 5 jets. The modeling of the W+ jets background, which at 
present is entirely based on the VECBOS Monte Carlo, will be tested using the large sample 
of Z+ jets. 

In summary, while the ultimate precision of the top quark mass measurement is not 
fully known at this time, it is likely that an accuracy of order of a few GeV/c 2 will be 
achievable. In conjunction with a W mass measurement with a precision of tens of MeV/c 2 , 
this measurement will be sensitive to physics beyond the Standard Model, and will provide 
useful information on the value of the Higgs mass. In Fig. |1(J8| we show what the M top vs. 
Myj/ measurements might look like by the year 2000. If the experimental point in the Myy 
vs. M top plane was to fall outside the region allowed by the Standard Model, this kind of 
measurement would provide indirect evidence for physics beyond the Standard Model. 




FIG. 108. The expected correlation between the masses of the top, the W and the Higgs. We also 
show, at an arbitrary point, results of measurements of the W mass to 50 MeV/c 2 and the top 
mass to 4 GeV/c 2 . 



C. Probing the Wtb vertex 

The structure of the Wtb vertex can be probed by studying top decay and/or production 
of single top quarks. Up to this point we have mostly discussed tt pair production. However, 
as mentioned in Section [TV], top quarks in pp collisions can also be produced singly, in the 
Drell-Yan process qq — > W* — ► tb, (see Fig. 17), and in the IV-gluon fusion process qg — ► tbq', 



139 



or qb — > tg', (see Fig. 18). The expected cross-sections for Drell-Yan and VT-gluon fusion 
single top production at the Tevatron are displayed in Fig. These cross-sections are 
smaller than the strong pp — > ti production cross-section, but they are sizable enough that 
single top production is expected to be observable at the Tevatron. Production of tW is 
also allowed, but its cross-section is much smaller. 

The potential for studying single top production has attracted a lot of attention for the 
future top physics program at the Tevatron. The cross-section is proportional to the top 
quark width, T(t — > Wb), and within the context of the Standard Model, proportional to 
the square of the CKM-matrix element |V#,|. Note that the top lifetime is too short to be 
directly measured, (see Fig. furthermore, the top width is very hard if not impossible 
to measure from the reconstruction of the Breit-Wigner, since the experimental resolution 
is one order of magnitude worse than the width itself (see Fig. |95[). Hence, only indirect 
measurements of T(t — > Wb) and |Vt&| can be performed. 

Assuming three generations and unitarity of the CKM matrix, the value of |V#,| is ex- 
pected to be very near unity, in the range 0.9988—0.9995 (Montanet et al, 1994). It is clearly 
very interesting to test this result. An additional potentially interesting measurement using 
the single top sample would be a comparison of the rates for pp —>■ tX and pp — > iX, which 
can be used to search for CP violating effects in the top sector. 

The possibility of extracting a single top-quark signal has been examined by many au- 
thors (Yuan, 1990; Cortese and Petronzio, 1991; Jikia and Slabospitsky, 1992; Ellis and 
Parke, 1992; Carlson and Yuan, 1993; Stelzer and Willenbrock, 1995; Amidei and Brock, 
1996). Because of the high multi-jet QCD background, only events with t — > Wb followed by 
W — > lu are useful. The signature then is a lepton, missing energy, and two 6-jets for Drell- 
Yan qq — > W* — > tb production, and one or two 6-jets + one light quark jet for VT-gluon 
fusion, (qg — > tbq' or qb — > tq', see Fig. 18). Just as in the ti search, the main background 
is from W+ jets production (see Section |V.E.1|) , and 6-tagging must be used to reduce 
this background to a manageable level. Rejection of the tagged WQQ background can be 
achieved by requiring the mass of the lepton, neutrino, and 6-jet to reconstruct to the known 
top mass. Pair production of ti also constitutes a significant background to observation of 
single top production. Since the jet multiplicity in ti events is higher, the optimal sample 
in which to isolate the single top signal seems to be that of events with one lepton + Er+ 
two and only two jets. 

A study of the expected signal and background for single top production at the Tevatron 
shows that a signal-to-background of order 1 to 2 can be achieved. The number of signal 
events for |V#,| rs 1 would be of order 120 per fb _1 , resulting in a statistical uncertainty 
in the cross-section measurement of order 17% in a 1 fb _1 data set (Amidei and Brock, 
1996). An additional study, optimized for the detection of the Drell-Yan qq — > W* — ► tb 
process, suggests that with the requirement of a double 6-tag the signal can be isolated, (see 



Fig. |109| Stelzer and Willenbrock, 1995). Indications from this study are that, assuming 



that | Vtb | is indeed close to unity, a 3 fb _1 exposure would yield a 20% measurement of the 
cross-section, and therefore a 10% measurement of If on the other hand no signal is 

seen in 3 fb _1 , then \V t b\ < 0.60 at the 95% confidence level. The authors of this study also 
suggest that the Drell-Yan process may be more useful in extracting |Vtb| than the H^-gluon 
fusion process. This is because expectations for the cross-section for pp — > W* — > tb are well 
understood, and can be also normalized to the observed rate of pp — > W* — > lu. In contrast, 
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the calculation of VF-gluon fusion suffers from uncertainties in the higher order corrections 
as well as in the input gluon density. 



Tevatron : 
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FIG. 109. Expected observed cross-section for pp — ► W* — > tb , t — > Wb, W — > Iv as a function of 
the invariant mass of the Wb. Also shown are expectations for the most important backgrounds 
(tbj denotes VF-gluon fusion). This analysis is based on the W+ 2 jets sample. From Stelzer and 
Willenbrock, 1995. 

The CKM matrix element \V t b\ can also be measured from the ratio of branching ratios 
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The 6-tag can differentiate between t — > Wb and t — > VTrf or W^s. Therefore, this ratio can 
be measured by comparing the number of single- and double-6-tagged tt — > lepton + jets 
events and by measuring the tagging rate in tt — > dilepton events. A preliminary analysis 
by the CDF collaboration finds (Yao, 1995): 
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and sets a (not-yet-very-interesting) limit \Vtb\ > 0.022 at the 95% confidence level. The 
statistical sensitivity for this branching ratio measurement is projected to be 3% for a 1 fb _1 
exposure (Amidei and Brock, 1996). 

Top decays also provide a unique opportunity to test the the structure of the charged 
weak current at the t — > Wb vertex. Because the top quark is so heavy, it is possible that 
new physics may manifest itself at the Wtb vertex. The most general form of the Wtb 
interaction is (Kane, Ladinsky, and Yuan, 1992) : 
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where P± = |(1 ± 75) and ia^ = — |[7 M , 7"]. 

The quantities and f R parametrize the strength of the left-handed and right-handed 
weak charged current. The / 2 's can be interpreted as giving rise to an anomalous weak 
magnetic moment. In the Standard Model at tree level f[ = 1 and ff = f 2 L — f 2 — 0. There 
is obviously no direct experimental information on these form factors, although consistency 
with the measured branching ratio for b — > 57 constrains f R to be at most a few percent 
(Fujikawa and Yamada, 1994). 

The polarization of the W in the t — > Wb decay probes the values of the form factors. 
Denoting the left-handed, right-handed, and longitudinal polarization states of the W by A_, 
A + , and Ao respectively, the expected relative polarizations of the W boson after averaging 
over the top and bottom polarization states, are (Kane, Ladinsky, and Yuan, 1992) : 

A_ = \tf + Pf?\ 2 
A+ = \f? + Pfi\ 2 

Ao = llff + Pttl 2 + \\f 2 L + PA R \ 2 



with (3 = -^ £ . Therefore, the Standard Model predicts A+ = 0, and the fraction of longitu- 
dinally polarized VF-bosons in top decays is ^(3 2 . These polarizations can be measured from 
the angle of emission of the lepton in VT^-decays. Studies indicate that the statistical accu- 
racies in the measurements of Ao and A + in an exposure of 1 fb _1 at the Tevatron should be 
5% and 2% respectively (Amidei and Brock, 1996). These form factors can also be probed 
by measuring the single top production cross-section (Carlson, Malkawi, and Yuan, 1994; 
Malkawi and Yuan, 1994). 

The longitudinal polarization state of the W is directly connected with the breaking of 
electroweak symmetry, since it arises from the Goldstone boson degree of freedom. This kind 
of study therefore provides a rather unique and particularly interesting test of the Standard 
Model. 
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D. Further tests of the Standard Model and searches for new physics in the top quark sector 



The pp — > ti cross-section can be calculated in QCD. Its measurement tests the predic- 
tive power of QCD and is sensitive to new physics. Considerable theoretical interest in the 
subject was triggered by the initial measurement of the ti cross-section by the CDF collab- 
oration (F. Abe et al, 1994a). The measured value was higher than expected, although still 
consistent with the QCD calculation, within the large experimental uncertainties. The more 
recent higher statistics measurement by both CDF and DO are in better agreement with the 
calculation (see Fig. |92|). 

The precision of the ti cross-section measurement depends on the accuracy of the lu- 
minosity normalization (3.5%), the background estimate, and the acceptance calculation. 
The uncertainty in the acceptance calculation is mostly due to the uncertainty in the 6-tag 
efficiency and the uncertainties in the modeling of ti production (e.g. the effects of gluon 
radiation). The high statistics ti data samples that will be collected at the Tevatron will pro- 
vide several handles to reduce this uncertainty. It is not clear what the ultimate systematic 
uncertainty on the ti cross-section will be; our guess is that a precision of 10%, comparable 
to the uncertainty in the QCD calculation of ti production, should be achievable. 
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FIG. 110. Expected invariant mass of the tt pair : (a) including the contribution of a color octet 
vector meson in the top color model (from Hill and Parke, 1994). The different curves show the 
expectation as the vector meson mass varies from 600 GeV/c 2 to infinity, (b) including the effect 
of a 475 GeV/c 2 color-octet technipion in multiscale walking technicolor (from Lane, 1995). The 
dotted line is the QCD prediction; the dashed line is the technipion contribution; and the solid line 
is the sum of the two. 



Large enhancements to the ti cross-section and, more dramatically, resonances in the ti 
invariant mass spectrum are expected in a number of models, (see Fig. |110| ). These could 
be due to color octet vector mesons (Hill and Parke, 1994) in models where the electroweak 
symmetry breaking is realized via top condensation (Hill, 1991; Martin, 1992a and 1992b), 
or to technipions (Eichten and Lane, 1994) in multiscale models of walking technicolor (Lane 
and Eichten, 1989; Lane and Ramana, 1991). 

Measurements of the ti invariant mass distributions are possible in the lepton + jets 
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channel. Preliminary CDF (Yao, 1995) and DO (Narain, 1996) analyses are consistent with 
QCD expectations, (see Fig. |1 1 1|) . With the high statistics that will become available, these 
models will be critically tested in the future. 
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FIG. 111. Reconstructed invariant mass of the tt pair compared to the background and tt + 
background expectations. From the CDF collaboration, Yao, 1995. The integrated luminosity is 
« 100 pb" 1 . 



Chromoelectric and chromomagnetic dipole moments of the top quark affect the structure 
of the tt-gluon vertex, and hence affect the tt production cross-section and tt transverse mo- 
mentum (Atwood, Kagan, and Rizzo, 1994; Rizzo, 1994; Cheung, 1995; Haberl, Nachtmann, 
and Wilch, 1995). A top-quark chromomagnetic dipole moment can occur in composite and 
technicolor models, with magnitude of order M^ op /A 2 , where A is the characteristic scale for 
new physics. As can be seen from Fig. [TTI3, the present measurement of the tt cross-section is 
already accurate enough to probe the scale A m 200 GeV. A chromoelectric dipole moment 
would be CP-violating, and could arise from large couplings between top quarks and Higgs 
bosons in the multi-Higgs doublet model (Atwood, Aeppli, and Soni, 1992; Brandenburg 
and Ma, 1993; Cheung, 1995; Haberl, Nachtman, and Wilch, 1995). 

CP violation effects in tt production would manifest themselves in different polarizations 
for the t and the i (Kane, Ladinsky, and Yuan, 1992; Schmidt and Peskin, 1992, Kao, 
Ladinsky, and Yuan, 1994). Because the top quark lifetime is short, the top decays before 
hadronization and polarization information is preserved in the t — > Wb decay. Assuming 
Standard Model Wtb couplings, the polarization of the top quark is then analyzed by the 
polarization of the W from the top decay. Examples of CP violating observables that can be 
studied are differences in rates between ttfi, and t R iji, or t L t R and t R ti, where the subscripts 
L and R denote left- and right-handed polarizations respectively. 

In the single top production process, top quarks are almost 100% longitudinally polarized, 
since they are produced through the weak interaction (Carlson and Yuan, 1993). Searches 
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FIG. 112. Contours of tt cross-section (M top = 176 GeV/c 2 ) as a function the chromoelectric (k) 
and chromomagnetic (k) dipole moments of the top quark. 

for CP violation in the top quark decay can then also be carried out in the t — > W + b — > l + vb 
decay by studying the quantity o • (P b x P ; ), where a is the top polarization vector (Kane, 
Ladinsky, and Yuan, 1992; Grzadowski and Gunion, 1992). 

Precise measurements of polarization effects in the top sector will however be rather dif- 
ficult. These measurements require very high statistics and good control of the systematics. 
Most likely, only very large asymmetries will be accessible experimentally. 

Physics beyond the Standard Model can also give rise to exotic decays of the top quark. 
The measurement of the ratio of branching ratios 

Br{ti — > / + jets) 
~ Br(tt -> 11+ jets) 

is quite generally sensitive to decays different from t — > WX, provided that the new top 
quark decay mode includes jets. More model-dependent searches for new top quark decays 
can also be carried out. 

One example of a non-Standard Model decay of the top quark is the decay t — > H + b, 
which can occur in models with two Higgs doublets, see Sec [VI. Gj . The H + would then 
decay into the heaviest fermion pairs, cs or rv. The branching ratios depend on the ratio of 
vacuum expectation values for the two doublets, tan/3, (see Fig. |113|) . If tan/3 is large, the 
signature for this decay mode would be an excess of lepton + hadronic r + 6-tag events; if 
tan/3 is small, one could search for the H + — > cs peak in the invariant mass distribution of 
lepton + jets events. Other possible exotic decays of the top quark which can be searched for 
include flavor changing neutral current decays such as t — > Zc (Han, Peccei, and Zhuang, 
1995) and t — > 7c, and decays into a supersymmetric top quark (stop) and a neutralino 
(Mrenna and Yuan, 1995). Some preliminary limits on these FCNC have been presented by 
CDF (LeCompte, 1995). 
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FIG. 113. (a) Expected branching ratios of t — * H + b in the two Higgs model for Mt op = 180 
GeV/c 2 as a function of log(tan/3). Solid line: M H = 70 GeV/c 2 ; dashed line: M H = 150 GeV/c 2 . 
(b) Expected branching ratio for H + — > rv as a function of log(tan/3). 



XI. CONCLUSION 

The evidence for the existence of the top quark from the CDF and DO collaborations 
at the Tevatron is persuasive. The mass of the top quark is M top = 175 ± 8 GeV/c 2 , in 
agreement with expectations from precision electroweak measurements. This mass is a factor 
of forty higher than the mass of the second heaviest fundamental fermion, and is of the same 
order of magnitude as the scale for electroweak symmetry breaking. The high mass of the top 



quark is a somewhat of a surprise; in Fig. |114| we show the evolution of the top mass limits 



and measurements since the discovery of the companion 6-quark. It is interesting to notice 
that theoretical arguments based on local supersymmetry from the early eighties, when the 
experimental lower limit on the top quark mass was only approximately 20 GeV/c 2 , favored 
a rather high top mass, see Section J|. It is however still far from clear whether the high 
value of the top mass is an accident, or a consequence of physics at a higher mass scale. 

The properties of the top quark will be studied much more precisely at the upgraded 
Tevatron starting in 1999. There is a possibility that effects beyond the Standard Model will 
manifest themselves in the top sector. If that is the case, the CDF and DO collaborations are 
well positioned to observe them. In the more distant future, the LHC pp collider at CERN, 
with the higher energy and luminosity, will serve as a tt factory. 
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FIG. 114. A brief history of the search for the top quark. Here we show the limits from searches in 
e + e~ (Section VI. A) and pp (Section VLB). The square data points are predictions from neutral 
current studies, mostly at LEP and SLC (Section III). We also show the limit from neutrino 
experiments that were obtained prior to the commissioning of LEP and SLC (Amaldi et al, 1987). 
The circular data points show the recent measurements of the top mass (Section IX. A), as well as 
the UA1 1984 result M top = 40 ± 10 GeV/c 2 , which turned out to be wrong (Section VLB). Just 
for fun, we also show the upper limit on the top mass from early measurements of the W leptonic 
branching ratio at the SppS (Section VI. C). More recent measurements of this quantity at both 
the Tevatron and the SppS are consistent with a high top quark mass. 
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